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CHAPTER I 


INTRODUCTION 

Biodiversity  is  declining  at  an  alarming  rate.  Although  it  is  difficult  to  quantify 
the  historical  rate  of  species  loss,  scientists  generally  agree  that  current  extinction  rates 
far  exceed  those  found  in  the  geologic  record  (Probst  and  Crow  1991). 

The  impact  of  human  alteration  of  the  natural  landscape  on  species  habitat  is  a 
key  destructive  force  leading  to  this  decline  of  biodiversity.  The  adverse  impact  of 
human  activities  is  not  limited  to  simply  the  destruction  of  a  species  habitat. 
Fragmentation  of  the  habitat  that  remains  can  also  have  profound  effects  depending  on 
when  and  how  the  habitat  is  fragmented.  Species  that  are  found  in  the  remaining  habitat 
fragments  often  become  isolated  from  other  populations  and  suffer  from  adverse  effects 
of  changes  in  the  surrounding  landscape.  Hence,  habitat  fragmentation  is  a  major 
concern  relating  to  biological  diversity  (Harris  1984,  Saunders  et  al.  1991)  and  is  an 
example  of  how  changes  in  the  specific  spatial  parameters  of  a  habitat  within  a  landscape 
can  be  important  to  species  survival.  The  effect  of  habitat  fragmentation  on  biological 
diversity  is  an  important  consideration  in  the  conservation  of  the  Earth's  resources 
because  of  detrimental  effects  on  biodiversity  and  the  distribution  and  abundance  of 
individual  species.  However,  the  degree  to  which  a  given  species  is  affected  by 
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habitat  fragmentation  is  dependent  on  a  complex  interaction  of  the  habitat  requirements 
of  the  species  and  the  shape,  size,  and  makeup  of  the  fragmented  habitat. 

Conservation  of  the  biological  diversity  of  a  landscape  would  be  facilitated  if 
there  was  a  way  to  determine  the  impact  of  habitat  changes  on  species  of  interest. 
Presently  there  is  no  good  way  to  assess  or  quantify  the  impact  of  landscape  changes  to 
habitat  suitability,  even  for  species  of  conservation  priority.  The  research  presented  here 
attempts  to  provide  this  capability  for  conservation  priority  Neotropical  and  resident 
songbirds  found  in  fragmented  landscapes  of  the  southeastern  United  States,  utilizing 
existing  data  sources,  geographical  information  system  and  spatial  analysis  software. 
This  research  is  now  feasible  due  to  recent  innovations  in  computer  and  remote  sensing 
technologies  that  have  reduced  computation  expenses,  provided  new  sources  of  spatially 
extensive  data,  and  provided  the  impetus  for  the  development  of  new  spatial  analysis 
software. 

Existing  land  use  and  land  cover  spatial  data  are  widely  available  and  can 
potentially  be  used  to  evaluate  the  impact  of  changes  in  the  landscape  on  species  of 
conservation  priority.  The  scale  of  the  data  varies  and  selection  of  a  data  source  is 
complicated  by  the  fact  that  there  is  no  single  best  scale  for  determining  landscape 
changes  to  habitat  suitability.  The  scale  depends  on  the  question  asked,  the  species  and 
habitats  involved,  and  the  processes  believed  to  be  important  (Wiens  1989). 

Two  existing  and  publicly  available  datasets  provide  an  opportunity  for 
development  of  methods  in  this  area.  They  are  the  USGS  Land  Use  Land  Cover  (LULC) 
and  the  Breeding  Bird  Survey  (BBS)  data  sets.  The  USGS  LULC  dataset  is  a  kilometer- 
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resolution  remotely  sensed  dataset  that  utilizes  a  uniform  classification  scheme  for  the 
entire  country.  The  BBS  is  a  broad  scale  bird  population  survey  dataset  that  covers  the 
United  States  and  Canada.  The  BBS  is  the  only  survey  of  this  scope  and  apparently  has 
never  been  used  in  quantifying  landscape  habitat  suitability  for  bird  species. 

Objective 

The  objective  of  this  research  was  to  use  the  BBS  and  USGS  LULC  data  to 
determine  if  kilometer-resolution  horizontal  spatial  pattern  metrics  are  suitable  indicators 
of  habitat  suitability  for  Neotropical  and  resident  songbirds.  The  study  area  covers  a 
block  of  the  Southeastern  USA,  which  included  90%,  30%,  and  30%  of  Alabama, 
Georgia,  and  Florida  respectively  (Figure  1).  It  includes  15  bird  species  listed  by 
Southeast  Partners  in  Flight  as  a  conservation  priority  (Table  1)  (Hunter  1998).  This 
study  focused  on  determining  the  utility  of  existing  data  in  predicting  bird  abundance  and 
evaluating  the  sensitivity  of  predictive  models  to  varied  size  of  landscape  units  analyzed. 


Figure  1 .  Study  area  boundary,  dates  of  photography  from  which  USGS  land 
use/land  cover  data  were  interpreted,  and  breeding  bird  survey  routes. 
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Table  1 .  Conservation  priority  bird  species  under  study  in  this  research,  the  American 
Ornithologist  Union  (AOU)  number  and  the  species  preferred  habitat. 


Species  Name 

AOU  No. 

Habitat 

Prairie  Warbler 

6730 

Abandoned  fields,  cut  over  burned  over  woods,  woodland  margins 

Northern  Bobwhite 

2890 

Abandoned  fields,  brushy  areas,  hedgerows,  thickets,  woodland 
margins,  open  woods 

Field  Sparrow 

5630 

Favor  scattered  saplings  or  shrubs  in  weedy  habitats,  overgrown 
fields,  woods  margins,  hedgerows,  and  thickets 

Loggerhead  Shrike 

6220 

Open  country,  open  fields,  pastures,  cultivated  fields  where  there 
are  scattered  trees  for  nesting  and  telephone  wires  or  fences  for 
perching 

Eastern  Kingbird 

4440 

Open  country,  prefer  areas  with  scattered  trees,  fences,  telephone 
wires,  and  fields 

White-eyed  Vireo 

6310 

Dense  thickets,  especially  where  moist,  common  habitats  are 
stream  side  shrubbery,  swamp  borders  and  openings,  willow 
thickets,  and  damp  tangles 

Eastern  Wood-pewee 

4610 

Open  to  medium-growth  forests  and  woodlots,  favor  neither  pines 
nor  hardwoods 

Yellow-throated  Vireo 

6280 

Wide  variety  of  woodlands,  favor  mature  deciduous  trees  in  fairly 
open  setting,  especially  where  moist,  avoid  pure  coniferous 
forests,  open  hardwoods,  woodland  borders,  extensive  forest 

Hooded  Warbler 

6840 

Primarily  in  deciduous  forests  but  also  in  mixed  forests,  favor 
moist  forests  with  a  fairly  dense  understory,  such  as  bottomlands, 
rich  woods,  and  ravines.  Sometimes  in  the  deciduous  understory 
of  mature  pine  forests,  interior  species 

Wood  Thrush 

7550 

Deciduous  or  mixed  forests  with  a  fairly  well-developed 
deciduous  understory,  especially  where  moist,  interior  species 

Yellow-billed  Cuckoo 

3870 

Deciduous  forests,  bottomland  woods,  woodland  thickets,  other 
hardwood  forest,  generally  avoid  coniferous  woods,  extensive 
forest 

Kentucky  Warbler 

6770 

Rich,  moist  deciduous  forests,  bottomland  forests,,  seldom  in 
conifers,  area  sensitive,  interior  specie 

Brown-headed  Nuthatch 

7290 

Pines,  prefer  mature,  open  pinewoods 

Prothonotary  Warbler 

6370 

Bottomlands  forests,  almost  always  near  standing  water,  prefer 
swamps  that  are  somewhat  open  with  scattered  dead  stumps, 
extensive  willow  thickets  near  lakes  or  ponds,  interior  species 

Brown-headed  Cowbird 

4950 

Open  woods,  margins,  thickets,  farmyards  and  residential  areas 

CHAPTER  II 


LITERATURE  REVIEW 

An  extensive  literature  review  was  conducted  to  determine  the  current  state-of-art 
in  the  assessment  of  the  impact  of  landscape  changes  on  habitat  suitability.  The  review 
was  useful  in  gamering  information  on  the  promising  techniques  utilized  in  past  studies 
as  well  as  valuable  background  information  that  provides  context  for  the  current  effort. 
The  discussion  of  the  literature  is  presented  in  three  parts.  First,  information  is  presented 
that  is  applicable  to  studies  relating  species  to  landscape  structure.  Second,  a 
comprehensive  review  of  spatial  pattern  metrics  found  in  the  literature  is  presented, 
including  information  on  the  specific  metrics  used  in  this  research.  Third,  a  detailed 
review  of  those  landscape  level  studies  of  associations  of  horizontal  spatial  pattern  to  bird 
abundance  that  are  specifically  pertinent  to  this  study  is  provided.  The  definition  of  some 
of  the  more  specialized  terms  used  in  this  discussion  can  be  found  in  Table  2. 

Relating  Species  to  Landscape  Structure 

Habitat  patches  are  distributed  within  landscapes  and  their  spatial  patterns  may 
exert  influence  on  the  abundance,  distribution,  and  dynamics  of  vertebrate  populations 
inhabiting  those  landscapes  (Wiens  1976,  1989).  The  spatial  configuration  of  a 


Table  2.  Specialized  Terms  and  Definitions 
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Fragmentation  -  an  external  disturbance  that  alters  the  large  patch,  thus  creating  isolated 
or  tenuously  connected  patches  of  the  original  habitat,  interspersed  with  an  extensive 
mosaic  of  other  habitat  types  (Wiens  1989).  There  are  two  components  associated  with 
fragmentation:  reduction  in  the  total  area  of  a  habitat  type,  and  apportionment  of  the 
remaining  habitat  into  smaller,  more  isolated  patches  (Saimders  et  al.  1991,  Harris  1984). 
Landscape  -  a  heterogeneous  land  area  composed  of  an  interacting  mosaic  of  habitat 
patches  (Forman  and  Godron  1986). 

Landscape  ecology  -  the  study  of  the  structure,  function  and  change  in  the  landscape 
composed  of  interacting  ecosystems.  It  provides  the  basis  for  studies  directed  at 
investigating  associations  dealing  with  landscape  structure  (Forman  and  Godron  1986). 
Landscape  scale  or  Broad  scale  -  refers  to  a  large  area  such  as  a  landscape.  Landscape 
varies  in  size  down  to  a  few  kilometers  in  diameter  (Forman  and  Godron  1986). 
Landscape  structure  -  composition  and  spatial  configuration  of  the  patches  within  a 
landscape.  Landscape  composition  considers  the  presence  and  amount  of  each  patch 
type  within  the  landscape  but  does  not  address  the  location  of  patches  within  the 
landscape.  Spatial  configuration  refers  to  the  physical  distribution  or  spatial  character  of 
patches  within  the  landscape. 

Local  scale  or  Fine  scale  -  refers  to  a  small  area  such  as  a  patch.  It  includes  between- 
plot  and  within-plot  studies.  Localized  areas  of  a  few  meters  or  hundreds  of  meters 
across  are  at  a  finer  scale  than  a  landscape  (Forman  and  Godron  1986). 

Patch  -  a  nonlinear  surface  area  differing  in  appearance  from  its  surroundings. 

Spatial  scale  -  refers  to  the  spatial  size  or  the  ratio  of  length  on  a  map  to  true  length. 
Spatial  scale  of  datasets  involve  both  grain  and  extent.  Grain  refers  to  the  resolution  of 
the  data,  that  is  the  area  represented  by  each  data  unit.  Extent  refers  to  the  overall  size 
of  the  study  area  (Turner  et  al.  1989). 
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species  habitat  affects  populations  by  influencing  the  patterns  of  movement  of 
individuals,  interactions  among  individuals,  and  exposure  to  factors  associated  with  the 
edge  of  habitats  (Wiens  1976,  1989).  Landscape  ecology  attempts  to  provide  quantitative 
measures  of  landscape  structure.  Quantitative  measurements  of  landscape  structure  are  a 
necessity  in  the  empirical  investigation  of  landscape  structure  and  bird  species 
associations. 

The  hypothesis  that  landscape  structure  plays  an  important  role  in  the  habitat 
suitability  for  populations  has  come,  in  part,  from  field  studies  on  forest  fragmentation 
(Saunders  et  al.  1991).  Fragmentation  of  habitats  has  received  much  attention,  especially 
in  agricultural  areas  where  natural  vegetation  is  broken  into  small,  isolated  patches  as 
land  is  converted  for  human  use.  The  majority  of  work  by  ecologists  has  focused  on  the 
fragmentation  of  forests  (e.g.  Harris  1984),  possibly  because  forest  habitat  loss  is  so 
visible  and  the  recovery  period  so  long  (Wiens  1989). 

How  different  bird  species  respond  to  habitat  fragmentation  is  determined  by  then- 
individual  spatial  requirements  and  affinity  for  edge  versus  interior  locations  in  habitat 
patches.  A  difference  in  individual  area  requirements  is  the  most  obvious  factor  in  their 
response.  Minimal-area  requirements  among  species  have  been  documented  in  numerous 
studies  (e.g.  Lynch  1987,  Robbins  et  al.  1989).  In  most  cases,  the  determination  of 
minimal  area  requirements  is  based  on  the  presence  or  absence  of  a  species  in  habitat 
fragments  of  different  sizes.  A  species  response  to  fragmentation  is  also  based  on  its 
affinity  for  edge  versus  interior  habitats.  As  the  size  of  a  habitat  patch  is  reduced  by 
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fragmentation,  the  proportion  of  the  fragment  that  adjoins  other  habitat  types  (edge) 
increases  (Forman  and  Godron  1986).  Studies  based  on  the  eastern  deciduous  forests  of 
North  America  show  that  the  abundance  of  vertebrate  species  associated  with  forest 
interiors  usually  declines,  while  the  abundance  of  vertebrate  species  specializing  on  forest 
edges  increases  in  response  to  forest  fragmentation  caused  by  agricultural  development 
and  urbanization  (Whitcomb  et  al.  1981,  Robbins  et  al.  1989,  Terborgh  1989).  Changes 
in  vegetation,  food  resources,  predation,  brood  parasitism,  and  competition  have  been 
noted  as  causes  of  the  observed  vertebrate  community  changes  (Kroodsma  1982, 
Brittingham  and  Temple  1983).  In  North  America,  neotropical  migratory  songbirds  are 
thought  to  be  particularly  vulnerable  to  habitat  fragmentation,  since  they  breed  primarily 
in  extensive  stands  of  mature,  floristically  diverse  forest  (Rosenberg  et  al.  1999). 
Extensive  stands  of  mature  forest  protect  from  the  proximity  to  a  habitat  edge,  especially 
for  the  single-brooded  low-nesting  neotropical  migrants  (Whitcomb  et  al.  1981).  Since 
nest  parasitism  by  Brown-headed  Cowbirds  decreases  with  distance  away  from  the  forest 
edge,  and  since  Brown-headed  Cowbird  breeding  is  concentrated  around  the  period  when 
neotropical  migrants  are  nesting,  these  species  are  especially  vulnerable  to  nesting  failure 
as  a  result  of  Brown-headed  Cowbird  activity. 

Most  studies  on  fragmentation  have  employed  a  patch-centered  sampling  scheme 
in  which  independent  forest  patches,  not  landscapes,  were  sampled  (e.g.  Rosenberg  and 
Raphael  1986,  Lehmkuhl  and  Ruggiero  1991).  Based  on  the  relationships  developed 
between  species  abundance  or  richness  and  a  number  of  patch  characteristics,  such  as 
patch  size,  inferences  were  made  about  how  landscape  structure  affects  wildlife 
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populations.  It  is  unclear  however,  whether  relationships  derived  at  the  patch  level  can 
be  extrapolated  to  the  landscape  level. 

Focusing  exclusively  on  fragmentation  of  habitat  patches,  however,  neglects  the 
fact  that  it  is  often  the  structure  of  an  entire  landscape  mosaic  rather  than  the  size  or 
shape  of  individual  patches  that  is  important  to  birds  (Wiens  1989).  The  likelihood  that 
dispersal  can  occur  between  fragments  is  influenced  by  the  configuration  of  the 
fragments  and  the  landscape  mosaic  in  which  they  are  located.  Fragments  of  habitat  are 
not  surrounded  by  totally  inhospitable  environments,  as  are  oceanic  islands.  Even  if 
patches  are  completely  isolated  from  areas  of  similar  habitat,  the  dynamics  of  their 
populations  may  be  influenced  by  features  of  the  surrounding  habitats  or  distances  to 
other  patches  of  the  same  habitat  (Wiens  1989).  Thus,  the  spatial  juxtaposition  of 
habitats  is  an  important  landscape  consideration  and  can  best  be  examined  by  landscape- 
level  studies  which  have  not  been  adequately  addressed  in  previous  studies  of  habitat 
suitability. 

The  more  recent  studies  found  in  the  literature  show  a  trend  towards  broadening 
the  spatial  scope  of  study.  This  has  entailed  going  from  an  emphasis  on  local-scale 
processes  to  the  inclusion  of  processes  influencing  plant  and  animal  populations  that 
occur  at  a  variety  of  spatial  scales,  including  the  landscape  scale.  One  advantage  of 
landscape-level  studies  of  species-habitat  associations  is  the  expanded  spatial  dimension 
and  the  connections  to  other  landscapes.  Studies  at  these  broader  scales  can  be  used  to 
understand  the  spatial  connections  present,  and  to  more  effectively  manage  the  broader 
ecosystem.  It  should  be  emphasized  the  limited  geographic  scope  of  most  forest 
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fragmentation  studies  stems  from  the  inherent  difficulty  in  surveying  large  land  areas 
in  a  limited  time  frame  (Rosenberg  et  al.  1999). 

Spatial  Pattern  Metrics 

A  landscape  is  distinguished  by  the  spatial  relationships  between  its  component 
parts.  Landscape  structure  is  determined  by  both  the  composition  and  configuration  of  a 
given  landscape.  Composition  considers  the  presence  and  amount  of  each  patch  type 
within  the  landscape  but  does  not  address  the  location  of  patches  within  the  landscape. 
Configuration  refers  to  the  physical  distribution  or  spatial  character  of  patches  within  the 
landscape  (McGarigal  and  Marks  1994).  Spatial  pattern  metrics  are  used  to  quantify  the 
horizontal  landscape  structure.  They  can  be  grouped  as:  1)  area  metrics,  2)  core  area 
metrics,  3)  patch  density  and  size  metrics,  4)  edge  metrics,  5)  shape  metrics,  6)  diversity 
metrics,  and  7)  interspersion/juxtaposition  metrics.  Since  the  metrics  typically  are  highly 
inter-correlated,  several  studies  have  used  principal  component  analysis  to  develop 
independent  factors  (McGarigal  and  McComb  1995,  Pearson  1993,  Rosenberg  et  al. 
1999).  The  specific  metrics  employed  in  this  study  were  selected  based  on  their 
successful  use  in  past  studies  and  documentation  in  the  literature.  A  discussion  of  each 
metric  group  is  presented  below.  The  mathematical  equations  used  to  calculate  the 
specific  metrics  used  in  this  study  are  also  provided.  Table  3  provides  definitions  of 
symbols  and  subscripts  used  in  the  equations.  Three  sample  sites  used  in  this  study  were 
selected  to  provide  graphic  examples  of  the  metrics.  The  three  sample  sites  were  selected 


Table  3.  List  of  Symbols  and  Subscripts 


disjunct  core  areas 


=  1,  ...  ,n  patches,  within  specified  neighborhood 


=  Area  (m2)  of  patch  ijs  within  specified  neighborhood  (m)  of  patchj. 


=  Core  area  (m  )  of  patch  ij  based  on  the  specified  buffer  width  (m 


=  Perimeter  (m)  of  patch  i 


=  Length  (m)  of  edge  of  patch  ij  adjacent  to  patch  type  (class)  k.  _ 


=  Total  length  (m)  of  edge  in  landscape;  includes  landscape  boundary  and  background 
edge  segments  if  the  user  decides  to  treat  boundary  as  edge;  otherwise,  only  boundary 
segments  representing  true  edge  are  included, _ 


=  Total  length  (m)  of  edge  in  landscape;  includes  entire  landscape  boundary  and 
background  edge  segments  regardless  of  whether  the  represent  true  edge, _ 


=  Total  length  (m)  of  edge  in  landscape  between  patch  types  (classes)  i  and  k;  includes 

landscape  boundary  segments  representing  true  edge  only  involving  patch  type  i. _ 

=  Total  length  (m)  of  edge  in  landscape  between  patch  types  (classes)  i  and  k;  includes  all 
landscape  boundary  and  background  edge  segments  involving  patch  type  i,  regardless  of 
whether  they  represent  true  edge. 


=  Total  length  (m)  of  edge  in  landscape  between  patch  types  (classes)  i  and  k;  includes 
entire  landscape  boundary  and  background  edge  segments,  regardless  of  whether  they 
represent  true  edge.  _ 


=  Dissimilarity  (edge  contrast  weight)  between  patch  types  i  and  k _ 


=  Total  number  of  patches  in  the  landscape,  excluding  any  background  patches. _ 


=  Total  number  of  patches  in  the  landscape  that  have  nearest  neighbors.  _ 


=  Number  of  patches  in  the  landscape  of  patch  type  (class)  i. _ 


=  Number  of  patches  in  the  landscape  of  patch  type  (class)  i  that  have  nearest  neighbors. 


=  Number  of  disjunct  core  areas  in  patch  ij  based  on  specified  buffer  width. _ 


=  Number  of  patch  types  (classes)  present  in  the  landscape,  excluding  the  landscape 

border  if  present.  _ _ _ 

=  Number  of  patch  types  (classes)  present  in  the  landscape,  including  the  landscape 

border  if  present. _ _ _ 

=  Number  of  patch  types  (classes)  present  in  the  landscape.  _ 

=  Distance  (m)  from  patch  ij  to  nearest  neighboring  patch  of  the  same  type  (class),  based 

on  edge-to-edge  distance. _ _ _ 

=  Number  of  adjacencies  (joins)  between  pixels  of  patch  types  (classes)  I  and  k. _ 

=  Proportion  of  the  landscape  occupied  by  patch  type  (class)  i. _ _ 
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as  they  have  varying  degrees  of  forest  fragmentation  and  illustrate  a  range  in  value  for 
the  metrics. 

1)  Composition 

Area  metrics  quantify  landscape  composition.  The  area  of  the  patches  comprising 
a  landscape  is  perhaps  the  single  most  important  piece  of  information  contained  in  the 
landscape  (McGarigal  and  Marks  1994).  Patch  area  has  ecological  utility,  as  shown  by 
the  considerable  evidence  that  bird  species  occurrence,  abundance  and  richness  are 
strongly  correlated  with  patch  size  (Robbins  et  al.  1989).  Knowing  how  much  of  the 
landscape  is  comprised  of  a  particular  patch  type  is  an  important  measure  in  ecological 
applications.  For  example,  in  forest  fragmentation  there  is  often  a  quantitative  loss  of 
habitat.  Thus,  in  studies  of  forest  fragmentation,  it  is  important  to  know  how  much 
habitat  exists  within  the  landscape.  Also,  even  though  many  vertebrate  species  that 
specialize  on  a  habitat  have  minimum  area  requirements  (e.g.  Robbins  et  al.  1989),  not  all 
species  require  the  habitat  to  be  in  one  contiguous  patch.  The  Northern  Spotted  Owls 
have  minimum  area  requirements  for  late-seral  forest  that  varies  geographically,  but 
individuals  use  late-seral  forest  that  may  be  distributed  among  many  patches  (Forsman  et 
al.  1984).  Thus,  late-seral  forest  area  might  be  a  good  index  of  habitat  suitability  within 
landscapes  the  size  of  Northern  Spotted  Owl  home  ranges  (Lehmkuhl  and  Raphael  1993). 


It  is  often  desirable  to  quantify  area  in  relative  terms  as  a  percentage  of  the  total 
landscape  area.  Relative  area  or  percent  of  the  landscape  can  be  determined  with  the 
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equation: 

n 

z°; 

%LAND  =  ^ —  (1 00)  (1) 

A 

Figure  2  depicts  3  sample  landscapes  that  vary  in  the  amount  and  pattern  of 
agricultural  land,  deciduous  forest,  evergreen  forest,  mixed  forest,  and  forested  wetlands. 
Site  A  is  the  most  fragmented  site,  relative  to  the  original  forest,  while  site  C  is  largely 
forested.  Site  B  is  intermediate.  The  dynamics  of  some  ecological  processes  are  likely  to 
be  quite  different  in  the  3  landscapes.  For  example,  populations  of  organisms  associated 
with  forest  land  are  likely  to  be  much  smaller  in  landscape  A  and  perhaps  subject  to  a 
higher  probability  of  local  extinction  than  in  either  landscape  B  or  C. 

Forest  area  is  particularly  important  to  vertebrate  abundance.  Robbins  et  al. 

(1989)  found,  from  15  variables,  percent  forest  most  often  related  to  bird  relative 
abundance.  Robinson  et  al.  (1995)  found  strong  correlations  between  percent  forest 
cover  and  nesting  success  for  most  species  studied,  more  so  than  the  correlation  with 
percent  forest  interior,  and  mean  forest  patch  size. 

2)  Core  Area 

Core  area  is  the  area  within  a  patch  beyond  a  specified  edge  or  buffer  distance. 

The  core  area  indices  integrate  into  a  single  measure  the  affects  of  patch  area,  patch 
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shape,  and  edge  effects.  Thus,  even  though  they  quantify  landscape  composition  they 
are  also  affected  by  landscape  configuration.  The  primary  significance  of  core  area  of 
patches  in  a  landscape  appears  to  be  related  to  the  edge  effect.  Some  birds,  for  example, 
are  adversely  affected  by  predation,  competition,  and  brood  parasitism  along  forest  edges. 
Core  area  has  been  found  to  be  a  better  predictor  of  habitat  quality  than  patch  area  for 
these  forest  interior  specialists  (Temple  1986).  Core  area  is  affected  by  patch  shape, 
while  patch  area  is  not  affected  by  shape.  Therefore,  even  if  a  patch  is  large  enough  to 
support  a  specific  species,  it  may  not  contain  enough  suitable  core  area  to  support  the 
species. 

Core  area  metrics  may  be  useful  in  the  study  of  habitat  fragmentation  since 
fragmentation  affects  both  habitat  area  and  configuration  (McGarigal  and  Marks  1994). 
At  the  same  time,  these  indices  confound  the  effects  of  habitat  area  and  configuration. 

As  an  example,  if  the  percent  of  a  landscape  that  is  comprised  of  core  area  is  small,  it 
indicates  that  there  is  little  core  area  available,  but  it  does  not  discriminate  between  a 
small  amount  of  the  patch  type  and  a  large  amount  of  the  patch  type  in  a  highly 
fragmented  configuration.  Thus,  core  area  indices  are  usually  best  interpreted  in 
conjunction  with  other  area  indices  that  provide  a  better  overall  picture  of  the  landscape 
structure  (McGarigal  and  Marks  1994). 

Core  area  indices  are  affected  by  the  specified  edge  width.  Since  the  edge  width 
dictates  the  core  area,  the  index  is  meaningful  only  if  the  specified  edge  width  is  relevant 
to  the  research  study.  Also,  the  utility  of  the  core  area  metrics  as  compared  with  the  area 
metrics  is  dependent  on  the  minimum  patch  size  and  edge  width  used.  If  the  minimum 
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patch  size  is  large  but  a  very  small  edge  width  is  used,  the  core  area  and  area  metrics 
will  be  nearly  identical.  Thus  the  core  area  indices  will  be  relatively  insensitive  to 
differences  in  patch  size  and  shape. 

The  total  core  area  index  (TCAI)  quantifies  core  area  for  forest  land  as  a 
percentage  of  total  forest  land.  The  equation  to  calculate  TCAI  is  listed  below. 

TCAI  =  — — (100)  (2) 

i>, 

y=i 

The  TCAI  represents  the  landscape  along  a  continuum  from  most  to  least 
fragmented  (McGarigal  and  Marks  1994).  Figure  3  depicts  3  sample  landscapes  that  vary 
in  the  amount  and  pattern  of  core  forest  land  based  on  a  250  meter  edge  width  from  the 
forest  edge.  According  to  this  index,  only  about  14%  of  the  forest  land  in  landscape  A  is 
considered  interior  habitat;  thus,  the  remaining  86%  is  edge  habitat.  Forest  land  in 
landscape  B  is  50%  interior  habitat  and  50%  edge,  while  72%  of  the  forest  land  in 
landscape  C  is  interior  and  the  remaining  28%  is  edge.  Without  any  other  information,  it 
could  be  deduced  that  the  forest  land  in  landscape  A  is  highly  fragmented.  It  is  useful  to 
know  the  percent  of  the  landscape  comprised  of  forest  land  to  know  if  the  TCAI  for 
landscapes  B  and  C  comprise  a  large  or  small  part  of  the  landscape.  For  example, 
knowing  that  landscapes  B  and  C  are  composed  of  approximately  61%  and  77%  forest 
respectively,  it  is  apparent  that  the  50%  and  72%  of  these  areas  that  are  considered  core 
area  represent  substantial  land  areas. 
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Core  area  is  important  to  vertebrate  abundance  or  occurrence.  Temple  (1986) 
found  total  core  area  to  be  highly  correlated  to  bird  species  abundance,  more  so  than  the 
correlation  with  the  total  area  of  forest. 

3)  Patch  Density 

Metrics  representing  the  number  or  density  of  patches  and  the  average  size  of 
patches  are  not  spatially  explicit  measures  but  are  usually  best  considered  as  representing 
landscape  configuration.  Patch  density  (PD)  has  the  same  utility  as  the  number  of 
patches  index,  except  that  patch  density  presents  the  number  of  patches  on  a  per  unit  area 
basis  allowing  comparisons  among  landscapes  of  varying  size.  The  number  or  density  of 
patches  of  a  particular  habitat  type  may  affect  a  variety  of  ecological  processes.  For 
example,  in  species  that  are  exclusively  associated  with  a  single  habitat  type,  the  patch 
density  may  influence  the  number  of  subpopulations  in  a  spatially  dispersed  population. 
The  number  of  subpopulations  could  influence  the  dynamics  and  persistence  of  the 
metapopulation  (Gilpin  and  Hanski  1991).  The  number  of  patches  or  patch  density  can 
affect  the  stability  of  species  interactions  and  opportunities  for  coexistence  in  predator- 
prey  and  competitive  systems  (Kareiva  1987).  Also,  the  number  of  patches  or  patch 
density  illustrates  habitat  subdivision  which  may  affect  the  spread  of  disturbances  across 
a  landscape  (Franklin  and  Forman  1987).  For  example,  patch  types  that  are  more 
subdivided  might  be  more  resistant  to  disturbances  such  as  disease  and  fire  than 
contiguous  patch  types;  however,  they  may  have  higher  rates  of  disturbance  for 
disturbances  such  as  windthrow  (McGarigal  and  Marks  1994). 


The  equation  used  to  calculate  PD  is  listed  below. 
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PD  =  —  (10,000)(100)  (3) 

A 

Figure  4  depicts  3  sample  landscapes  that  vary  in  the  amount  and  pattern  of  forest 
land.  PD  indicates  that  the  forest  land  is  more  subdivided  in  landscape  A  than  landscape 
C. 

Patch  density  is  important  to  bird  abundance.  McGarigal  and  McComb  (1995) 
found  a  component  defined  by  patch  density  to  be  significant  to  certain  species 
abundance,  while  components  depicting  patch  size  and  patch  shape  did  not  show  a 
significant  association  to  the  bird  species  abundance. 

4)  Patch  Size 

As  discussed  in  the  area  metrics,  the  area  of  the  patches  comprising  a  landscape 
is  one  of  the  most  important  pieces  of  information  contained  in  the  landscape  (McGarigal 
and  Marks  1994).  The  average  size  comprised  by  each  patch  type  is  also  important. 
Progressive  reduction  in  the  size  of  habitat  fragments  is  a  key  component  of  habitat 
fragmentation.  Therefore,  a  landscape  that  has  a  smaller  average  patch  size  than  another 
landscape  is  considered  more  fragmented. 

Mean  patch  size  represents  the  average  size  of  all  patches  of  a  particular  type 
within  the  landscape.  It  is  derived  from  the  number  of  patches  but  does  not  convey  any 
information  about  how  many  patches  are  present  in  the  landscape.  Since  a  mean  patch 
size  of  5  ha  could  represent  1  or  50  patches,  which  could  have  important  ecological 
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implications,  it  is  probably  best  interpreted  along  with  patch  density  (McGarigal  and 
Marks  1994). 

The  minimum  patch  size  and  the  extent  of  the  image  influence  the  mean  patch 
size.  Patches  are  often  subdivided  by  the  extent  boundaries  (McGarigal  and  Marks 
1994). 

The  equation  used  to  calculate  MPS  is  listed  below. 


n 


(4) 


Figure  4  depicts  3  sample  landscapes  that  vary  in  the  amount  and  pattern  of  forest 
land.  MPS  attempts  to  rank  the  3  landscapes  with  respect  to  forest  land  fragmentation, 
with  A  being  most  fragmented  and  C  being  the  least  fragmented.  The  examples  used  in 
Figure  4  have  dramatic  differences  in  MPS.  By  interpreting  the  MPS  with  patch  density, 
it  is  evident  that  landscape  A  was  very  fragmented,  with  small  numerous  patches  while 
landscape  C  was  not  very  fragmented,  having  a  smaller  number  of  large  patches. 

The  average  patch  size  of  forest  patches  is  important,  particularly  to  wildlife 
management.  Howell  et  al.  (2000)  found  a  high  correlation  between  forest  interior  bird 
species  abundance  and  mean  size  of  forest  patches,  more  so  than  percent  forest  and  edge 
density. 

5)  Edge  Density 

Edge  metrics  are  considered  to  best  represent  landscape  configuration,  even  though  they 
are  not  spatially  explicit.  The  amount  of  edge  in  the  landscape  is  pertinent  to  various 
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ecological  criteria,  such  as  that  of  wildlife  abundance.  The  forest  edge  effect  is  due  to 
the  landscape  having  differences  in  wind  and  light  intensity  that  reaches  the  forest  patch 
which,  in  turn  alters  the  microclimate  and  disturbance  rates  (e.g.,  Ranney  et  al.  1981). 

The  patch  shape  and  orientation  and  the  adjacent  land  cover  affect  the  proportion  of  the 
patch  that  is  affected  by  edge.  A  large  round  patch  has  a  minimal  amount  of  edge  habitat 
while  a  large  convoluted  patch  may  be  primarily  edge  habitat.  It  is  generally  accepted 
that  edge  effects  are  viewed  from  an  organism  perspective,  since  species  usually  have 
either  an  affinity,  adversity,  or  are  unaffected  by  edge  (McGarigal  and  Marks  1994). 

Studies  have  suggested  that  changes  in  vegetation,  predation,  brood  parasitism, 
and  competition  along  forest  edges  has  resulted  in  the  reduction  of  populations  of  various 
vertebrate  species  dependent  upon  forest  interiors  (e.g.  Robbins  et  al.  1989,  Wilcove 
1985). 

Edge  density  (ED)  standardizes  edge  to  a  per  unit  area  basis  so  that  comparisons 
among  landscapes  of  varying  size  can  be  made.  The  index  is  affected  by  the  resolution  of 
the  image.  The  greater  the  detail  with  which  the  edges  are  delineated,  the  greater  the 
apparent  edge  length.  For  example,  in  images  with  a  coarse  resolution,  edges  may  mostly 
appear  as  straight  lines,  while  with  finer  resolutions,  the  edges  may  appear  as  more 
convoluted  lines.  Edge  metrics  developed  from  images  with  different  resolutions  should 


not  be  compared. 


The  equation  used  to  calculate  ED  is  listed  below. 
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IX 

ED  =  -^L_  (10,000)  (5) 

Figure  4  shows  3  sample  landscapes  that  vary  in  the  amount  and  pattern  of  forest 
land.  ED  is  highest  in  landscape  A.  For  a  species  that  requires  forest  land  edge  habitat, 
ED  might  be  used  to  model  habitat  suitability.  In  this  case,  landscape  A  would  be  most 
suitable  while  landscape  C  would  be  least  suitable. 

Edge  density  is  important  in  determining  response  by  vertebrates  to  forest 
fragmentation  or  the  edge  and  interior  of  forests.  Rosenberg  et  al.  (1999)  found  a  strong 
relationship  between  the  chance  of  detecting  Brown-headed  Cowbirds  and  potential  nest 
predators  and  highly  fragmented  sites  that  had  high  edge  density  or  large  amounts  of 
edge.  This  response  of  Brown-headed  Cowbirds  and  potential  nest  predators  was 
opposite  that  of  Tanagers. 

6)  Mean  Shape  Index 

Shape  metrics  quantify  landscape  configuration  in  terms  of  the  complexity  of  the 
patch  shape.  Patch  shape  is  possibly  as  important  as  patch  size;  however,  there  is 
relatively  little  information  about  the  effects  of  shape  on  the  ecosystem  (Forman  and 
Godron  1986).  Patch  shape  is  important  in  the  dispersal  and  foraging  of  organisms.  For 
example,  birds  flying  over  woods,  are  more  apt  to  find  a  long  narrow  clearing  that  is 
oriented  perpendicular  to  their  direction  of  movement,  while  they  may  miss  a  round 
clearing  (Forman  and  Godron  1986).  The  long  narrow  clearing  parallel  to  their 
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movement  may  also  be  missed.  Patch  shape  has  been  shown  to  influence  inter-patch 
processes  such  as  small  mammal  migration  (Buechner  1989).  However,  the  primary 
significance  of  the  shape  of  patches  in  a  landscape  appears  to  be  related  to  the  edge 
effect. 

One  method  of  quantifying  patch  shape  is  by  assessing  the  complexity  of  the 
patch  shape  compared  to  a  standard  circular  shape.  The  shape  index  is  minimum  if  the 
patch  is  circular  and  increases  as  the  patches  become  noncircular.  This  shape  index  is 
widely  applicable  and  used  in  landscape  ecological  research  (Forman  and  Godron  1986). 
The  mean  shape  index  (MSI)  measures  the  average  perimeter-to-area  ratio  for  a  specific 
patch  type.  This  shape  index  is  limited  in  the  same  way  as  the  edge  indices  in  reference 
to  the  differences  between  how  lines  are  portrayed  in  vector  and  raster  images.  The 
perimeter-to-area  ratio  method  in  assessing  shape,  is  insensitive  to  differences  in  form 
and  structure  or  patch  morphology  (McGarigal  and  Marks  1994).  Patches  may  have 
different  shapes,  but  may  have  identical  areas  and  perimeters  and  shape  indices.  This 
index  is  not  useful  as  a  measure  of  patch  morphology,  but  is  best  considered  as  a  measure 
of  overall  shape  complexity.  This  index  represents  the  average  patch  shape  for  a  specific 
habitat,  and  thus  does  not  necessarily  fully  describe  the  shapes  of  patches  in  a  landscape 
if  the  distribution  of  patch  shapes  is  complex. 

The  equation  used  to  calculate  MSI  is  listed  below. 
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Figure  4  shows  3  landscapes  that  vary  in  the  amount  and  pattern  of  forest  land. 
The  MSI  values  for  all  3  landscapes  are  greater  than  1,  indicating  the  average  patch  shape 
in  all  3  landscapes  is  noncircular.  The  forest  land  patches  in  landscapes  B  and  C  (least 
fragmented)  are  most  irregular,  while  the  patches  in  landscape  A  (most  fragmented)  are 
least  irregular  in  shape.  These  results  indicate  that  human-induced  fragmentation  in 
landscape  A  caused  a  simplification  in  patch  shapes  compared  to  the  geometrically 
complex  patch  shapes  found  in  the  more  natural,  unaltered  landscapes  B  and  C. 

The  shape  of  forest  stands  and  openings  is  particularly  important  for  wildlife, 
because  habitat  conditions  and  shape  are  often  closely  interrelated  (Marcot  and  Meretsky 
1983).  Palmer  et  al.  (2000)  even  used  MSI  in  a  highly  dynamic  streambed  landscape  and 
found  MSI  to  be  linked  to  faunal  abundance.  Garrabou  et  al.  (1998)  found  MSI  to  be  one 
of  the  most  suitable  indices  for  describing  its  spatial  pattern  in  rocky  benthic 
communities. 

7)  Diversity 

The  diversity  metrics  quantify  landscape  composition.  Diversity  measures  have 
been  used  in  a  variety  of  ecological  applications,  such  as  measures  of  plant  and  animal 
species  diversity.  They  are  commonly  used  in  community  ecology.  Diversity  measures 
are  influenced  by  richness  and  evenness,  with  richness  pertaining  to  the  number  of  patch 
types  present  and  evenness  pertaining  to  the  distribution  of  area  among  different  types. 
Some  indices  are  more  sensitive  to  richness  while  others  are  more  sensitive  to  evenness. 
The  actual  species  composition  of  a  community  is  not  provided  in  the  diversity  indices. 
For  example,  a  community  could  have  high  species  diversity  and  be  comprised  primarily 


24 


of  common  or  undesirable  species.  Or,  a  community  could  have  low  species  diversity 
and  be  comprised  of  rare  species. 

Two  common  diversity  indices  are  the  Simpson's  diversity  index  and  Shannon's 
diversity  index.  Simpson's  diversity  index  (SIDI)  is  less  sensitive  to  richness  and 
therefore  places  more  weight  on  the  common  species.  Shannon's  diversity  index  is  more 
sensitive  to  richness  and  thus  places  more  weight  on  rare  species.  A  wide  variety  of 
diversity  indices  have  been  used  to  measure  landscape  composition  (O'Neill  et  al.  1988, 
Turner  1990).  SIDI  was  used  in  this  study. 

The  value  of  SIDI  represents  the  probability  that  any  patch  type  selected  at 
random  would  be  different  types.  The  higher  the  value,  the  greater  the  likelihood  that 
any  two  randomly  drawn  patches  would  be  different  patch  types  or  greater  diversity. 

The  equation  used  to  calculate  SIDI  is  listed  below. 

rn 
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Figure  5  shows  3  landscapes  that  vary  in  composition  and  pattern.  The  SIDI  represents 
the  landscapes  along  a  continuum  from  most  to  least  diverse.  The  diversity  of  landscapes 
A  and  B  is  similar,  while  landscape  C  is  more  diverse.  In  landscape  A,  SIDI  indicates 
that  there  is  a  54%  probability  that  2  randomly  chosen  patches  would  represent  different 
patch  types. 

Species  diversity  of  a  landscape  patch  appears  to  be  mainly  determined  by  habitat 
diversity  and  the  disturbance  regime  (Forman  and  Godron  1986).  In  using  SIDI  to 
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determine  effects  of  landscape  fragmentation  on  bird  communities,  it  was  found  that 
continuous  forest  landscapes  supported  more  species  than  did  fragmented  landscapes 
McIntyre  (1995). 

8)  Interspersion  and  juxtaposition 

Interspersion  and  juxtaposition  metrics  quantify  landscape  configuration.  Each 
patch  is  evaluated  for  adjacency  with  all  other  patch  types.  The  interspersion  index 
measures  the  extent  to  which  patch  types  are  interspersed,  not  necessarily  dispersed. 
Higher  values  result  from  landscapes  in  which  the  patch  types  are  well  interspersed  or 
equally  adjacent  to  each  other.  Lower  values  characterize  landscapes  in  which  the  patch 
types  are  poorly  interspersed.  The  index  is  not  directly  affected  by  the  number,  size, 
contiguity,  or  dispersion  of  patches.  The  interspersion  index  is  a  relative  index  that 
represents  the  observed  level  of  interspersion  as  a  percentage  of  the  maximum  possible 
given  the  total  number  of  patch  types  (McGarigal  and  Marks  1994). 

The  equation  used  to  calculate  UI  is  listed  below. 
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Figure  5  shows  3  landscapes  that  vary  in  composition  and  pattern.  The  UI  values  indicate 
that  the  interspersion  of  available  patch  types  is  greatest  in  landscape  C  and  least  in 
landscape  B.  This  probably  occurs  because  landscape  B  has  patch  types  that  are  very 
small,  thus  the  distribution  of  edge  lengths  among  the  unique  patch  types  is  more  uneven. 
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UI  is  used  in  a  variety  of  purposes.  Hunziker  (1999)  found  the  UI  measures 
correlate  significantly  with  landscape  preference  values  in  terms  of  illustrating  landscape 
fragmentation.  Kean  et  al.  (1999)  used  UI  to  measure  landscape  pattern  in  areas  that 
have  burned  and  found  fire  influences  landscape  pattern  by  creating  more  fragmented, 
and  disconnected  landscapes. 

Horizontal  Spatial  Pattern  Related  to  Bird  Abundance 

Songbirds  are  an  almost  ideal  subject  for  investigating  the  potential  use  of  spatial 
metrics  in  the  evaluation  of  habitat  suitability.  Songbirds  have  the  advantage  of  being 
widely  dispersed  and  are  known  to  be  sensitive  to  the  same  landscape  changes  that  can  be 
remotely  sensed,  and  for  which  data  are  widely  available.  In  the  relatively  sparse 
literature  existing  on  this  subject  there  is  general  agreement  that  area  metrics  (e.g., 
percent  forest  cover)  provide  a  simple  first  approximation  for  exploring  habitat 
suitability.  Several  studies  report  good  correlations  with  the  relative  amount  of  suitable 
habitat  area  and  bird  success  (Robinson  et  al.  1995,  Pearson  1993,  Robbins  et  al.  1989, 
Roberts  and  Norment  1999,  Knick  and  Rotenberry  1995).  However,  total  area  may  be 
less  important  than  its  distribution  within  the  landscape.  The  Robinson  et  al.  (1995) 
study  notes  that  nesting  success  of  the  birds  studied  was  positively  correlated  with  forest 
area  but  that  forest  area  was  autocorrelated  with  percent  forest  interior.  A  high  percent 
forest  interior  discourages  predatory  Brown-headed  Cowbirds  and  offers  a  better 
biological  explanation  of  the  bird  success  than  forest  area  alone.  Thus  core  forest  area 
has  become  a  popular  indicator  for  interior  forest  species  (Temple  1986,  Donovan  et  al. 
1997,  McGarigal  and  McComb  1995).  Some  Neotropical  migrant  forest-interior  birds 
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have  displayed  an  adverse  sensitivity  to  the  creation  of  edges  within  large  forest 
interiors,  indicating  that  they  are  limited  by  edge  zones  rather  than  forest  fragment  size 
(Germaine  et  al.  1999).  Proximity  to  a  habitat  edge  may  greatly  reduce  the  desirability  of 
a  portion  of  a  forest,  especially  for  the  single-brooded  low-nesting  neotropical  migrants 
(Whitcomb  et  al.  1981).  Donovan  et  al.  (1997)  found  that  the  extent  of  edge  effects  in 
nest  predation  in  the  central  United  States  was  related  to  the  degree  of  habitat 
fragmentation.  Rates  of  predation  on  artificial  nests  were  high  in  both  edge  and  core 
forest  habitat  in  highly  fragmented  landscapes,  and  low  in  both  habitats  in  unfragmented 
landscapes.  Significant  edge  effects  were  found  in  moderately  fragmented  landscapes, 
where  rates  of  predation  were  high  along  edges  and  low  in  the  core  habitat.  Heske  et  al. 
(1999)  also  found  that  medium-sized,  generalist  mammalian  predators  on  songbird  nests 
reach  their  highest  population  densities  in  fragmented  landscapes  with  abundant  edge 
habitat,  particularly  agricultural  edges.  The  increased  nest  predation  and  parasitism  in 
combination  with  isolation  from  other  forests  may  reduce  both  the  rates  of  return  by  adult 
birds  and  colonization  by  first-time  breeders  (Robbins  et  al.  1989).  These  points  argue 
for  the  consideration  of  patch  shape  and  interspersion,  as  well  as  patch  size  and  area.  It 
has  been  demonstrated  that  the  area  of  forest  interior  increases  at  a  steeper  rate  than  edge, 
for  nearly  square  or  circular  patches  (Levenson  1981).  Contrary  to  forest  interior  species, 
the  distribution  of  some  bird  species  may  be  insensitive  to  area  metrics  due  to  a 
dependence  on  edge  (Bolger  et  al.  1997,  Rosenberg  et  al.  1999). 

Since  many  landscape  metrics  are  highly  correlated  with  each  other,  principal 
components  analysis  has  been  used  to  deal  with  the  problems  associated  with 
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multicollinearity.  Various  studies  have  utilized  this  technique  in  addressing  the  issue 
of  autocorrelation  and  have  been  successful  in  developing  meaningful  factors  that  can  be 
used  to  describe  the  components  of  the  landscape  (McGarigal  and  McComb  1995, 
Osborne  1984,  Rosenberg  et  al.  1999). 

From  reviewing  these  analyses  it  is  clear  that  the  appropriate  approach  for 
quantifying  habitat  suitability  will  vary  from  species  to  species  depending  upon  their 
particular  environmental  adaptations.  Also,  there  may  be  cases  where  species  are 
insensitive  to  variation  in  the  landscape.  For  example,  generalist  species  utilize  a  variety 
of  habitat  types,  so  the  landscape  effect  is  weaker  than  for  more  specialized  species 
(Pearson  1993).  A  specific  example  is  the  Northern  Parulas  ( Parula  americana)  which 
are  relatively  insensitive  to  variation  in  landscape  structure  but  are  highly  sensitive  to 
arboreal  poison  ivy  (Robins  et  al.  1989).  While  landscape  structure  is  just  one  of  many 
possible  limitations  to  bird  success,  it  is  certainly  an  important  factor  and  is  perhaps  the 


one  most  under  human  control. 


Figure  2.  Landscape  composition  metrics  for  three  sample  landscapes  that  vary  in  the  amount  and  pattern  of  agriculture 
and  forest  land. 
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Figure  3.  Total  core  area  index  for  three  sample  landscapes  that  vary  in  the  amount  and  pattern  of  forest  land. 


Figure  4.  Spatial  pattern  metrics  for  three  sample  landscapes  that  vary  in  the  amount  and  pattern  of  forest  land. 


32 


Figure  5.  Interspersion  and  juxtaposition  index  (UI)  and  Simposon’s  diversity  index  (SIDI)  for  three 
sample  landscapes. 


CHAPTER  III 


METHODS 

Study  Design  Overview 

The  objective  of  this  research  was  to  use  the  BBS  and  USGS  LULC  data  to 
determine  if  kilometer-resolution  horizontal  spatial  pattern  metrics  are  suitable  indicators 
of  habitat  suitability  for  Neotropical  and  resident  songbirds.  The  mean  count  of  each  bird 
/route  was  calculated  using  BBS  data.  This  was  done  for  15  bird  species  and  3  bird 
groups  for  the  53  BBS  routes  shown  in  Figure  1.  The  landscape  structure  was  quantified 
using  a  suite  of  12  spatial  metrics  calculated  from  USGS  LULC  data  for  each  of  the  53 
BBS  routes.  The  metrics  for  all  routes  were  pooled  and  summarized  into  broad  landscape 
structure  components  using  principal  components  analysis  (PC A)  techniques.  Multiple 
regression  techniques  were  used  to  determine  if  kilometer-resolution  horizontal  spatial 
pattern  metrics  are  suitable  indicators  of  habitat  suitability. 

The  15  species  selected  for  inclusion  in  this  research  (Table  1)  were  identified  as 
conservation  priority  species  in  prior  research,  Partners  in  Flight  (Flunter  1998).  The 
designation  of  a  species  as  a  conservation  priority  species  indicates  that  its  situation  can 
be  considered  as  particularly  precarious  and  that  it  should  be  the  focus  of  conservation 
efforts.  According  to  the  BBS  data  the  incidence  of  most  of  these  species  has  been 
declining,  either  on  an  overall  basis  or  in  localized  areas  (Appendix  A).  Habitat  loss  or 
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the  fragmentation  of  existing  habitat  is  often  cited  as  a  primary  reason  for  this 
population  decline.  However,  this  belief  is  usually  based  on  intuitive  opinion  or  on  local 
level  studies.  While  habitat  loss  and  fragmentation  may  indeed  be  the  root  cause  of  the 
decline,  there  is  not  a  sufficient  empirical  basis,  at  the  landscape  level,  to  provide  a  useful 
understanding  of  this  relationship.  Presently  there  is  no  reliable  procedure  for  assessing 
or  quantifying  the  impact  of  landscape  changes  on  habitat  suitability,  even  for  these 
conservation  priority  species. 

This  research  effort  is  based  on  the  premise  that  conservationists  need  a  method 
of  determining  landscape  level  indicators  of  habitat  suitability  for  these  conservation 
priority  species  and  that  such  a  method  can  be  developed  using  the  existing  the  BBS  and 
USGS  LULC  data.  Specifically  conservationists  need  horizontal  spatial  pattern  metrics 
that  provide  indicators  of  habitat  suitability  for  these  bird  species. 

Out  of  the  over  30  breeding  or  resident  landbirds  listed  by  Partners  in  Flight  as 
conservation  priority  for  this  study  area,  14  were  selected.  The  Brown-headed  Cowbird, 
which  is  not  a  conservation  priority  specie,  was  also  included  since  it  has  been  studied 
and  found  to  be  associated  with  nesting  success.  All  the  conservation  priority  species 
were  assessed,  but  only  those  that  were  well  represented  on  the  BBS  routes  in  the  study 
area  were  included  in  this  study.  The  criterion  for  selection  was  the  birds  presence  on 
approximately  two  thirds  of  the  BBS  routes.  This  eliminated  species  that  are  typically 
poorly  sampled  by  the  BBS,  such  as  nocturnal  birds,  seabirds,  shorebirds,  and  raptors 
(Hepinstall  and  Sader  1997).  After  eliminating  the  species  that  were  not  well  represented 
or  typically  poorly  sampled,  14  bird  species  remained.  These  14  species  represent  habitat 
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preferences  that  range  from  forested  areas  to  successional  scrub.  The  majority  of  the 
species  are  Neotropical  migrants,  with  a  few  resident  bird  species. 

The  development  of  quantitative  data  describing  landscape  structure  required  the 
selection  of  landscape  spatial  pattern  metrics  and  the  determination  of  the  physical  extent 
on  which  these  metrics  were  to  be  calculated.  The  landscape  spatial  pattern  metrics 
selected  for  this  study  are  listed  in  Table  4.  They  are  a  representative  sample  of  metrics 
from  the  broad  categories  of  area,  core  area,  patch  density  and  size,  edge,  shape,  diversity 
and  interspersion  metrics.  These  landscape  spatial  pattern  metrics  were  calculated  for 
four  physical  extents  around  each  BBS  route.  Figures  6  and  7  show  a  typical  BBS  route 
and  the  four  extents  around  a  BBS  route,  respectively.  The  four  extents  were  defined  by 
a  narrow  (0.4  km)  and  wide  (10  km)  radius  buffer  around  each  full  BBS  route  as  well  as  a 
narrow  and  wide  buffer  around  each  BBS  internal  5  mile  sub-segment.  These  sub- 
segments  are  referred  to  as  partial  routes  (Figure  6). 

A  detailed  description  of  the  methods  used  in  this  research  is  presented  in  the 
remaining  sections  of  this  chapter.  Figure  8  shows  a  schematic  outline  of  the  sequence 
of  procedures  used  to  develop  the  quantitative  data  and  then  to  determine  the  efficacy  of 
kilometer-resolution  horizontal  spatial  pattern  metrics  as  indicators  of  habitat  suitability. 
The  procedures  in  the  schematic  outline  are  abbreviated  (i.e.  (PI)  indicates  the  first 
procedure).  These  abbreviations  are  used  in  the  text  so  that  one  can  easily  locate  the 
procedure  in  the  overall  project  design.  The  SPSS  software  was  used  for  all  statistical 
analyses  (SPSS  1998)  and  the  Environmental  Systems  Research  Institute  (ESRI  1998) 
ARC/INFO  GIS  software  was  used  for  all  GIS  database  development  and  analysis 
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procedures.  The  FRAGSTATS  Spatial  Pattern  Analysis  Program  for  Quantifying 
Landscape  Structure  was  used,  in  conjunction  with  ARC/INFO,  to  calculate  the  spatial 
pattern  metrics. 


Breeding  Bird  Survey  Data 

The  BBS  is  conducted  by  volunteers  along  secondary  roads,  randomly  selected 
within  degree  blocks  of  latitude  and  longitude.  Each  BBS  route  is  approximately  40-km 
long  and  consists  of  50  stops  spaced  0.8  km  apart  (Sauer  et  al.  1996).  All  birds  seen  or 
heard  within  a  0.4  km  radius  during  three  minutes  are  recorded.  Each  survey  route  is  run 
during  the  peak  of  the  breeding  season,  with  certain  guidelines  for  time  of  day  and 
weather  conditions  intended  to  reduce  biases  in  the  data  (Robbins  et  al.  1986,  Peteijohn 
and  Sauer  1993). 

The  records  for  each  BBS  route  included  each  species  American  Ornithologist 
Union  (AOU)  number  and  the  number  of  individuals  observed  for  the  first  ten  stops, 
second  ten  stops,  third  ten  stops,  fourth  ten  stops,  and  fifth  ten  stops,  and  for  the  full  BBS 
route  (Robbins  et  al.  1986).  It  is  these  5,  10-stop  summaries  or  5  mile  subsegments  that 
form  the  basis  of  the  5  partial  BBS  routes,  and  together,  the  full  BBS  route.  A  diagram  of 
a  typical  BBS  route  is  shown  in  Figure  6.  There  were  53  full  BBS  routes  used  in  this 
study:  33  in  Alabama,  9  in  Florida,  and  1  lin  Georgia.  And  with  each  of  the  53  full  BBS 
routes  broken  into  5  partial  routes,  this  resulted  in  265  partial  BBS  routes  in  the  study. 

Using  a  minimum  of  3  years  of  BBS  data,  the  mean  yearly  abundance  for  each 
species  was  developed  for  each  partial  and  full  BBS  route  (Figure  8,  Step  PI).  The  mean 
yearly  abundance  was  calculated  by  adding  the  abundance  of  the  species  for  each  year 
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that  data  were  available  (3-5  years)  for  each  route  or  partial  route  and  dividing  by  the 
number  of  years  of  data. 

The  5  years  of  data  included  in  this  study  correspond  to  the  year  of  the 
photography  from  which  the  USGS  LULC  were  interpreted,  as  well  as  the  two  years 
prior  to  and  two  years  after  the  date  of  the  photography.  The  additional  non-photograph 
years  of  data  were  included  to  reduce  undue  variation  from  extraneous  factors  such  as 
severe  weather.  Data  from  only  3  or  4  years  were  used  for  some  BBS  routes  as  not  all 
BBS  routes  had  data  for  the  appropriate  5  years.  Table  5  shows  the  years  (from  1970- 
1976)  of  BBS  data  used  to  develop  the  mean  count  of  birds/route  for  each  species. 

The  mean  count  data  were  used  to  place  each  bird  into  a  high,  medium,  and  low 
abundance  category  for  each  route.  The  birds  were  categorized  using  the  same  criterion 
as  was  used  in  the  development  of  range  maps  for  birds  across  the  South  (Hamel  1992). 
In  the  development  of  range  maps  for  birds  of  the  South,  Hamel  utilized  BBS  data  from 
1966-1985.  Hamel  determined  the  highest  mean  count  of  birds/route  across  the  South 
(Table  6).  Areas  with  at  least  30  %  of  the  peak  value  are  considered  High,  those  with  at 
least  10%  are  considered  Medium,  and  those  with  at  least  5%  of  peak  value  are 
considered  Low  Abundance.  This  same  methodology  was  applied  to  mean  count  of 
birds/route  developed  for  this  study.  Categorizing  bird  abundance  in  this  manner  allows 
for  the  differentiation  between  the  relative  abundance  of  a  given  species  in  different 
areas.  This  in  turn  provides  the  means  of  differentiating  horizontal  spatial  pattern  metric 
values  that  should  be  good  indicators  of  habitat  suitability  for  Neotropical  and  resident 
songbirds  based  on  sites  with  high  abundance. 
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Bird  species  were  also  grouped  into  three  categories  in  order  to  study  the 
horizontal  spatial  pattern  metric  values  that  are  suitable  indicators  of  habitat  utilization. 
The  groups  included:  edge  and  scrub  plus  field-edge  bird  species,  forest  interior/edge 
bird  species,  and  forest  interior  bird  species.  These  groups  were  based  on  habitat 

utilization  bird  groups  developed  by  Whitcomb  et  al.  (1981).  The  mean  count  of 
birds/route  for  these  groups  were  developed  by  adding  the  following  bird  species: 

Forest  Interior  Bird  Species 
Hooded  Warbler 
Kentucky  Warbler 
Forest  Interior/Edge  Bird  Species 
Wood  Thrush 
Prothonotary  Warbler 
Eastern  Wood  Pewee 
Yellow-throated  Vireo 
Yellow-billed  Cuckoo 
White-eyed  Vireo 
Brown-headed  Nuthatch 
Edge  and  Scrub  plus  Field-Edge  Bird  Species 
Prairie  Warbler 
Loggerhead  Shrike 
Field  Sparrow 
Eastern  Kingbird 
Northern  Bob  white 
Brown-headed  Cowbird 
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Physical  Extent  of  Spatial  Pattern  Analysis 

A  GIS  coverage  of  each  BBS  route  was  developed  using  USGS  DLG  data  at  a 
scale  of  1:100,000  (U.S.  Geological  Survey,  1989)  (Figure  8,  Step  P2).  Each  full  BBS 
route  was  divided  into  5  equal  parts  that  represented  the  5  partial  BBS  routes.  The  full 
BBS  route  and  the  5  partial  routes  served  as  the  basis  on  which  the  narrow  and  wide 
buffers  were  developed  (Figure  8,  Step  P3)  (Figure  7).  The  narrow  buffer  had  a  0.4  km 
radius,  the  distance  that  a  surveyor  can  hear  or  see  a  bird  at  a  BBS  stop.  The  wide  buffer 
had  a  10  km  radius,  a  width  used  in  prior  research  (Robinson  et  al.  1995).  As  used  in  this 
report,  the  term  “Extent  1”  is  a  0.4  km  radius  buffer  of  a  full  BBS  route;  “Extent  2”  is  a 
10  km  radius  buffer  of  a  full  BBS  route;  “Extent  3”  is  a  0.4  km  radius  buffer  of  each 
partial  BBS  route;  and  “Extent  4”  is  a  10  km  radius  buffer  of  each  partial  BBS  route.  The 
physical  area  of  each  extent  for  each  BBS  route  varied  somewhat,  since  the  BBS  routes 
had  different  shapes.  Typically,  the  area  for  Extent  1  was  3100  ha  and  was  comprised  of 
an  average  of  50  patches;  Extent  2  was  100,000  ha  and  was  comprised  of  an  average  of 
370  patches;  Extent  3  was  675  ha  and  was  comprised  of  an  average  of  13  patches;  and 
Extent  4  was  45,000  ha  and  was  comprised  of  an  average  of  174  patches.  Figure  7 
illustrates  each  extent.  Note  that  there  is  a  small  overlap  in  the  buffer  areas  for  Extent  3 
and  a  large  overlap  in  the  buffer  areas  for  Extent  4. 


US  Geological  Survey  Land  Use  Land  Cover  Data 

The  physical  extents  defined  above  were  overlain  USGS  LULC  data.  The  USGS 
provides  LULC  data  in  digital  format  in  1:250,000  quadrangles,  each  covering  an  area  of 
1  degree  of  latitude  by  2  degrees  of  longitude.  The  USGS  LULC  quadrangles  in  the 


40 


study  area  were  interpreted  from  NASA  U2/R8-57  high  altitude  aerial  photography 
taken  in  1972,  1973,  and  1974  (US  Geological  Survey  1990).  Figure  1  shows  the  outline 
and  name  of  the  USGS  1 :250, 000-scale  quadrangles  and  the  estimated  date  of  the 
photography  from  which  they  were  interpreted.  Twelve  1:250,000  quadrangles  were 
used  to  incorporate  all  BBS  routes. 

The  USGS  LULC  data  represent  a  national  classification  scheme  that  has 
achieved  widespread  acceptance  and  is  being  used  in  a  number  of  operational  mapping 
programs  (Avery  and  Berlin  1992).  Table  7  lists  the  patch  types  in  Level  I  and  II  of  the 
USGS  LULC  dataset  (Fegeas  et  al.  1983).  Each  polygon  in  the  LULC  dataset  is  referred 
to  in  this  report  as  a  "patch".  The  patch  type  is  a  descriptive  term  attached  to  each  patch, 
or  basic  unit  that  makes  up  the  landscape  (Urban  et  al.  1987).  All  patch  types  have  a 
minimum  polygon  size  of  16  ha,  except  urban,  built-up  land,  water,  confined  feeding 
operations,  other  agricultural  land  and  strip  mines,  quarries,  and  gravel  pits,  which  have  a 
minimum  polygon  size  of  4  ha. 

For  each  of  the  twelve  1 :250,000  quadrangles  of  digital  LULC  data,  two  different 
ESRI  Arc/Info  GIS  coverages  were  developed.  In  the  first  coverage,  the  original  Level  II 
patch  types  were  retained  (Figure  8,  Step  P4).  In  the  second  coverage,  the  deciduous 
forest,  mixed  forest,  evergreen  forest,  and  forested  wetland  patch  types  were  reclassified 
into  a  general  patch  type  called  forest,  with  all  line  segments  depicting  the  original 
boundaries  between  the  different  forest  patch  types  removed  (Figure  8,  Step  P5).  Once 
the  two  GIS  coverages  were  developed  for  each  USGS  LULC  quadrangle,  the  4  physical 
extents  defined  above,  were  extracted  from  the  coverages  (Figure  8,  Step  P6). 
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Spatial  Pattern  Metrics 

FRAGSTATS  Spatial  Pattern  Analysis  Program  for  Quantifying  Landscape 
Structure  (McGarigal  and  Marks  1994)  was  used  to  calculate  the  landscape  spatial  pattern 
metrics  separately  for  each  of  the  4  extents  around  each  BBS  route  (Figure  8,  Step  P7). 
Table  4  lists  pertinent  information  about  each  metric,  including  whether  it  describes 
composition  or  configuration,  its  range  in  value,  a  description  of  how  the  metric  is 
calculated,  and  the  patch  type  on  which  the  calculation  was  based. 

As  shown  in  Table  4,  the  composition  metrics,  percent  agriculture,  percent 
deciduous  forest,  percent  evergreen  forest,  percent  mixed  forest,  and  percent  forested 
wetland,  were  based  on  specific  level  II  patch  types.  The  configuration  metrics,  mean 
patch  size,  total  core  area  index,  patch  density,  edge  density,  and  mean  shape  index,  were 
based  on  the  generalized  forest  patch  type.  The  metrics  developed  for  the  entire 
landscape  mosaic  where  all  patch  types  were  considered  simultaneously  were  Simpson’s 
diversity  index  (SDI)  and  the  interspersion/juxtaposition  index  (UI).  Table  4  shows  that 
these  metrics  were  based  on  all  level  II  patch  types. 

The  configuration  metrics,  MPS,  TCAI,  PD,  ED,  and  MSI,  were  developed  on  the 
generalized  forest  patch  type.  This  was  done  for  several  reasons.  First,  one  of  the  intents 
of  the  study  was  to  determine  if  there  was  an  association  of  bird  species  with  landscape 
forest  fragmentation,  not  fragmentation  of  each  specific  forest  patch  type.  Secondly,  had 
the  specific  forest  patch  types  been  used,  this  would  have  resulted  in  a  total  of  27  metrics. 
Summarizing  this  many  metrics  into  meaningful  components  would  have  proven 
difficult.  Finally,  had  the  specific  forest  categories  been  used,  many  individual  study 
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sites  would  have  been  eliminated  from  subsequent  statistical  analyses.  This  was  due  to 
the  fact  that  "0"  is  not  a  valid  value  for  configuration  metrics  such  as  PD,  MPS,  and  MSI. 
These  metrics  are  always  greater  than  "0"  if  the  patch  type  is  present.  If  the  patch  type  is 
not  present,  a  "no  data"  value  is  produced.  A  "no  data"  value  would  eliminate  the  site 
from  statistical  analyses.  Since  forested  wetlands  and  deciduous  forest  patch  types  were 
absent  from  many  of  the  study  sites  (Table  8),  these  sites  would  have  been  eliminated 
from  subsequent  statistical  analyses. 

For  the  area  metrics,  though,  “0”  is  a  valid  value.  Hence,  the  composition  metrics 
were  based  on  the  more  specific  forest  patch  types  as  well  as  the  cropland  and  pasture 
patch  type.  The  literature  documents  examples  of  using  specific  as  well  as  general  forest 
categories  to  develop  landscape  spatial  pattern  metrics  (e.g.  Robinson  et  al.  1995, 
McGarigal  and  McComb  1995);  however,  the  decision  is  dependent  on  the  research 
objective. 


Database  Statistical  Analysis 

Two  statistical  approaches  were  used  in  this  research  study:  principal  components 
analysis  (PCA)  and  stepwise  multiple  regression.  The  SPSS  was  used  for  all  statistical 
analyses  (SPSS  1998).  Many  of  the  12  landscape  spatial  pattern  metrics  are  known  to  be 
highly  correlated  with  one  another.  A  principal  components  analysis  was  used  to 
simplify  the  structure  of  the  spatial  metric  data  sets  by  reducing  them  to  a  smaller  set  of 
uncorrelated  variables  that  accounted  for  a  large  part  of  the  variation  in  the  original  data 
set.  A  PCA  is  often  useful  in  reducing  the  amount  of  data  and  is  helpful  in  interpreting 
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the  results.  Principal  components  analysis  often  reveals  relationships  that  were  not 
previously  suspected,  and  so  allows  interpretations  that  would  not  ordinarily  result 
(Johnson  and  Wichem  1992). 

Four  PCAs  were  conducted,  one  for  each  of  the  4  study  extents  (Figure  8,  Step 
P8).  The  resulting  principal  component  matrix  was  then  transformed  by  Varimax 
Rotation.  This  transformation  maximizes  the  correlation  between  the  PCs  and  original 
variables,  thus  facilitating  interpretation  (SAS  Institute  1988). 

The  second  statistical  analysis  method  used  was  stepwise  multiple  regression. 
Stepwise  multiple  regression  regresses  one  variable  on  a  set  of  variables  in  an  exploratory 
way,  to  obtain  a  minimum  of  unexplained  residual  variance  in  terms  of  the  smallest 
number  of  variables  from  the  data  set  (Osborne  1984).  The  bird  species  abundance  data 
were  regressed  against  the  calculated  principal  component  scores  resulting  from  the  PCA 
to  determine  if  kilometer-resolution  horizontal  spatial  pattern  metrics  are  suitable 
indicators  of  habitat  suitability  for  Neotropical  and  resident  songbirds.  A  separate 
multiple  regression  analysis  was  conducted  on  the  4  extents  by  regressing  each  bird 
species  on  the  PCs  (Figure  8,  Step  P9)  and  each  habitat  utilization  group  of  bird  species 
on  the  PCs.  This  resulted  in  a  total  of  72  multiple  regression  analyses. 

For  models  with  the  highest  R2  value,  logistic  regression  was  used  to  examine  the 
relationship  between  the  factors  and  the  probability  of  medium-high  species'  mean 
count/route.  In  this  analysis,  the  dependent  variable,  assumes  a  value  of  0  if  the  species 
was  on  a  route  with  low  or  accidental  abundance,  and  1  if  it  was  on  a  route  with  medium 


or  high  abundance. 


Climate  Data 
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An  analysis  of  background  climate  data  for  the  study  area  was  conducted.  Earthlnfo 
Environmental  Database  Summary  of  the  Day  climate  data  were  used  to  develop 
background  climate  data  for  the  study  timeframe  and  a  longer  period  of  record  timeframe 
(Earthlnfo  1998).  These  data  were  used  as  a  means  to  determine  if  the  weather  during 
the  study  timeframe  was  atypical.  This  information  is  listed  in  Appendix  B. 


Figure  6.  Diagram  of  a  full  and  partial  BBS  route. 
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Figure  7.  Diagram  of  a  BBS  route  and  the  4  physical  study  area  extents  formed  by 
buffering  the  entire  and  partial  routes  with  a  0.4  km  and  10  km  radius  buffer. 
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Figure  8.  Schematic  outline  of  sequence  of  procedures  used  to  assess  bird  and  landscape 
structure  associations. 
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Table  4.  Landscape  spatial  pattern  metrics  used  in  this  study. 


Metric 

Nomenclature 

Type  of  Metric 

Range  in  Value 

Calculation 

Patch  Type 
Calculation 
Based  On: 

Percent 
Agriculture  (% 
Ag) 

Composition 

(area) 

0  <  %Ag  <=  100 

Sum  of  the  areas  of  all  patches  of  the 
corresponding  patch  type,  divided  by 
total  landscape  area,  multiplied  by  100 

Cropland  and 
pasture 

Percent 
Deciduous 
Forest  (%  Dec) 

Composition 

(area) 

0  <  %Dec  <=  1 00 

Sum  of  the  areas  of  all  patches  of  the 
corresponding  patch  type,  divided  by 
total  landscape  area,  multiplied  by  100 

Deciduous 

forest 

Percent 
Evergreen 
Forest  (%  Ev) 

Composition 

(area) 

0  <  %Ev  <=100 

Sum  of  the  areas  of  all  patches  of  the 
corresponding  patch  type,  divided  by 
total  landscape  area,  multiplied  by  100 

Evergreen 

forest 

Percent  Mixed 
Forest  (%  Mx) 

Composition 

(area) 

0  <  %Mx  <=100 

Sum  of  the  areas  of  all  patches  of  the 
corresponding  patch  type,  divided  by 
total  landscape  area,  multiplied  by  100 

Mixed  forest 

Percent 
Forested 
Wetland 
(%  Fw) 

Composition 

(area) 

0  <  %Fw  <=  100 

Sum  of  the  areas  of  all  patches  of  the 
corresponding  patch  type,  divided  by 
total  landscape  area,  multiplied  by  100 

Forested 

wetland 

Mean  Patch 
Size  (MPS) 

Configuration 

(size) 

MPS  >  0,  without  limit.  Range 
limited  by  the  minimum  patch  size 
and  extent  of  image 

Sum  of  the  areas  of  all  patches  of  the 
corresponding  patch  type,  divided  by  the 
number  of  patches  of  the  same  type 

Forest 

Patch  Density 
(PD) 

Configuration 

(density) 

PD  >  0,  without  limit 

Number  of  patches  of  the  corresponding 
patch  type  divided  by  total  landscape 
area 

Forest 

Edge  Density 
(ED) 

Configuration 

(edge) 

ED  >=  0,  without  limit 

Sum  of  the  lengths  of  all  edge  segments 
involving  corresponding  patch  type, 
divided  by  the  total  landscape  area 

Forest 

Mean  Shape 
Index  (MSI) 

Configuration 

(shape) 

MSI  >=  1,  without  limit.  MSI  =  1 
when  all  patches  of  the 
corresponding  patch  type  are 
circular  MSI  increases  without  limit 
as  the  patch  shapes  become  more 
irregular 

Sum  of  the  patch  perimeter  divided  by 
the  square  root  of  patch  area  for  each 
patch  of  the  corresponding  patch  type, 
adjusted  by  a  constant  to  adjust  for  a 
circular  standard,  divided  by  the  number 
of  patches  of  the  same  type 

Forest 

Total  Core 

Area  Index 
(TCAI) 

Configuration 
composition 
(core  area) 

0  <=  TCAI  <  100  TCAI  =  0  when 
none  of  the  patches  of  the 
corresponding  patch  type  contain 
any  core  area  TCAI  approaches 

100  when  the  patches  of  the 
corresponding  patch  type,  because 
of  size,  shape,  and  edge  width, 
contain  mostly  core  area 

TCAI  equals  the  percentage  of  a  patch 
type  in  the  landscape  that  is  core  area 
based  on  a  250  m  (Robinson  et  al.  1995, 
Temple  1986)  specified  edge  width 

Forest 

Simpson’s 
Diversity  Index 
(SID!) 

Composition 

(diversity) 

0  <=  SIDI  <  1  SIDI  =  0  when  the 
landscape  contains  only  1  patch 
SIDI  approaches  1  as  the  number 
of  different  patch  types  (patch 
richness)  increases  and  the 
proportional  distribution  of  area 
among  patch  types  becomes  more 
equitable  Higher  value  the  greater 
likelihood  only  two  randomly  drawn 
patches  would  be  different  patch 
types  (greater  diversity) 

SIDI  equals  1  minus  the  sum,  across  all 
patch  types,  of  the  proportional 
abundance  of  each  patch  type  squared 

All  level  II 
patch  types 
simultaneously 

Interspersion 

and 

Juxtaposition 
Index  (IJI) 

Configuration 

(interspersion) 

0  <  IJI  <+  100  Approaches  0 
when  distribution  of  adjacencies 
among  unique  patch  types 
becomes  increasingly  uneven. 
Equals  100  when  all  patch  types 
are  equally  adjacent  to  all  other 
patch  types. 

IJI  equals  the  observed  interspersion 
over  the  maximum  possible  interspersion 
for  the  given  number  of  patch  types 

All  level  II 
patch  types 
simultaneously 
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Table  5.  Years  of  BBS  data  used  for  each  BBS  route  (from  1970-1976)  to  develop  the 
mean  yearly  abundance  for  each  species.  An  “x”  signifies  that  the  BBS  data 
were  available  for  that  year. 
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Table  5.  (contin’d)  Years  of  BBS  data  used  for  each  BBS  route  (from  1970-1976) 
to  develop  the  mean  yearly  abundance  for  each  species.  An  “x”  signifies 
that  the  BBS  data  were  available  for  that  year. 


Table  6.  Highest  mean  count/route  for  15  bird  species.  Developed  from  BBS 
data  from  1966-1985  (Hamel  1992). 


Species  Name 

AOUNo. 

Highest  mean  count/route 

Prairie  Warbler 

6730 

66 

Northern  Bobwhite 

2890 

277 

Field  Sparrow 

5630 

108 

Loggerhead  Shrike 

6220 

32 

Eastern  Kingbird 

4440 

36 

White-eyed  Vireo 

6310 

249 

Eastern  Wood-pewee 

4610 

27 

Yellow-throated  Vireo 

6280 

16 

Hooded  Warbler 

6840 

24 

Wood  Thrush 

7550 

53 

Yellow-billed  Cuckoo 

3870 

46 

Kentucky  Warbler 

6770 

20 

Brown-headed 

Nuthatch 

7290 

19 

Prothonotary  Warbler 

6370 

50 

Brown-headed 

Cowbird 

4950 

94 

Table  7.  Patch  types  in  Level  I  and  II  of  the  USGS  LULC  dataset. 
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LEVEL  1 

LEVEL  II 

Urban  or  built-up  land 

Residential 

Commercial  and  services 

Industrial 

Transportation,  communication,  utilities 

Industrial  and  commercial  complexes 

Mixed  urban  or  built-up  land 

Other  urban  or  built-up  land 

Agricultural  land 

Cropland  and  pasture 

Orchards,  groves,  vineyards,  nurseries  and  ornamental  horticultural  areas 

Confined  feeding  operations 

Other  agricultural  land 

Rangeland 

Herbaceous  rangeland 

Shrub  and  brush  rangeland 

Mixed  rangeland 

Forest  land 

Deciduous  forest  land 

Evergreen  forest  land 

Mixed  forest  land 

Water 

Streams  and  canals 

Lakes 

Reservoirs 

Bays  and  estuaries 

Wetland 

Forested  wetland 

Nonforested  wetland 

Barren  land 

Dry  salt  flats 

Beaches 

Sandy  areas  not  beaches 

Bare  exposed  rock 

Strip  mines,  quarries,  gravel  pits 

Transitional  areas 

Table  8.  Number  of  study  sites  in  the  four  extents  containing  deciduous  forest, 
evergreen  forest,  mixed  forest,  or  forested  wetland. 


Patch  Type 

Number  of  Study  Sites 

Extent  1 

Extent  2 

Extent  3 

Extent  4 

Deciduous  forest 

24 

43 

67 

183 

Evergreen  forest 

43 

53 

152 

260 

Mixed  forest 

52 

53 

234 

265 

Forested  wetland 

18 

40 

46 

163 

CHAPTER  IV 


RESULTS 

Introduction 

The  analysis  encompassed  four  discrete  steps:  extracting  abundance  values  for 
each  bird  species  from  each  landscape  unit  (BBS  route),  dividing  landscape  metrics  into 
appropriate  categories,  summarizing  landscape  metrics  into  PCA  axes  to  remove 
autocorrelation,  and  exploring  the  sensitivity  of  bird  abundance  to  the  various  PCA  axes. 

Extracting  Abundance  Values  for  Each  Bird  Species  from  Each  Landscape  Unit 

Each  BBS  route  is  approximately  40-km  long  and  consists  of  50  stops  spaced  0.8 
km  apart  (Keller  and  Scallan  1999).  The  records  for  each  BBS  route  included  the  number 
of  individuals  of  each  species  of  bird  observed  for  the  first  ten  stops,  second  ten  stops, 
third  ten  stops,  fourth  ten  stops,  and  fifth  ten  stops,  and  for  the  full  BBS  route  (Robbins  et 
al.  1986).  It  is  these  five  10-stop  summaries  or  5  mile  subsegments  that  form  the  basis  of 
the  5  partial  BBS  routes,  and  together,  the  full  BBS  route.  There  were  53  full  BBS  routes 
used  in  this  study.  And  with  each  of  the  53  full  BBS  routes  broken  into  5  partial  routes, 
this  resulted  in  265  partial  BBS  routes  in  the  study.  Since  the  partial  routes  provided  a 
finer  level  of  detail,  these  data  were  used  as  well  as  the  more  typical  full  BBS  route  data. 
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The  species  included  in  this  study  were  identified  as  conservation  priority 
species  in  1998.  The  data  that  were  used  in  this  study  are  based  on  BBS  data  from  1970- 
1976.  Even  though  the  species  had  not  yet  been  listed  as  conservation  priority  in  the 
1970-1976  timeframe,  the  average  of  the  53  mean  count  of  birds/route  for  each  species 
were  low,  with  the  exception  of  the  Northern  Bobwhite  (Figure  9).  The  more  recent 
timeframe  included  in  Figure  9  (1994-1998),  shows  that  the  mean  yearly  count  of 
birds/route  are  low,  with  the  Northern  Bobwhite  decreasing  substantially  from  the  1970- 
1976  timeframe.  Based  on  the  53  full  BBS  routes  the  mean  count  of  birds/route  for  each 
species  were  categorized  into  accidental,  low,  medium,  and  high  abundance  categories 
(Hamel  1992).  These  groupings  are  based  on  the  same  groupings  developed  by  Hamel 
and  used  for  range  maps  for  birds  across  the  South  (Hamel  1992).  The  groupings  are 
useful  in  determining  the  routes  that  have  high  or  medium  abundance  and  presumably  the 
most  suitable  habitat.  Based  on  Figure  10,  the  majority  of  BBS  routes  have  a  low  or 
accidental  abundance  for  1 1  of  the  species  under  study.  The  remaining  4  species  have  a 
medium  to  high  abundance  on  the  majority  of  routes. 

The  relative  abundance  data  from  the  partial  BBS  routes  was  also  examined. 

Since  there  is  no  accepted  method  for  categorizing  the  partial  routes  into  abundance 
categories,  a  table  is  used  to  show  the  distribution  of  each  species  into  groups  of  mean 
yearly  count  of  birds/partial  route  (Table  9).  For  six  of  the  bird  species,  there  were  zero 
counts  for  over  50  %  of  the  partial  routes.  For  all  bird  species  studied  (except  the 
Northern  Bobwhite),  well  over  50%  of  the  partial  routes  had  zero  counts  or  less  than  or 
equal  to  1  mean  yearly  count  of  birds/partial  route.  The  Northern  Bobwhite  was  the 
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exception,  since  it  had  a  much  higher  mean  yearly  count  of  birds/partial  route  than  the 
other  species. 


Dividing  Landscape  Metrics  Into  Appropriate  Categories 

The  landscape  structure  was  quantified  using  a  suite  of  12  spatial  metrics 
calculated  from  USGS  LULC  data  for  each  of  the  four  extents  around  each  BBS  route. 
As  used  in  this  report,  the  term  “Extent  1”  is  a  0.4  km  radius  buffer  of  a  full  BBS  route; 
“Extent  2”  is  a  10  km  radius  buffer  of  a  full  BBS  route;  “Extent  3”  is  a  0.4  km  radius 
buffer  of  each  partial  BBS  route;  and  “Extent  4”  is  a  10  km  radius  buffer  of  each  partial 
BBS  route.  The  extents  are  illustrated  in  Figure  7. 

The  variation  of  the  metric  values  for  the  4  extents  indicates  that  the  extent  of  the 
analysis  area  substantially  influences  some  landscape  spatial  pattern  metric  values  (Table 
10).  The  narrow  width  extents  resulted  in  mean  patch  sizes  and  core  area  values  that 
were  universally  smaller,  while  also  producing  patch  density  values  that  were 
substantially  larger,  and  edge  density,  and  interspersion  values  that  were  slightly  larger. 
There  were  less  noticeable  differences  in  the  MSI  and  SIDI  values.  These  findings  are 
also  illustrated  by  the  distribution  of  the  values  of  PD,  MPS,  ED,  MSI,  TCAI,  SIDI,  and 
DI  for  all  sites  for  the  4  extents,  which  are  shown  in  Figures  11-17,  respectively.  Again, 
the  pattern  of  difference  is  particularly  evident  in  the  distribution  of  values  for  PD,  MPS, 
TCAI  between  Extents  1  and  3  versus  Extents  2  and  4. 

The  composition  metrics  were  not  as  affected  by  spatial  extent.  Figure  18  shows 
that  the  cropland  and  pasture,  deciduous  forest,  evergreen  forest,  mixed  forest,  and 
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forested  wetland  patch  types  used  in  this  study  were  quite  similar  in  percent  for  the  4 
extents.  For  all  extents,  these  patch  types  comprised  approximately  90  %  of  the 
landscape. 

Summarizing  Landscape  Metrics  into  PCA  Axes  to  Remove  Autocorrelation 

A  principal  components  analysis  was  used  to  simplify  the  structure  of  the  spatial 
metric  data  sets  by  reducing  them  to  a  smaller  set  of  uncorrelated  variables  that  account 
for  a  large  part  of  the  variation  in  the  original  data  set.  Principal  components  analyses 
performed  to  describe  landscape  structure  produced  similar  components  for  the  4 
different  extents  (Table  1 1).  PCs  are  interpreted  based  on  the  pattern  and  strength  of  the 
loadings.  Based  on  the  component  loadings,  the  3  PCs  in  the  study,  were  universally 
interpreted  as  configuration  (PCI),  composition  (PC2),  and  diversity/interspersion  (PC3). 
These  PCs  provided  a  good  summary  of  all  the  original  variables. 

The  results  for  the  PC  analyses  for  the  4  extents  is  listed  in  table  1 1 .  The 
eigenvalue,  percent  variance,  and  cumulative  variance  for  each  component  for  each 
extent  are  listed  at  the  top.  The  eigenvalue  provides  a  means  of  identifying  the  number  of 
components  retained  in  each  analysis.  Principal  components  with  eigenvalues  >  1.0  with 
meaningful  interpretations  based  on  the  pattern  and  strength  of  the  loadings  were 
retained.  Three  components  were  retained  for  all  extents.  The  percent  variance  explains 
the  amount  of  variance  in  the  dataset  explained  by  each  principal  component.  The 
cumulative  percent  variance  simply  provides  a  cumulative  total  of  the  percent  variance 
explained  by  two  or  more  PCs.  PCI  accounted  for  about  35-37  percent  of  the  variance  in 
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the  landscape  datasets  for  Extents  1-4,  while  PCI  and  PC2  together  explained  about 
50-56  percent  of  the  variance,  and  all  three  PCs  combined  explained  about  65-74  percent 
of  the  variance  in  the  landscape  datasets. 

The  lower  part  of  Table  1 1  lists  the  component  loadings.  The  component 
loadings  are  the  correlations  between  the  principal  component  and  each  original  variable. 
The  interpretation  of  component  loadings  was  based  on  the  largest  loadings  for  each 
component.  Generally,  this  included  all  correlations  greater  than  0.5.  These  are  listed  in 
bold  type  in  the  table.  The  correlations  between  the  principal  component  and  each 
original  variable  are  either  negative  or  positive.  This  is  depicted  in  Table  1 1  with  a  or 
"+"  sign  in  front  of  the  loading. 

In  each  case,  PCI  had  high  component  loadings  for  MPS,  TCAI,  ED,  PD,  and 
MSI.  These  are  all  measures  that  help  to  describe  the  configuration  of  a  forested 
landscape.  Therefore,  PCI  was  interpreted  as  an  overall  configuration  component, 
contrasting  a  fewer  number  of  more  complex  shaped  large  forest  patches  with  large  core 
areas  with  a  relatively  small  amount  of  edge  in  largely  forested  landscapes  (positive 
loading),  with  many  small  more  simplistically  shaped  patches  with  small  core  areas,  with 
a  relatively  large  amount  of  edge,  in  mostly  agricultural  landscapes  (negative  loading). 
PC2  had  high  component  loadings  for  percent  mixed  forest,  evergreen  forest,  and 
forested  wetlands.  These  are  all  measures  that  describe  the  habitat  composition  of  a 
landscape.  Therefore,  PC2  was  interpreted  as  a  composition  component,  contrasting 
mostly  evergreen  or  forested  wetland  landscapes  (positive  loadings)  with  mixed  forest 
landscapes  (negative  loading).  That  these  variables  loaded  separately  on  their  own 
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principal  component,  indicates  that  the  habitat  composition  varies  independently  of  the 
configuration  of  the  forest  in  a  landscape.  Finally,  PC3  had  high  component  loadings  for 
SDDI,  UI,  and  deciduous  forest.  These  are  measures  that  describe  diversity  and 
interspersion  of  the  landscape.  Therefore,  PC3  was  interpreted  as  a 
diversity/interspersion  component,  contrasting  diverse  and  interspersed  landscapes 
(positive  loading)  with  landscapes  that  are  not  diverse  and  interspersed. 

Based  on  the  principal  components  analysis,  it  was  concluded  that  similar 
landscape  structure  components  exist  for  all  spatial  extents;  the  configuration  of  the  forest 
varies  independently  of  the  habitat  composition;  and  both  the  configuration  of  the  forest 
and  the  habitat  composition  vary  independently  of  the  diversity  and  interspersion  of  the 
landscape. 

Exploring  the  Sensitivity  of  Bird  Abundance  to  the  Various  PC  A  Axes 

How  the  landscape  structure  components  are  associated  with  the  observed 
abundance  of  the  conservation  priority  bird  species  was  examined  using  the  results  of  the 
PCA.  A  stepwise  multiple  regression  analysis  was  conducted  in  which  the  3  principal 
components  were  included  as  new  variables  to  model  their  associations  with  the 
abundance  of  each  bird  species.  The  mean  yearly  bird  count/route  or  partial  route  for 
each  of  the  fifteen  bird  species  and  the  habitat  utilization  bird  groups  was  regressed  on 
the  calculated  PCA  scores  for  the  three  PCs.  This  resulted  in  60  multiple  linear 
regression  models  for  the  individual  bird  species  and  12  for  the  habitat  utilization  bird 
groups.  The  results  of  the  stepwise  multiple  regression  analyses  indicated  which 
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components  showed  the  strongest  relationships  with  relative  abundance  of  bird  species. 
The  R2  values  represent  the  proportion  of  the  variation  in  relative  abundance  that  can  be 
explained  by  variation  in  the  components  included  in  the  model.  Based  on  the  4  extents, 
the  models  explained  3-50%  of  the  variation  in  abundance  of  the  individual  species.  The 
models  developed  for  Extent  1  explained  8-50%  of  the  variation  in  abundance  of  the 
individual  species;  those  developed  for  Extent  2  explained  14-39%;  those  developed  for 
Extent  3  explained  4  -  28%;  and  those  developed  for  Extent  4  explained  3-19%. 

A  summary  of  each  stepwise  multiple  regression  model  is  listed  in  Table  12.  The 
constant  and  coefficients  are  not  listed  in  Table  12,  rather  a  “+”  or  is  used.  A  "+" 
represents  a  positive  association  with  the  component,  while  a  represents  a  negative 
association  with  the  component.  The  partial  R2  values  are  fisted  for  each  component  in 
the  model,  as  well  as  the  R2  value  for  the  full  model.  If  there  is  no  data,  this  means  the 
PC  was  excluded  from  the  model  in  the  stepwise  multiple  regression  analysis.  All 
variables  included  in  the  model  have  a  P  <  0.05. 

In  spite  of  the  relatively  high  degree  of  unexplained  variance,  the  models 
generated  in  this  study  provide  a  methodology  for  determining  which  kilometer- 
resolution  horizontal  spatial  pattern  metrics  provide  indicators  of  habitat  suitability  for 
the  species  under  study.  These  models  are  not  suitable  for  predicting  bird  abundance, 
however;  the  methods  can  be  used  when  assessing  the  relative  importance  of  areas  for 
conservation  efforts  and  in  stratifying  the  areas  into  categories  of  suitability  for  the 
species  to  be  conserved  or  assessing  the  impacts  of  alternative  management  plans  that 
could  alter  or  remove  habitat  for  bird  species.  The  following  summary  of  results  pertains 
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to  the  extent  1  analysis.  The  overall  findings  from  extent  1  resulted  in  the  best  overall 
results  and  best  illustrate  the  relationship  between  bird  abundance  and  kilometer- 
resolution  horizontal  spatial  pattern  metrics. 

The  forest  configuration  and  habitat  composition  components  (PCI  and  PC2) 
were  found  to  be  the  most  frequent  significant  predictors  of  relative  abundance  for  the  15 
bird  species.  The  diversity/interspersion  component  (PC3)  was  a  significant  predictor  of 
abundance  for  only  2  birds.  Some  general  observations  can  be  made  about  the  way  that 
the  relative  abundance  of  individual  species  were  predicted  by  the  different  components. 
Specifically,  a  suite  of  species  that  included  the  White-eyed  Vireo,  Yellow-throated 
Vireo,  Hooded  Warbler,  Wood  Thrush,  Kentucky  Warbler,  Brown-headed  Nuthatch,  and 
the  Prothonotary  Warbler  were  positively  associated  with  the  configuration  component. 
This  appears  to  support  the  categorization  of  these  birds  as  forest  interior  to  forest 
interior/edge  species,  i.e.  preferring  relatively  fewer  patches  with  less  edge,  that  have  a 
relatively  large  forest  core  area  and  mean  forest  patch  size.  Another  set  of  species,  the 
Prairie  Warbler,  Northern  Bobwhite,  Field  Sparrow,  Loggerhead  Shrike,  and  the  Eastern 
Kingbird,  appear  to  be  edge  species.  These  birds  responded  negatively  to  the 
configuration  component,  inferring  they  prefer  edge  dominated  sites  with  many  patches 
that  have  a  smaller  mean  patch  size  and  core  area.  Not  surprisingly,  the  Brown-headed 
Cowbird  appears  to  be  a  generalist  in  that  it  was  not  affected  by  forest  configuration,  and 
did  not  seem  to  prefer  forest  interior  or  edge  type  habitats.  The  Yellow-Billed  Cuckoo 
was  also  not  sensitive  to  habitat  configuration,  a  finding  that  is  consistent  with  previous 
landscape  scale  research  (Howell  et  al.  2000).  In  general,  with  species  grouped  by 
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habitat  utilization,  those  species  grouped  as  interior  species  or  interior/edge  species 
responded  positively  to  the  configuration  component,  while  those  species  grouped  as 
edge  and  scrub  plus  field-edge  species  appear  to  be  respond  negatively  to  the 
configuration  component. 

A  suite  of  species,  including  the  Wood  Thrush,  Kentucky  Warbler,  Eastern 
Wood-Pewee,  Prairie  Warbler,  Brown-headed  Cowbird,  Yellow-billed  Cuckoo,  and  the 
Field  Sparrow  appear  to  be  negatively  associated  with  the  habitat  composition 
component.  This  implies  that  these  species  prefer  areas  with  relatively  large  amounts  of 
mixed  forest  and  small  amounts  of  evergreen  forest  or  forested  wetland.  The 
Prothonotary  Warbler  was  the  only  species  that  responded  positively  to  the  habitat 
composition  component,  preferring  areas  with  relatively  large  amount  of  forested  wetland 
or  evergreen  forest  and  small  amount  of  mixed  forest.  This  information  is  consistent  with 
that  developed  at  the  patch  level,  except  for  the  Field  Sparrow  and  Prairie  Warbler 
(Hamel  1992).  On  the  basis  of  the  results  of  this  analysis,  the  individual  species  were 
lumped  together  into  habitat  utilization  groups  for  another  stepwise  multiple  regression 
analysis. 

Habitat  Utilization  Groups 

The  interior  bird  species  had  an  R2  value  of  38  %  for  Extent  1 .  They  were 
positively  associated  with  forest  configuration  and  negatively  associated  with  habitat 
composition,  with  forest  configuration  explaining  about  28%  of  the  variance  and  habitat 
composition  explaining  about  10%  of  the  variance.  The  model  for  the  interior  bird 
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species  indicated  that  the  birds  were  more  likely  to  be  found  in  mixed  forest 
landscapes,  with  large  forest  patches  and  forest  core  areas,  and  a  fewer  number  of  patches 
with  more  complex  shapes  and  edge. 

The  interior/edge  bird  species  had  an  R2  value  of  42%  for  Extent  1 .  They  had  the 
same  associations  with  the  landscape  as  the  Interior  bird  species,  that  is  they  were 
positively  associated  with  the  forest  configuration  component  and  negatively  associated 
with  the  habitat  composition  component.  Forest  configuration  explained  about  20%  of 
the  variance  and  habitat  composition  explained  about  16%  of  the  variance.  The  model 
for  the  interior/edge  species  indicated  that  the  birds  were  more  likely  to  be  found  in 
mixed  forest  landscapes,  with  large  forest  patches  and  forest  core  areas,  and  a  fewer 
number  of  patches  with  more  complex  shapes  and  edge. 

The  edge  and  scrub,  field-edge  bird  species  had  an  R2  value  of  36%  for  Extent  1. 
They  had  a  negative  association  with  the  forest  configuration  component  and  a  negative 
association  with  the  habitat  composition  component.  Forest  configuration  explained 
about  26%  of  the  variance  and  habitat  composition  explained  about  10%  of  the  variance. 
The  model  for  the  edge  and  scrub,  field-edge  species  indicated  that  the  birds  were  more 
likely  to  be  found  in  mixed  forest  landscapes,  with  a  larger  number  of  smaller  forest 
patches  and  forest  core  areas  that  have  more  simplistic  shapes  and  more  edge. 

Individual  Species  with  Highest  R2  Values 

Out  of  the  15  species  studied,  five  species  had  models  with  R2  values  from  37- 
50%  (Table  12  and  13).  This  still  leaves  a  great  amount  of  unexplained  variance. 
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However,  given  that  existing  datasets  were  used  and  the  other  limitations  of  the  current 
study,  these  findings  show  promise  for  further  investigation.  The  five  species  with  the 
highest  R2  values  were  the  Wood  Thrush,  Hooded  Warbler,  White  Eyed  Vireo,  Kentucky 
Warbler,  and  Prothonotary  Warbler.  The  10  remaining  species  had  lower,  highly  variable 
R2  values.  For  these  10  cases,  the  large  amount  of  unexplained  variance  limits  a 
meaningful  discussion  of  the  models  for  those  particular  species.  Further  examination  of 
the  other  5  species  with  relatively  higher  R2  values  is  informative.  The  form  of  the  model 
is  shown  below  as  Equation  9.  Table  13  lists  the  coefficients  for  the  statistical  models  for 
each  of  these  species.  The  habitat  suitability  needs  for  these  bird  species  were  also 
assessed  using  logistic  regression. 

Abundance  =  C0  +  Cx  (PCI)  +  C2  (PC2)  +  C3  (PC3)  (9) 

The  statistical  models  for  the  individual  species  are  best  understood  from  the  perspective 
of  each  species'  natural  history  and  habitat  preferences. 


Wood  Thrush 

The  model  for  the  Wood  Thrush  had  an  R2  =  50%.  In  the  multiple  regression 

2 

analysis  in  this  study,  the  Wood  Thrush  responded  primarily  to  habitat  composition  (pR 
=  44%)  and  to  a  much  lesser  degree  forest  configuration  (pR2  =  6%).  In  the  logistic 
regression  analysis,  the  Wood  Thrush  responded  significantly  only  to  the  habitat 
composition  of  the  landscape.  As  the  habitat  increased  in  mixed  forest  and  decreased  in 
evergreen  and  forested  wetlands,  the  probability  of  medium  to  high  Wood  Thrush 
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abundance  increased  (Figure  19).  It  appears,  from  this  study,  that  the  Wood  Thrush 
may  not  be  affected  so  much  by  the  core  area  and  size  of  the  forest  patches,  assuming  that 
the  forest  fragments  are  large  enough  to  include  an  average  territory  size.  These 
landscape  level  findings  are  consistent  with  previous  research  conducted  at  the  patch 
level.  In  previous  research  it  was  found  that  forest  area  was  a  significant  predictor  of 
relative  abundance  for  the  Wood  Thrush,  with  the  predicted  probability  of  occurrence 
increasing  as  the  area  of  forest  increased  (Robbins  et  al.  1989).  The  Wood  Thrush  is  a 
ground-foraging  Neotropical  migrant  that,  based  on  studies  at  the  patch  level,  prefers 
interior  and  edges  of  deciduous  and  mixed  forest,  especially  well-developed,  upland, 
mesic  ones  (James  et  al.  1984).  In  contrast,  in  a  recent  study  in  which  both  vegetation 
variables  and  landscape  variables  were  used,  the  Wood  Thrush  responded  only  to 
multiple  vegetation  variables  (positive  association  with  stems  less  than  2  cm,  logs,  and 
stems  greater  than  50  cm),  suggesting  a  requirement  for  second  growth  forest  or  a 
developed  understory  within  a  mature  forest  setting.  It  did  not  respond  to  the  landscape 
variables  (MPS,  TCAI,  ED,  %Forest)  (Howell  et  al.  2000).  This  apparent  disparity  could 
be  related  to  the  10-kilometer  radius  extent  on  which  Howell  calculated  the  landscape 
metrics.  Since  the  home  range  for  the  Wood  Thrush  is  documented  as  .15  kilometers 
(Roth  et  al.  1996),  the  10  kilometer  extent  appears  to  be  excessive.  In  this  dissertation 
research,  the  best  findings  for  the  Wood  Thrush  were  based  on  extent  1  (.4  km  radius 
buffer  on  the  entire  BBS  route).  In  Robbins  et  al.  (1989),  the  extent  on  which  the 
landscape  variables  were  calculated  was  a  2  km  radius  buffer.  It  appears  in  Howell  et  al. 
(2000),  that  the  vegetation  variables,  which  were  calculated  on  a  5  meter  radius,  are  more 
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important  in  explaining  the  abundance  of  the  bird  when  the  landscape  variables  were 
calculated  on  an  excessive  extent.  This  does  not  preclude  the  association  between  the 
bird  abundance  and  landscape  variables. 

Hooded  Warbler 

The  model  for  the  Hooded  Warbler  had  an  R2  of  45%.  In  this  study,  the  Hooded 
Warbler  responded  positively,  at  a  landscape  level,  primarily  to  forest  configuration. 
Forest  configuration  explained  about  35%  of  the  variance  and  diversity/interspersion 
explained  about  10%  of  the  variance.  In  the  logistic  regression  analysis,  the  Hooded 
Warbler  was  significantly  affected  by  the  forest  configuration  and  diversity/interspersion 
components.  The  probability  of  finding  medium  to  high  abundance  of  the  Hooded 
Warbler  increased  with  increasing  forest  patch  size  and  core  area  size  and  decreasing 
patch  numbers  and  edge(Figure  20),  and  increased  with  lower  landscape  diversity  and 
interspersion  and  percent  deciduous  forest(Figure  21).  These  landscape  level  findings  are 
consistent  with  the  existing  patch  level  knowledge  about  the  species.  For  example,  the 
Hooded  Warbler  is  a  small  migratory  songbird  that,  in  its  breeding  range,  inhabits 
deciduous  and  mixed  hardwood  forests.  It  favors  moist  forests  with  a  fairly  dense 
understory  (Ogden  and  Stutchbury  1994).  It  is  considered  a  forest-interior  species 
because  it  is  restricted  to  larger  woodlots  (Hamel  1992).  There  is  one  conflict  with 
existing  patch  level  knowledge,  in  that  the  Hooded  Warbler  is  known,  from  a  patch 
habitat  scale,  to  prefer  deciduous  forest,  but  as  a  part  of  the  diversity/interspersion 
component,  the  bird  is  shown  not  to  prefer  deciduous  forest.  Since  the 
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diversity/interspersion  component  is  confounded  with  the  percent  deciduous  forest 
variables,  it  is  difficult  to  determine  if  the  association  with  this  component  is  based 
primarily  on  the  diversity  and  interspersion  of  the  landscape  or  the  percent  of  the 
landscape  composed  of  deciduous  forest.  Nevertheless,  by  far  the  most  important 
landscape  variable  to  the  Hooded  Warbler  is  the  overall  configuration  of  the  forest,  not 
the  landscape  diversity  and  interspersion. 

White-Eved  Vireo 

The  White-Eyed  Vireo  model  had  an  R2  of  41%.  In  this  study,  the  White-Eyed 
Vireo  responded,  at  a  landscape  level,  about  equally  to  forest  configuration  and  landscape 
diversity  and  interspersion.  Forest  configuration  explained  about  20  %  of  the  variance 
and  diversity/interspersion  explained  about  21%  of  the  variance.  In  the  logistic 
regression  analysis,  the  White-Eyed-Vireo  was  significantly  affected  by  forest 
configuration  and  diversity/interspersion  components.  The  probability  of  finding 
medium  to  high  abundance  of  the  White-Eyed-Vireo  increased  with  increasing  forest 
patch  size  and  core  area  size  and  decreasing  number  of  patches  and  edge  (Figure  20),  and 
increased  with  lower  landscape  diversity  and  interspersion  (Figure  21).  It  is  surprising 
that  there  was  a  statistically  significant  positive  relationship  with  the  forest  configuration 
component,  since  the  bird  is  not  considered,  at  a  patch  level,  to  be  a  forest  interior  bird. 

In  its  breeding  range,  it  is  common  in  its  preferred  habitat  of  dense  secondary  deciduous 
scrub,  streamside  thickets,  wood  margins,  and  overgrown  pastures  (Hamel  1992,  Graber 
et  al.  1985).  It  would  be  expected  that  the  bird  would  have  had  a  negative  relationship 
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with  the  forest  configuration  component.  The  White-Eyed  Vireo  was  only  one  of  two 
birds  that  had  an  association  with  the  diversity  and  interspersion  of  the  landscape. 

There  are  several  possible  explanations  for  the  apparently  contradictory  findings 
for  the  White-Eyed  Vireo.  The  most  obvious  problem  is  that  the  USGS  LULC  dataset 
lacks  a  habitat  category  that  coincides  with  that  which  the  bird  has  been  found  to  prefer  at 
the  local  scale.  The  bird  prefers  a  habitat  of  secondary  deciduous  scrub,  overgrown 
pastures  and  abandoned  farmland,  or  wood  margins,  and  these  categories  are  not  depicted 
on  the  USGS  LULC  dataset.  Therefore,  it  is  not  surprising  that  the  results  are 
inconsistent  with  existing  local  habitat  knowledge.  Also,  since  the  home  range  of  this 
bird  is  not  known,  the  extent  to  which  the  landscape  metrics  were  calculated  could  be 
either  too  small  or  too  large.  Lastly,  since  the  BBS  is  a  roadside  survey  that  has  several 
potential  sources  of  bias,  the  bird  counts  for  this  bird  could  also  be  contributing  to  the 
apparently  contradictory  results. 

Prothonotarv  Warbler 

The  Prothonotary  Warbler  model  had  an  R2  of  37%.  From  a  landscape  context, 
the  Prothonotary  Warbler  was  most  associated  with  the  habitat  composition  component, 
and  to  a  lesser  degree  the  forest  configuration  component.  Habitat  composition  explained 
about  24%  of  the  variance  with  forest  configuration  explaining  about  13%  of  the 
variance.  In  the  logistic  regression  analysis,  the  Prothonotary  Warbler  was  significantly 
affected  by  habitat  composition  and  forest  configuration.  The  probability  of  finding 
medium  to  high  abundance  of  the  Prothonotary  Warbler  increased  significantly  with 
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decreasing  mixed  forest  cover  and  increasing  evergreen  and  forested  wetland  cover 
(Figure  19),  and  increased  with  increasing  forest  patch  size  and  core  area  size  and 
decreasing  patch  numbers  and  edge  (Figure  20).  These  landscape  level  findings  are 
consistent  with  findings  at  the  patch  level,  in  that  in  breeding  areas,  it  inhabits 
bottomland  hardwood  forests  and  other  forested  wetlands.  It  exhibits  area  sensitivity, 
avoiding  forests  <  100  ha  in  area  (Kahl  et  al.  1985).  Other  key  features  of  its  breeding 
habitat  are  presence  of  water  near  wooded  areas  with  suitable  cavity  nest  sites,  low 
elevation,  flat  terrain,  shaded  forest  habitats  with  sparse  understory,  and  in  some  places, 
presence  of  bald  cypress,  all  of  which  were  not  included  in  this  study  (Kahl  et  al.  1985, 
Petit  1999,  Robbins  et  al.  1989). 

Kentucky  Warbler 

The  model  for  the  Kentucky  Warbler  had  an  R2  of  38%.  In  this  study,  the 
Kentucky  Warbler  responded,  at  a  landscape  level,  primarily  to  habitat  composition  and 
to  a  much  lesser  degree  forest  configuration.  The  habitat  composition  component 
explained  the  greater  part  of  the  variance  at  31%,  with  forest  configuration  explaining 
about  7%.  In  the  logistic  regression  analysis,  the  Kentucky  Warbler  was  significantly 
affected  only  by  forest  composition.  The  probability  of  finding  medium  to  high 
abundance  of  the  Kentucky  Warbler  increased  significantly  with  increasing  mixed  forest 
cover  and  decreasing  evergreen  and  forested  wetland  cover  (Figure  19).  In  previous 
research  conducted  at  the  patch  level,  forest  area  was  a  significant  predictor  of  probability 
of  occurrence  of  the  Kentucky  Warbler  (Robbins  et  al.  1989).  It  is  commonly  found  in 


rich,  moist,  deciduous  forests  in  the  southeastern  United  States  and  is  rarely  observed 
in  agricultural  habitats  (McDonald  1998)  or  in  conifers  (Hamel  1992).  Also,  a  study  in 
Missouri  indicated  that  large  blocks  of  suitable  habitat  (>  500  ha)  are  necessary  for 
successful  breeding  (Gibbs  and  Faaborg  1990)  while  the  minimum  tract  size  has  also 
been  listed  as  45  ha  (Hamel  1992).  Although  the  species  may  have  minimum  are 
requirements,  it  does  not  necessarily  have  to  be  present  in  1  contiguous  patch. 
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Figure  9.  Mean  yearly  abundance  of  each  bird  species  across  all  routes. 


Figure  10.  Distribution  of  birds  in  terms  of  high,  medium,  low,  and  accidental  abundance. 
Abundance  categories  determined  by  first  determining  the  highest  mean  count/route  for 
each  bird  species  across  its  range  (Table  6)  (Hamel  1992).Routes  with  at  least  30  %  of 
the  peak  value  are  considered  high,  10%  are  considered  medium,  5%  of  peak  value  are 
considered  low,  and  less  than  5%  are  considered  accidental  abundance  (Hamel  1992). 
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Figure  12.  D: 
Extents  1-4. 
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Upper  Value  of  Range 
Total  Core  Area  Index  (TCAI) 


Figure  15.  Distribution  of  the  value  for  total  core  area  index  for  each  study  site  for 
Extents  1-4. 


Upper  Value  of  Range 
Simpson's  Diversity  Index  (SID I) 

Figure  16.  Distribution  of  the  value  for  Simpson's  diversity  index  for  each  study  site  for 
Extents  1-4. 
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Principal  Component  2  -  Habitat  Composition 

Figure  19.  Probability  of  medium  to  high  abundance,  based  on  logistic  regression,  for  the 
Kentucky  Warbler,  Prothonotary  Warbler,  and  Wood  Thrush,  in  relation  to  habitat 
composition  (PC2). 


Figure  20.  Probability  of  medium  to  high  abundance,  based  on  logistic  regression,  for  the 
Hooded  Warbler,  White-eyed  Vireo,  and  Prothonotary  Warbler,  in  relation  to  forest 
configuration  (PCI). 


77 


Figure  21.  Probability  of  medium  to  high  abundance,  based  on  logistic  regression,  for 
the  Hooded  Warbler  and  White-eyed  Vireo,  in  relation  to  landscape 
diversity/interspersion  (PC3). 


Table  9.  Frequency  of  the  mean  yearly  count  of  birds/route  for  each  species  over  all  partial  BBS  routes. 
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Species 

Name 

Prairie  Warbler 

Northern  Bobwhite 

Field  Sparrow 

Loggerhead  Shrike 

Eastern  Kingbird 

White-eyed  Vireo 

Eastern  Wood-pewee 

Yellow-throated  Vireo 

Hooded  Warbler 

Wood  Thrush 

Yellow-billed  Cuckoo 

Kentucky  Warbler 

Brown-headed 

Nuthatch 

Prothonotary  Warbler 

Brown-headed 

Cowbird 

Table  11.  Results  of  the  principal  components  analyses  for  the  4  extents.  Correlations  greater  than  0.5  are  listed  in  bold 
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Table  12.  Summary  of  multiple  linear  regression  models  for  each  bird  species  for  each  extent 


Tablel3.  Coefficients  for  the  statistical  models  for  each  of  the  five  species 
with  R2  values  greater  than  36%. 


Species  Name 

AOU 

C0  (Constant) 

O 

p 

C2  (PC2) 

C3  (PC3) 

Model 

White-eyed  Vireo 

6310 

7.46 

2.72 

0.00 

-2.78 

41% 

Hooded  Warbler 

6840 

1.43 

0.99 

0.00 

-0.52 

45% 

Wood  Thrush 

7550 

9.16 

1.76 

-4.27 

0.00 

50% 

Kentucky  Warbler 

6770 

1.43 

0.43 

-0.88 

0.00 

38% 

Prothonotary  Warbler 

6370 

1.86 

1.02 

1.42 

0.00 

37% 

Interior  species 

2.87 

1.42 

-0.86 

0.00 

38% 

Interior/Edge  species 

33.46 

7.68 

-6.75 

-4.38 

42% 

Edge  and  scrub,  field- 
edge  species 

71.93 

-15.93 

-9.93 

0.00 

36% 

CHAPTER  V 


DISCUSSION 

Introduction 

The  intent  of  this  study  was  to  determine  if  landscape  spatial  patterns  from  broad  scale 
data  are  useful  indicators  of  habitat  suitability  for  the  species  under  study.  Resulting 
models  for  five  of  the  bird  species  correctly  classified  between  one-third  and  one-half  of 
the  land  area.  The  unexplained  variance  (>50%)  of  these  models  probably  reflects  the 
limitations  of  this  study  and  the  data  sources  used.  Nevertheless,  the  statistically 
significant  study  findings  should  make  it  clear  that,  in  addition  to  the  habitat  patch 
preferences  and  natural  history  traits  of  the  birds,  managers  must  consider  landscape 
features  in  conservation  and  management  activities.  Applied  over  broad  areas,  the 
methodology  and  models  can  substantially  contribute  to  conservation  efforts. 

Spatial  Pattern  Metrics  and  Habitat  Suitability  Assessment 

This  study  has  challenged  the  common  view  that  BBS  data  can  only  be  used  for 

trend  analysis.  The  results  of  the  study  show  that  there  are  significant  benefits  to  utilizing 

BBS  and  USGS  LULC  data  to  determine  which  kilometer-resolution  horizontal  spatial 

pattern  metrics  are  suitable  indicators  of  habitat  suitability  for  conservation  birds.  The 

bird-to-landscape  habitat  associations  found  in  this  study  are  generally  consistent  with 

patch  level  habitat  preferences.  The  models  generated  have  low  to  moderate  R2  values 
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and  hence  may  not  be  appropriate  for  predicting  abundance.  However,  they  are 
appropriate  for  assessing  broad  scale  habitat  suitability  and  landscape  change  as  reasons 
for  population  declines,  especially  for  the  five  species  with  R2  >  35%  (W ood  Thrush, 
Hooded  Warbler,  White-Eyed  Vireo,  Prothonotary  Warbler,  and  Kentucky  Warbler). 
This  incorporation  of  spatial  pattern  metrics  from  broad  scale  data  should  be  quite  useful 
in  habitat  suitability  assessment. 

Incorporation  of  spatial  pattern  metrics  into  a  habitat  suitability  assessment  is 
premised  upon  the  generation  of  appropriate  metrics.  This  research  is  based  on  the 
assertion  that  including  a  comprehensive  set  of  landscape  metrics  is  the  best  way  to 
produce  an  all-encompassing  quantitative  description  of  the  landscape.  Landscape 
metrics,  such  as  composition,  patch  size,  edge  density,  patch  shape,  and  core  area,  are 
highly  correlated,  and  yet  measure  significantly  different  aspects  of  the  landscape 
(Forman  and  Gordon  1986).  The  individual  importance  of  these  variables  leads  to  the 
theory  that  a  landscape  is  best  described,  not  by  its  parts,  but  as  a  comprehensive  whole. 
For  example,  misleading  findings  may  result  if  one  includes  only  the  percent  forest  cover 
metric  in  the  analysis  of  two  landscapes  with  equivalent  amounts  of  forest  cover.  The 
two  landscapes  could  differ  in  all  other  landscape  measures  (e.g.  patch  size,  edge  density, 
patch  shape,  core  area)  due  to  the  patch  arrangement.  Without  other  metrics  that  would 
provide  this  additional  information,  the  use  of  the  percent  forest  metric  could  result  in 
misleading  findings.  This  contention  is  supported  by  the  findings  of  Rosenberg  et  al. 
(1999),  that  inclusion  of  only  a  single  variable,  such  as  mean  forest  patch  size,  may  yield 
misleading  conclusions  about  how  birds  may  respond  to  habitat  changes.  Rosenberg  et 
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al.  (1999)  and  others  have  also  been  successful  in  describing  the  landscape  and  finding 
bird-to-landscape  associations  with  a  comprehensive  set  of  metrics  (McGarigal  and 
McComb  1995,  Osborne  1984,  Roberts  and  Norment  1999).  In  related  research, 
Mladenoff  et  al.  (1997)  was  also  most  successful  in  using  a  combination  of  landscape 
pattern  metrics  when  differentiating  landscape  pattern  from  USGS  LULC  data.  The 
argument  for  the  utilization  of  a  comprehensive  set  of  metrics  is  not  necessarily 
weakened  by  the  apparent  success  of  other  studies  that  have  focused  only  on  individual 
metrics  or  a  small  subset  of  metrics  with  relatively  strong  results.  For  example,  Kruess 
and  Tschamtke  (1994)  used  only  an  isolation  measure  to  study  the  fragmentation  of 
habitats  in  the  agricultural  landscape  and  its  effect  on  the  number  of  parasitized  insect 
species  (R2  =  69%).  Others  have  been  similarly  successful  (R2  >  40%)  using  only  one  or 
two  metrics  (Temple  1986,  McIntyre  1995,  Galli  et  al.  1976,  Robinson  1995,  Van  Dorp 
and  Opdam  1987).  Since  these  studies  have  focused  on  a  single  metric  or  a  few  metrics, 
it  is  impossible  to  say  if  the  results  could  have  been  improved  with  a  more 
comprehensive  set  of  metrics. 

Inclusion  of  a  comprehensive  set  of  landscape  metrics  is  complicated  by  their 
high  degree  of  correlation.  In  many  instances,  the  search  for  a  metric  or  metrics  that 
provide  the  best  results  in  bird-to-landscape  habitat  associations  has  centered  around  the 
fact  that  many  of  the  landscape  metrics  are  highly  correlated  and  should  not  be  used 
together  in  multiple  regression  analyses  because  of  their  multicollinearity.  However, 
with  the  use  of  PCA,  the  interrelationships  among  a  large  number  of  variables  can  be 
explained  in  terms  of  their  common  underlying,  uncorrelated  principal  components  (Hair 
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et  al.  1992).  This  is  a  powerful  concept  that  allows  the  inclusion  of  the  influence  of  all 
these  important  variables  in  spite  of  their  high  degree  of  correlation.  The  uncorrelated 
principal  components  found  in  the  current  study  and  in  Rosenberg  et  al.  (1999)  and 
McGarigal  and  McComb  (1995)  imply  that  the  components  measure  different 
'dimensions'  in  the  data  (Manly  1990).  In  the  current  study,  the  'dimensions'  included 
habitat  composition,  forest  configuration,  and  landscape  diversity/interspersion.  Had 
only  a  single  variable  been  assessed,  this  comprehensive  description  of  the  landscape 
would  be  replaced  by  a  narrower  and  perhaps  misleading  viewpoint.  From  the  PCA 
(Table  11),  the  component  loadings  or  correlations  of  the  components  and  the  original 
variables  show  that,  in  the  forest  configuration  component,  the  MPS,  TCAI,  PD,  MSI, 
and  percent  agriculture  variables  all  have  about  equal  loadings.  This  shows  the  relative 
equal  importance  of  all  the  metrics  in  providing  an  overall  description  of  the  landscape. 
Figure  22,  provides  an  example  of  how  misleading  it  can  be  if  only  individual  metrics  are 
used  to  describe  a  landscape.  Figure  22  shows  two  landscapes  with  the  same  MPS,  but 
with  landscape  B  representing  a  more  fragmented  landscape  than  landscape  A.  The 
inclusion  of  the  other  forest  configuration  metrics  provides  a  better  overall  description  of 
the  landscape. 

The  analysis  of  spatial  pattern  metrics  led  to  several  generalizations  that  may  be 
useful  in  addressing  habitat  suitability  issues  at  a  broad  scale.  For  the  five  best  species 
models,  forest  configuration  was  a  significant  variable  in  all  5  models,  landscape 
composition  was  significant  in  3  of  the  models  and  diversity/interspersion  was  significant 
in  only  2  of  the  models.  When  composition  and  configuration  were  the  components  in 
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the  species  models,  composition  always  explained  a  greater  part  of  the  variance. 
Typically,  diversity/interspersion  explained  the  smallest  part  of  the  variance  in  the 
models.  Clearly,  considerable  focus  should  be  on  composition  and  configuration  when 
addressing  habitat  suitability  issues. 

Some  species  may  require  further  study  to  clarify  bird-to-landscape  associations. 
For  example,  the  White-Eyed  Vireo  model  showed  a  preference  for  landscapes  with  large 
forest  patches  with  large  core  areas,  less  edge  and  more  naturally  shaped  patches.  This  is 
similar  to  what  forest  interior  species,  such  as  the  Hooded  Warbler,  were  found  to  prefer. 
This  is  contrary  to  what  was  expected  since  the  White-Eyed  Vireo  is  widely  known  to 
prefer  dense  secondary  deciduous  scrub,  streamside  thickets,  wood  margins,  and 
overgrown  pastures  (Hamel  1992,  Graber  et  al.  1985),  and  forest  dominated  habitats  yet 
with  a  large  component  of  scrub/shrub  (personal  communication  Hunter  2000).  However, 
since  the  preferred  habitat  types  of  the  White-Eyed  Vireo  were  not  specifically  depicted 
in  the  USGS  LULC  data  (e.g.  scrub/shrub),  these  findings  may  be  misleading  and  should 
be  further  studied  with  data  that  address  the  species  habitat  needs.  This  also  applies  to 
other  species,  such  as  the  Prairie  Warbler,  Northern  Bobwhite,  and  Field  Sparrow,  that 
prefer  open  grassy  stands,  and  the  Loggerhead  Shrike  and  Eastern  Kingbird,  that  prefer  a 
more  open  savanna  type  habitat.  In  this  study,  these  species  showed  a  preference  for 
landscapes  with  a  high  percentage  of  agricultural  land  and  smaller  forest  patches  with 
less  core  area  with  more  edge  and  unnatural  man  induced  shapes.  Again,  these  findings 
are  probably  the  result  of  the  lack  of  scrub/shrub  and  grassland  habitat  classes  for  these 
species,  thus  limiting  the  bird-to-landscape  associations  for  these  birds.  Future  research 
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should  address  these  edge  and  scrub,  and  field-edge  species  with  pertinent  data.  One 
example  of  a  dataset  that  could  be  used  includes  the  USGS  National  Land  Cover  Data  set 
which  was  developed  from  30  m  Thematic  Mapper  (TM)  data.  This  dataset  was 
developed  in  the  1990's  and  includes  classifications  for  categories  such  as  shrubland, 
grasslands/herbaceous,  as  well  as  specific  forest  and  agricultural  categories.  Future 
research  would  provide  a  comparison  of  the  applicability  of  the  finer  grained  datasets  for 
the  purposes  of  habitat  suitability  assessment. 

Conservation  Implications 

It  has  been  observed  that  conservation  and  resource  management  should  be  based 
at  least  as  much  on  science  as  on  politics  and  economics.  It  has  also  been  noted  that  it  is 
important  to  continue  to  study  bird  populations  in  order  to  conserve  them,  as  information 
on  how  populations  respond  to  landscape  change  is  available  for  very  few  species  (Wiens 
1994).  The  current  research  is  part  of  an  effort  to  address  this  need.  The  results  of  the 
current  study,  as  well  as  others  studies  that  have  focused  on  the  landscape  scale,  have 
varied  dramatically  between  species,  implying  that  separate  conservation  guidelines  for 
each  species  must  be  generated  (Rosenberg  et  al.  1999). 

This  research  and  others  (e.g.  Donovan  et  al.  1997,  Howell  et  al.  2000,  Dettmers 
and  Bart  1999)  supports  the  position  that  a  good  first  approximation  for  assessing  habitat 
suitability  for  certain  species  is  to  assess  habitat  characteristics  at  the  landscape  scale. 
With  the  use  of  GIS  and  remotely  sensed  data,  relatively  inexpensive  mapping 
applications  can  be  developed  at  the  landscape  scale.  Assessment  of  the  amount  and 
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configuration  of  habitat  over  large  areas,  in  a  manner  that  is  neither  labor  intensive  nor 
prohibitively  time-consuming,  holds  many  benefits  for  mangers  and  researchers  working 
on  large-scale  issues  (Dettmers  and  Bart  1999).  This  is  fortunate,  since  detailed  data  on 
plant  species  composition  or  structure  is  seldom  available  for  large  areas. 

A  first  approximation  assessment  of  habitat  suitability  can  be  developed  at  the 
landscape  scale  by  using  spatial  pattern  metrics  developed  from  kilometer  scale 
landcover  data.  This  is  useful  in  determining  the  relative  importance  of  areas  for 
conservation  efforts  and  in  assessing  the  impact  of  alternative  management  plans  that 
could  alter  or  remove  habitat  for  bird  species.  For  example,  results  of  the  current  study 
and  others  show  that,  at  the  landscape  scale,  the  Wood  Thrush  was  primarily  associated 
with  habitat  composition  and  forest  configuration  (Fauth  1997,  Robbins  et  al.  1989). 
Therefore,  as  a  first  step  in  assessing  suitable  habitat  for  the  Wood  Thrush,  habitat 
composition  and  forest  configuration  metrics  developed  from  broad  scale  data  can  be 
examined  for  a  regional  area.  Pertinent  sites  from  the  broad  scale  first  approximation 
assessment  could  be  further  assessed  in  terms  of  habitat  vegetation  and  structure 
variables.  For  example,  Wood  Thrush  abundance  is  known  to  be  sensitive  to  habitat 
vegetation  and  structure  variables,  such  as  stems  less  than  2  cm,  number  of  logs,  and 
stems  greater  than  50  cm  (Howell  et  al.  2000). 

This  study  provides  the  basis  for  an  improved  understanding  of  the  influence  of 
landscape  structure  on  bird  species  population  changes.  Past  attempts  to  explain 
population  declines  of  migrants  detected  by  BBS  data  have  been  complicated  by  a  variety 
of  factors,  including  a  lack  of  information  on  existing  vegetation  and  on  possible  long- 
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term  changes  in  vegetation  (James  et  al.  1992).  The  methods  and  models  developed  in 
the  current  study  are  useful  in  determining  the  role  of  habitat  fragmentation  in  causing 
population  declines  shown  in  the  BBS  data.  The  models  generated  for  the  forest  interior 
and  forest  interior/edge  species  (R2>  35%)  may  be  directly  applicable  to  determining  the 
degree  to  which  landuse  changes  around  BBS  routes  affect  habitat  suitability  and  may 
offer  partial  explanations  for  the  long-term  population  changes  in  these  species.  It  is 
widely  recognized  that  population  declines  are  at  least  partially  attributable  to  habitat 
fragmentation.  Findings  from  the  current  study,  as  well  as  other  recent  studies  (e.g. 
Rosenberg  et  a.  1999,  Donovan  et  al.  1997,  Howell  et  al.  2000,  McGarigal  and  McComb 
1995,  Bolger  et  al.  1997)  provide  empirical  evidence  that  habitat  suitability,  as  quantified 
by  landscape  metrics,  is  associated  with  bird  species  abundance.  Future  research  of  this 
type  could  help  refine  conservation  guidelines  to  address  how  fragmentation  of  the 
landscape  affects  bird  species  populations. 

Limitations  of  the  Study  Approach 

A  potential  limitation  of  this  study  is  the  quality  of  the  existing  BBS  data.  A 
potential  bias  is  introduced  in  the  data  set  by  the  fact  that  the  survey  process  most 
frequently  counts  species  characteristic  of  roadsides,  with  the  species  that  sing  loudest 
and  most  frequently  being  the  easiest  to  detect.  Also,  volunteers  vary  in  their  ability  to 
hear,  identify,  and  estimate  the  abundance  of  birds  (Sauer  and  Droege  1990).  In  the 
current  study,  this  makes  it  difficult  to  determine  if  there  is  indeed  little  or  no  association 
in  the  models  with  low  R2  values,  or  if  the  species  is  just  not  well  documented  in  the  BBS 
data.  A  further  criticism  of  the  survey  is  that  the  BBS  is  a  roadside  survey,  and  habitat 
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changes  along  roadsides  may  not  be  representative  of  regional  habitat  changes. 
Additionally,  many  habitats  are  not  well  covered,  and  the  species  that  specialize  in  those 
habitats,  such  as  wetland  birds,  are  poorly  sampled. 

It  was  initially  thought  that  the  partial  BBS  route  data  might  provide  a  finer  level 
of  detail  and  perhaps  improve  the  analytical  results  in  this  study.  However,  the  lack  of 
bird  counts  on  the  partial  routes  for  most  of  the  species  seemed  to  confound  rather  than 
improve  the  results.  In  fact,  60%  of  the  birds  studied  were  not  counted  at  all  on  at  least 
30%  of  the  partial  routes.  Also,  the  lack  of  an  exact  geographic  location  of  each  partial 
route  required  an  assumption  of  spatial  pattern  metrics  that  coincided  with  5  equally 
divided  pieces,  which  may  not  have  represented  the  true  locations  used  by  the  volunteers. 
Prior  research  supports  the  concept  that  the  BBS  data  for  entire  routes  be  used.  For 
example,  population  trend  analyses  have  typically  been  based  on  yearly  counts  for  entire 
routes  (Robbins  et  al.  1986),  and  comparisons  between  observers  found  almost  identical 
species  total  for  the  entire  50  stops  but  not  for  the  individual  stops  (Robbins  et  al.  1986). 

A  second  limitation  of  this  study  is  the  USGS  LULC  classification  scheme.  The 
classification  scheme  of  the  USGS  LULC  data  was  not  developed  with  the  habitat 
requirements  of  specific  wildlife  species  in  mind  (Hepinstall  and  Sader  1997).  The 
classification  scheme  appears  adequate  for  the  forest  interior  and  forest  interior/edge 
species,  but  is  deficient  for  the  species  that  prefer  scrub/shrub  or  other  categories  not 
differentiated  in  the  USGS  LULC  data.  This  lack  of  specificity  may  explain  the  relatively 
weak  ability  of  the  composition  component  incorporating  the  classification  variables  to 
account  for  the  variance  in  bird  species  abundance  that  prefer  the  scrub/shrub  habitat. 
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Also,  local  scale  vegetation  characteristics  were  not  included  in  this  study  and  hence 
this  study  does  not  provide  managers  a  perspective  on  the  relative  importance  of  local 
versus  landscape  scale  habitat  for  the  bird  species.  Local  scale  habitat  typically  includes 
variables  to  describe  the  local  vegetation  structure,  such  as  number  of  living  stems  in 
different  size  classes,  percent  canopy  cover,  and  percentage  of  downed  logs.  Studies 
such  as  Knick  and  Rotenberry  (1995),  Bolger  et  al.  (1997),  Howell  et  al.  (2000),  and 
Pearson  (1993),  have  incorporated  both  landscape  and  local  habitat  variables.  As  might 
be  expected,  the  strength  of  the  results  for  some  of  the  species  in  the  studies  rose  with  the 
addition  of  the  local  habitat  variables. 

A  third  limitation  of  this  study  is  the  uncertainty  involved  in  selecting  the  best 
extent  for  use  in  spatial  pattern  metric  calculation.  Both  0.4  km  and  10  km  extents  were 
used  in  the  analysis,  with  the  0.4  extent  apparently  providing  the  best  results.  The 
inclusion  of  the  wide  extent  was  supported  by  prior  findings  indicating  that  a  10  km 
radius  best  explains  the  distribution  of  cowbirds  and  encompasses  the  average  distance 
moved  by  female  cowbirds  during  the  breeding  season  (Howell  et  al.  2000).  In  contrast 
to  studies  by  Robinson  et  al.  (1995),  and  Howell  et  al.  (2000),  the  current  study  did  not 
provide  strong  results  for  the  10  km  extent.  The  designation  of  the  10  km  radius  makes 
intuitive  sense  for  the  Robinson  et  al.  (1995)  study,  which  actually  looked  at  nest 
parasitism  by  the  Brown-headed  Cowbird.  However,  for  studies  such  as  the  current  study 
and  Howell  et  al.  (2000),  the  rationale  for  the  designation  of  a  10  km  radius  is  not  as 
meaningful,  since  the  specific  purpose  was  not  to  assess  nest  predation  or  brood 
parasitism  by  the  Brown-headed  Cowbird. 
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A  fourth  limitation  of  this  study  is  the  interpretation  of  the  results  of  the 
statistical  analysis.  Care  must  be  exercised  in  interpreting  results  of  the  stepwise 
regression  and  logistic  regression.  While  the  regression  models  cannot  be  interpreted  as 
cause  and  effect,  they  can  demonstrate  the  relative  importance  of  each  of  the  components 
included  in  the  models.  The  landscape  features  selected  by  a  species  may  be  masked  by 
averaging  spatial  pattern  measures.  All  BBS  routes  used  in  this  study  included  different 
habitat  types  with  different  patch  distributions.  The  spatial  pattern  metrics  are  means  per 
route  or  partial  route  and  provide  an  overall  picture  of  the  landscape  structure;  however, 
variability  within  a  route  may  be  masked  by  the  mean.  This  statistical  technique  does 
however  provide  an  important  first  means  of  assessing  variation  in  species  abundance  at 
various  spatial  extents  and  provides  an  overall  picture  of  landscape  habitat  characteristics. 
It  can  be  expected  that  extensive  detailed  studies  of  each  species  would  reveal  finer 
mechanisms  of  habitat  selection. 

Suggested  Refinements  in  the  Method 

The  BBS  data  were  utihzed  in  this  study  with  the  full  realization  that  there  are 
several  potential  sources  of  bias  in  the  data  set.  In  spite  of  these  limitations  the  unique 
geographic  and  temporal  extent  of  the  BBS  data  provides  a  strong  argument  for  its 
utilization.  It  is  recommend  that  the  species  with  models  with  R2  >  35%  should  be 
studied  further,  particularly  for  areas  that  have  temporal  BBS  and  landscape  data.  If 
applicable  to  the  species,  additional  landscape  variables,  such  as  elevation  and  isolation 
variables  such  as  nearest  neighbor  distance,  could  be  added  to  the  analysis.  The  models 
generated  could  then  be  applied  to  the  landscape  data  for  different  timeframes  to 
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determine  if  the  information  can  be  used  to  explain  population  declines  as  seen  in  the 
BBS  in  terms  of  habitat  suitability.  Nevertheless,  the  R2  values  will  still  probably  be 
moderate.  If  stronger  associations  are  to  be  expected,  the  bird  census  data  will  need  to  be 
designed  and  selected  for  the  particular  species  at  hand. 

The  analytical  technique  used  in  this  study  could  also  be  refined  by  improvements 
in  the  categorization  of  the  landcover  data  set.  Categories  that  reflect  the  specific 
landscape  requirements  of  the  birds  under  study  should  improve  the  results  of  the 
analysis  by  capturing  the  pertinent  variation.  This  is  particularly  needed  for  the  species 
that  prefer  scrub/shrub  or  other  categories  not  differentiated  on  the  USGS  LULC  data  set. 
The  current  category  in  the  USGS  LULC  data  set  most  closely  affiliated  with  these  birds 
is  the  category  of  rangeland  (land  where  the  potential  natural  vegetation  is  predominantly 
grasses,  grasslike  plants,  forbs,  or  shrubs).  However,  this  category  makes  up  less  than  1 
%  of  the  landcover/landuse  in  the  study  area.  Most  of  the  rangelands  in  the  US  are  in  the 
western  range  (Anderson  et  al.  1976).  Also,  because  the  same  12  independent  variables 
of  landscape  structure  were  used  for  all  15  bird  species,  various  important  variables  such 
as  elevation,  standing  water,  and  local  scale  vegetation  variables  that  strongly  influence 
the  distribution  of  some  birds,  were  not  included.  For  these  reasons,  the  R2  values  are 
relatively  low  for  all  species,  and  the  statistical  models  generated  have  limited  utility  in 
predicting  bird  abundance. 

The  method  of  analysis  can  be  further  refined  by  careful  consideration  of  the  scale 
and  extent  used  in  the  analysis.  If  a  study  considers  a  mixture  of  bird  species,  as  is 
addressed  in  this  research,  defining  the  area  or  extent  that  corresponds  to  a  single  scale 
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for  all  species  can  be  a  challenge.  There  is  apparently  no  single  best  scale  for 
investigating  avian  communities  and  landscape  structure.  Both  narrow  and  wide  spatial 
extents  were  used  in  this  study.  The  narrow  buffer  distance,  which  coincided  with  the 
BBS  design,  provided  the  best  statistical  results.  Ideally,  the  buffer  extent  used  should 
represent  the  landscape  over  which  the  birds  detected  are  breeding.  The  buffer  extents 
used  in  this  study  represented  a  compromise,  since  determining  the  extent  of  the  breeding 
area  is  dependent  on  a  complex  interaction  of  the  species  requirements  and  the  makeup  of 
the  landscape  itself.  The  selection  of  a  buffer  more  in  line  with  the  species  home  ranges, 
if  known,  is  more  meaningful.  For  example,  the  home  range  for  the  Wood  Thrush  is 
approximately  150m,  about  1.2  km  for  the  Prairie  Warbler,  about  1. 8-2.1  km  for  the 
Loggerhead  Shrike,  and  about  .8-1.6  km  for  the  Northern  Bobwhite  (Hamel  1992). 

These  figures  are  more  in  line  with  the  narrow  extent  of  .4  km  rather  than  10  km.  Since 
the  home  ranges  varied  for  the  species  in  the  current  study,  and  the  home  range  was 
unknown  for  over  50%  of  the  species  studied,  the  home  ranges  were  not  used  in  this 
study.  However,  in  future  studies,  it  may  be  better  to  concentrate  on  a  specific  species  or 
group  of  species  with  similar  habitat  needs  and  home  ranges.  Regardless,  the  10  km 
buffer  appeared  excessive  and  an  intermediate  buffer,  regardless  if  the  homerange  is 
known,  would  provide  a  better  reference  in  future  research. 

The  definition  of  the  extent  over  which  the  landscape  metrics  are  calculated 
determines  the  division  of  polygons  when  the  area  is  clipped  from  the  broader  database. 

In  this  study,  the  most  dramatic  effect  of  the  difference  in  the  narrow  and  wide  buffer 
extents  was  related  to  patch  size,  patch  density,  and  edge  density.  The  patches  in  the 
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narrow  buffer  areas  had  a  greater  chance  of  being  divided  when  the  area  was  clipped 
from  the  USGS  LULC  database.  This  is  illustrated  in  Figure  23,  which  shows  the  narrow 
and  wide  buffer  areas  around  a  partial  BBS  route.  It  demonstrates  how  one  large  forest 
patch  was  divided  into  smaller  patches  as  it  was  clipped  from  the  narrow  and  wide  buffer 
extents.  Turner  et  al.  (1989,  1990)  and  Wiens  (1989)  had  investigated  this  in  prior 
research,  concluding  that  quantitative  changes  in  measurement  across  spatial  scales  will 
differ  according  to  how  scale  is  defined.  Thus,  the  definition  of  scale,  in  terms  of  extent 
and  grain,  is  an  important  consideration  in  landscape  scale  studies.  Unfortunately,  there 
is  no  single  best  scale  for  investigation  of  all  avian  species  and  landscape  patterns.  The 
proper  scale  is  dependent  on  the  objective  of  the  analysis,  the  specific  species  and 
habitats  involved,  and  the  processes  that  are  believed  to  be  important.  One  errs  not  by 
advocating  that  a  particular  scale  may  be  useful  for  examining  bird  species  patterns,  but 
by  forcing  investigations  of  specific  patterns  into  a  scale  of  analysis  that  is  improper 
(Wiens  1989).  It  is  possible  that  the  methods  used  in  the  current  research  could  be 
improved  by  further  investigation  of  the  optimum  scale.  By  varying  the  scale  used  in  the 
study,  utilizing  a  finer  grained  dataset  and  an  intermediate  extent,  the  further  research 
would  likely  provide  a  more  unified  understanding  of  how  the  change  in  scale  affects  the 
associations.  This  in  turn,  might  shed  light  on  the  selection  of  the  scale  that  is 
appropriate  for  a  given  investigation  and  improve  the  method  of  analysis. 

Summary  and  Conclusions 

This  study  has  determined  that  landscape  spatial  patterns  from  broad  scale  data 
are  useful  indicators  of  habitat  suitability  for  the  species  under  study.  However,  the 
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multiple  regression  analysis  indicated  wide  variation  among  the  species  in  how  spatial 
pattern  metrics  affected  the  species  abundance.  For  five  of  the  species  examined,  the 
spatial  pattern  metric  models  generated  are  appropriate  for  the  assessment  of  the 
relationship  of  landscape  changes  and  population  declines. 

This  research  determined  that  a  comprehensive  set  of  landscape  metrics  provides 
a  way  to  produce  an  all-encompassing  quantitative  description  of  the  landscape.  The 
landscape  structure  was  quantified  using  a  suite  of  12  spatial  pattern  metrics  calculated 
from  USGS  LULC  data.  The  metrics  for  all  routes  were  pooled  and  summarized  into 
uncorrelated  landscape  structure  components  using  principal  components  analysis  (PC A) 
techniques.  Multiple  regression  techniques  showed  that  3  components,  landscape 
composition,  forest  configuration,  and  landscape  diversity/interspersion,  are  important  in 
the  abundance  of  some  of  the  bird  species. 

The  unexplained  variance  (>50%)  of  the  models  developed  in  this  study  probably 
reflects  the  inherent  variability  introduced  by  the  utilization  of  the  existing  sources  of 
bird  and  landscape  data.  The  existing  BBS  data,  which  may  be  biased  toward  roadside 
birds  that  sing  loudest  and  most  frequently,  may  also  be  confounded  by  unknown 
variability  among  volunteer  observers.  Also,  the  existing  classification  scheme  of  the 
USGS  LULC  data  was  not  developed  with  the  habitat  requirements  of  specific  wildlife 
species  in  mind  and,  as  such,  appears  adequate  for  some  of  the  forest  interior  and  forest 
interior/edge  species,  but  not  for  scrub/shrub  species.  Another  possible  source  of 
variability  is  the  suboptimization  of  the  extent  selected  for  use  in  the  spatial  pattern 


metric  calculation. 


It  can  be  concluded  from  this  study  that  spatial  pattern  metrics  from  broad  scale 
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data  are  appropriate  for  assessing  habitat  suitability;  however,  with  low  to  moderate  R 
values.  If  stronger  associations  are  to  be  expected,  bird  census  data  will  need  to  be 
designed  and  selected  for  the  particular  species  at  hand.  Also,  landcover  categories  that 
reflect  the  specific  landscape  requirements  of  the  birds  under  study  should  improve  the 
results  of  the  analysis  by  capturing  the  pertinent  variation.  This  will  require  a  finer  grain 
dataset  that  can  differentiate  the  specific  landcover  types.  The  analysis  might  also  be 
improved  if  an  intermediate  buffer  radius  was  used,  since  this  research  was  constrained 
by  the  physical  extent  of  the  landscape  study  area.  Ideally,  the  buffer  extent  used  should 
represent  the  landscape  over  which  the  birds  detected  are  breeding.  The  selection  of  a 
buffer  more  in  line  with  the  species  home  ranges,  if  known,  would  be  most  meaningful. 
In  future  studies,  it  may  be  better  to  concentrate  on  a  specific  species  or  group  of  species 
with  similar  habitat  needs  and  home  ranges. 

The  conservation  of  bird  populations  is  of  great  importance.  Conservation 
strategies  to  reverse  the  declines  of  forest  bird  populations  will  require  knowledge  of 
habitat  requirements  across  the  range  of  species.  It  has  been  concluded  from  this  study 
that  spatial  pattern  metrics  from  broad  scale  data  can  provide  a  good  first  approximation 
for  assessing  habitat  suitability.  This  provides  an  important  avenue  towards  proper 
assessment  of  habitat  suitability  and  towards  the  development  of  broad-scale 
conservation  management  practices. 


j 
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Figure  22.  Landscape  configuration  metrics  for  two  sample  landscapes  that  have  the  same  mean  patch  size. 


Figure  23.  One  5  mile  partial  breeding  bird  survey  route  with  a  0.4  and 
1 0  km  buffer  radius,  with  USGS  land  use/land  cover  data  shown  in  the 
background. 
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The  following  provides  a  description  of  the  preferred  habitat  of  the  15  species 
studied  in  this  research.  A  description  of  the  species  population  trends  is  also  discussed. 

Prairie  Warbler  (6730) 

This  species  is  often  listed  as  an  example  of  alarming  decline  among  Neotropical 
migrants.  The  Prairie  Warbler  breeds  in  shrubby  old  fields,  early-stage  regenerating 
forests,  and  other  early  successional  habitats.  Based  on  BBS  data  for  physiographic 
regions,  1 1  showed  decreases  and  2  showed  increases.  The  region  that  includes  this 
study  area  showed  declines  only  in  uplands.  Lowland  populations  were  stable,  and  even 
increased  in  certain  areas  (James  et  al.  1992).  Since  colonial  times,  deforestation  has 
created  extensive  breeding  habitat,  which  is  now  being  lost  to  urbanization  and 
reforestation.  Suppression  of  fires,  which  permits  forest  regeneration  may  play  a  factor 
in  habitat  loss  (Nolan  1999). 

Northern  Bobwhite  (2890) 

The  Northern  Bobwhite,  a  non-migratory  bird  that  is  considered  to  be  one  of 
North  America's  most  important  game  birds,  is  widely  distributed  throughout  the  eastern 
US  and  Mexico.  Its  preferred  habitat  includes  agricultural  fields,  grasslands,  open 
parklike  pine  and  mixed  pine-hardwood  forests.  According  to  BBS  data,  Christmas  Bird 
Count,  and  state  game  agencies,  there  is  evidence  of  widespread  decline  throughout  the 
US  (Brennan  1999).  Annual  declines  from  1966-1988  averaged  3%  in  eastern  US,  0.7% 
in  central  US,  and  1.8%  overall  in  US  (Droege  and  Sauer  1990).  The  declines  are 
primarily  attributed  to  habitat  loss  from  changing  land  uses  in  agriculture,  forestry,  and 
expanding  sububanization.  Specifically,  the  clean  farming  with  nearly  all  weeds 
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removed  and  high-density  pine  plantations  have  had  negative  habitat  results.  Also, 
changes  in  forestry  and  agricultural  land  use  have  resulted  in  fragmentation  of  the 
Northern  Bobwhite's  habitat  on  a  broad  landscape  scale  (Brennan,  1999).  Frequent 
vegetation  disturbance  from  prescribed  fire  and  mechanical  disturbances  is  needed  for 
maintaining  Bobwhite  populations  (Landers  and  Mueller  1986). 

Field  Sparrow  (5630) 

The  field  sparrow  breeds  in  scattered  samplings  or  shrubs  in  weedy  habitats, 
overgrown  fields,  wood  margins,  hedgerows,  and  thickets  of  the  eastern  US  and  southern 
Canada.  It  is  a  partial  migrant,  with  some  remaining  in  their  breeding  grounds  in  the 
winter  while  others  move  farther  south.  For  the  study  area,  it  is  considered  a  non- 
migratory  species.  According  to  the  BBS  data,  there  has  been  a  nationwide  decline  of 
3.4%  per  year  between  1966  and  1991.  Changes  in  land  use  are  the  primary  factors 
affecting  numbers.  The  presence  of  suitable  habitat  is  necessary  to  maintain  local 
populations  (Carey  et  al.  1994). 

Loggerhead  Shrike  (6220) 

Throughout  most  of  the  southern  part  of  its  range,  which  includes  the  study  area, 
the  loggerhead  shrike  is  resident.  It  inhabits  open  habitats  with  short  vegetation,  such  as 
grasslands,  pastures  with  fencerows,  mowed  roadsides,  agricultural  fields,  riparian  areas, 
and  open  woodlands.  According  to  the  BBS  data,  there  is  a  current  decrease  of  about 
3.5-5%/yr  across  the  range.  This  is  one  of  the  few  species  to  show  significant  declines  in 
most  states,  provinces,  and  physiographic  strata  of  the  continent  (Robbins  et  al.  1986). 
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Population  loss  can  be  partially  attributed  to  declines  in  hay  crops,  pasture  land,  and 
hedgerows  (Yosef  1996). 

Eastern  Kingbird  (4440) 

The  Eastern  Kingbird  is  the  most  widely  distributed  of  the  kingbirds  that  breed 
north  of  Mexico.  It  is  a  migratory,  aerial  hawking  insectivore  that  breeds  in  open 
environments,  usually  in  fields  with  scattered  shrubs  and  trees,  and  along  woodland  edges 
in  forested  regions.  According  to  the  BBS,  the  data  indicate  that  the  continental 
population  of  the  kingbird  didn't  change  significantly  between  1966-1991.  By  region, 
though,  several  significant  trends  appear:  no  significant  changes  in  abundance  in  eastern, 
central  or  w.  North  American  from  1966-1991,  but  from  1982-1991  abundance  in  e. 

North  America  declined  significantly,  but  also  increased  significantly  in  central  and 
western  regions.  Since  1982,  trends  have  shown  significant  recoveries  in  Alabama, 

South  Carolina,  and  Tennessee,  but  most  states  continue  to  show  declines.  Population 
declines  in  east  and  south  North  America  are  probably  caused  by  habitat  loss  resulting 
from  human  development  and  natural  plant  succession.  The  decline  of  small  farms  since 
the  1960s-1990s  has  resulted  in  the  loss  of  open  space.  Forest  succession  has  also 
resulted  in  the  loss  of  significant  habitat  (Murphy  1996). 

White-eyed  Vireo  (6310) 

The  preferred  habitat  of  the  White-Eyed  Vireo  is  dense  secondary  deciduous 
scrub,  wood  margins,  and  overgrown  pastures.  It  is  a  migratory  songbird  that  is  primarily 
detected  by  ear  than  by  eye.  Based  on  BBS  data,  significant  declines  in  numbers  over  its 
range  was  -0.72%  per  year  for  the  long-term  period  (1966-1988)  and  -2.99%  and  for  the 
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short-term  period  (1978-1988).  In  the  eastern  US  they  found  significant  short-term 
but  not  long-term  declines.  Nesting  areas  are  frequently  cleared  by  humans.  With 
significant  declines  in  the  White-eyed  Vireo  population  in  some  regions,  it  is  suggested 
that  some  measures,  such  as  promoting  suitable  scrub  habitat  either  by  leaving  open  areas 
to  grow  or  by  opening  some  forested  areas  through  partial  cutting,  be  addressed.  The 
measures  have  not  been  tested  (Hopp  et  al.  1995). 

Eastern  Wood-Pewee  (4610) 

The  Eastern  Wood-Pewee  is  a  migratory  species  that  breeds  in  wooded  habitats  of 
the  eastern  US  and  Canada.  It  breeds  in  virtually  every  type  of  wooded  habitat  in  the 
east,  from  urban  shade  trees,  roadsides,  woodlots,  and  orchards  to  mature  forest  (McCarty 
1996).  BBS  data  show  a  significant  decrease  in  population  from  1965-1993  (Robbins  et 
al.  1986).  Declines  of  35.6%  were  found  during  1966-1993  while  there  was  a  13.4% 
decline  from  1984-1993  (Price  et  al.  1995).  The  species  uses  both  edge  and  forest 
interior  for  breeding.  It  is  apparently  not  sensitive  to  forest  fragmentation  when  choosing 
breeding  sites  (Blake  and  Karr  1987).  Heavy  browsing  of  forests  by  white-tailed  deer 
may  play  a  part  in  its  decline. 

Yellow-throated  Vireo  (6280) 

The  Yellow-Throated  Vireo  is  a  migratory  species  that  breeds  in  the  eastern  US  in 
edge  habitats  of  both  bottomland  and  upland  deciduous  and  mixed  deciduous-coniferous 
forests.  Habitats  include  forest  edges  of  streams,  rivers,  swamps,  treefall  gaps,  and  roads, 
and  woodland  habitats  of  parks  and  towns  (Rodewald  and  James  1996).  The  species  has 
disappeared  or  decreased  at  several  small  forest  reserves  in  the  eastern  US  (Askins  et  al. 
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1990).  Overall  population  trends  from  the  BBS  data  showed  significant  increase  of 
1.1  %/year  during  1966-1994.  During  this  timeframe,  species  increased  1.3%/year  in  the 
eastern  part  of  the  range,  but  declined  0.9%/year  in  the  central  region.  While  this  species 
is  often  associated  with  forest  edge  habitat,  large  blocks  of  forest  are  necessary  for 
successful  breeding.  Extensive  clear  cutting  will  adversely  affect  the  species  (Rodewald 
and  James  1996). 

Hooded  Warbler  (6840) 

The  Hooded  Warbler  is  a  small  migratory  songbird  that  breeds  in  southernmost 
Canada  and  the  eastern  US  and  winters  in  Central  America.  It  prefers  mixed  hardwood 
forest  in  the  north  and  cypress-gum  swamps  in  the  south.  It  is  considered  a  forest-interior 
species  because  it  is  restricted  to  larger  woodlots.  It  is  declining  in  only  a  few  parts  of  its 
breeding  range  (Evans  Ogden  and  Stutchbury  1994).  According  to  BBS  data,  regionally 
there  is  much  variation  in  the  population  trends.  From  1966-1991,  Georgia  was  the  only 
state  to  show  a  significant  decrease.  For  the  Easter  Region  of  BBS,  there  was  a 
significant  long-term  increase  (1.58%/yr)  (Sauer  and  Droege  1989).  In  Canada  the 
species  is  considered  threatened  because  suitable  habitat  is  becoming  increasingly  scarce 
and  fragmented.  Since  the  Hooded  Warbler  is  an  area-sensitive  forest  songbird,  it  is 
threatened  on  breeding  grounds  primarily  by  forest  fragmentation  (Evans  Ogden  and 
Stutchbury  1994). 

Wood  Thrush  (7550) 

The  Wood  Thrush  is  a  Neotropical  migratory  bird.  Its  preferred  habitats  includes 
interior  and  edges  of  deciduous  and  mixed  forests,  especially  well-developed,  uplands.  It 
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is  more  likely  to  occur  in  larger-area  forests.  BBS  continent-wide  data  show 
significant  decline  in  abundance  of  1.7%  per  year  from  1966-1994  (Sauer  and  Droege 
1992).  Fragmentation  of  forests  is  possible  cause  for  decline  (Roth  et  al.  1996). 

Yellow-billed  Cuckoo  (3870) 

The  Yellow-billed  Cuckoo  is  a  migratory  bird  that  breeds  in  deciduous  forests, 
bottomland  woods,  woodland  thickets,  and  other  hardwood  forests.  It  generally  avoids 
coniferous  woods.  It  prefers  extensive  forests  (Hamel  1992).  From  1980  to  1994, 
populations  declined  in  most  states.  Eastern  US  populations  suffered  3.2%  annual 
decline.  The  species  is  sensitive  to  habitat  fragmentation  and  degradation  of  riparian 
woodlands  due  to  agricultural  and  residential  development  (Dobkin  1994). 

Kentucky  Warbler  (6770) 

The  Kentucky  Warbler  is  a  migratory  species  that  breeds  primarily  in  the 
southeastern  US.  It  prefers  rich,  moist,  deciduous  forests,  and  bottomland  hardwoods. 
Studies  of  forest  fragmentation  in  Missouri  indicate  that  blocks  of  suitable  habitat  (at 
least  500  ha)  are  necessary  for  successful  breeding  (Gibbs  and  Faaborg  1990).  It  is  rarely 
observed  in  agricultural  habitats.  BBS  data  suggest  that  since  1980  the  Kentucky 
Warbler’s  continentwide  population  has  been  slowly  decreasing,  however,  local  increases 
and  range  expansion  also  seem  to  be  occurring  (McDonald  1998).  Forest  management 
practices  that  promote  a  dense  understory  and  ground  cover  should  enhance  habitat  for 
this  species  (Bushman  and  Therres  1988). 
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Brown-headed  Nuthatch  (7290) 

The  Brown-headed  Nuthatch  is  a  non-migratory  bird  that  lives  in  pine  forest  of 
the  southeastern  US  (Withgott  and  Smith  1998).  According  to  BBS  data,  between  1966 
and  1996,  there  has  been  a  significant  -2.2%/year  decline  throughout  the  range  (Sauer  et 
al.  1997).  The  decline  is  primarily  attributed  to  the  altered  pine  forests  by  commercial 
logging.  Old-growth  pine  forest  is  almost  gone,  with  even-aged  stands  of  younger  pines 
replacing  them.  Fire  suppression  has  also  negatively  affected  habitat  suitability  allowing 
hardwood  to  flourish  (Engstrom  et  al.  1984). 

Prothonotary  Warbler  (6370) 

The  Prothonotary  Warbler  is  a  migratory  species  that  breeds  primarily  in  the 
southeastern  US.  It  prefers  bottomland  hardwood  forests  and  is  considered  an  interior 
species.  BBS  data  indicated  significant  overall  decreases  of -1.6%  annually,  between 
1966  and  1996.  Declines  are  attributed  primarily  to  loss  of  habitat.  Bottomland 
hardwood  forests  have  been  logged  or  converted  to  pasture  cropland  throughout  the 
southeastern  US  (Petit  1999). 

Brown-headed  Cowbird  (4950) 

The  Brown-headed  Cowbird  is  North  America’s  most  well  known  brood  parasite. 
For  the  study  area,  it  is  considered  a  non-migratory  species.  Its  preferred  habitat  includes 
open  woods,  margins,  thickets,  agricultural  and  residential  areas.  BBS  data  show  for  the 
timeframe  between  1965-1987  significant  increases  for  Georgia,  North  Carolina,  Iowa, 
North  Dakota,  Utah,  Colorado  and  significant  decreases  for  Minnesota,  Michigan, 
Wisconsin,  New  York,  Rhode  Island,  Ohio,  Ontario,  West  Virginia,  Tennessee,  New 
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Brunswick,  Oklahoma,  Texas  and  Oregon.  Because  of  forest  fragmentation  increases 
edge  and  provides  access,  forest  interior  species  have  become  exposed  to  the  Brown¬ 
headed  Cowbird  (Lowther  1993). 


Bird  Abundance 

Table  A1  lists  the  mean  yearly  abundance  for  each  species  for  each  full  BBS 
route,  while  Table  A2  lists  the  mean  yearly  abundance  for  each  species  for  each  partial 
BBS  route. 


Table  Al.  Mean  yearly  abundance  for  each  species  for  each  of  the  53  full 
BBS  routes. 
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Eg 

AL 

10 
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KB 
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Eg 

AL 

11 
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0.0 
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■& 

AL 
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0.2 
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AL 
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AL 

16 

HI 
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1.0 
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0.4 

H 

AL 

17 
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Eg 

AL 
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AL 
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AL 
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AL 
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AL 
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0.0 

12.0 

AL 
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IH 
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0.0 

0.0 
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Table  A1 .  (cont’d)  Mean  yearly  abundance  for  each  species  for  each  of 
the  53  full  BBS  routes. 
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Table  A2.  Mean  yearly  abundance  for  each  species  for  each  of  the  265  partial 
BBS  routes. 


Table  A2.  (cont’d)  Mean  yearly  abundance  for  each  species  for  each  of 
the  265  partial  BBS  routes. 
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Table  A2.  (cont’d)  Mean  yearly  abundance  for  each  species  for  each  of 
the  265  partial  BBS  routes. 
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Table  A2.  (cont’d)  Mean  yearly  abundance  for  each  species  for  each  of 
the  265  partial  BBS  routes. 
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Table  A2.  (cont’d)  Mean  yearly  abundance  for  each  species  for  each  of 
the  265  partial  BBS  routes. 


APPENDIX  B 

CLIMATE  DATA  ANALYSIS 
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Methods 
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The  weather  conditions  during  the  years  of  this  study  (1970  -  1976)  were 
compared  to  a  long-term  period  of  weather  record.  This  was  done  to  determine  if  the 
years  included  in  the  study  were  atypical  when  compared  to  the  long-term  period  of 
record.  Pertinent  weather  stations  were  selected  from  Earthlnfo  Environmental  Database 
Summary  of  the  Day  CDs.  This  was  done  separately  for  Alabama,  Georgia,  and  Florida 
and  then  combined  into  a  comprehensive  database.  The  climate  data  for  the  study 
timeframe  were  based  on  weather  stations  that  had  data  for  the  years  1970  to  1976  while 
the  long-term  climate  data  were  based  on  weather  stations  that  had  a  minimum  of  thirty 
years  of  record.  For  each  weather  station  and  timeframe,  the  following  data  were 
obtained  directly  or  were  calculated  from  data  on  the  CD:  latitude,  longitude,  minimum 
yearly  mean  temperature  in  degrees  Fahrenheit,  and  maximum  yearly  mean  temperature 
in  degrees  Fahrenheit.  The  number  of  weather  stations  used  for  the  following  four 
scenarios  was: 

174  minimum  yearly  mean  temperature  for  long-term  period  of  record 
123  minimum  yearly  mean  temperature  for  the  study  timeframe 
174  maximum  yearly  mean  temperature  for  long-term  period  of  record 
123  maximum  yearly  mean  temperature  for  the  study  timeframe 

A  GIS  was  used  to  provide  a  geographic  view  of  the  climate  conditions  in  the 
study  area  for  the  study  timeframe  and  a  longer  period  of  climate  record.  This  provided  a 
general  geographic  overview  of  the  climate  conditions  of  the  study  timeframe  and  the 
longer  period  of  climate  record. 
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Separate  GIS  coverages  were  developed  for  the  minimum  yearly  mean 
temperature  for  long-term  period  of  record,  minimum  yearly  mean  temperature  for  the 
study  timeframe,  maximum  yearly  mean  temperature  for  long-term  period  of  record, 
maximum  yearly  mean  temperature  for  the  study  timeframe.  All  weather  stations  within 
each  GIS  coverage  were  used  to  develop  a  surface  model.  First,  a  triangulated  irregular 
network  (TIN)  was  developed.  A  TIN  model  is  a  set  of  adjacent,  non-overlapping 
triangles  computed  from  irregularly  spaced  points.  Each  TIN  model  was  then 
transformed  into  a  lattice  file  with  a  grid  cell  size  of  250  meters.  The  grid  files  were  used 
to  illustrate  the  gradation  of  values  in  temperature.  Also,  to  show  the  differences  between 
the  two  timeframes,  the  values  for  the  study  timeframe  were  subtracted  from  those  of  the 
long-term  period  of  record  for  minimum  and  maximum  yearly  mean  temperature. 

Unlike  temperature,  precipitation  is  discontinuous  in  both  the  geographic  and 
temporal  variables.  As  a  result,  precipitation  data,  gathered  at  a  discreet  location  and 
over  discreet  intervals,  inherently  have  much  higher  degrees  of  variability.  A  general 
snapshot  of  the  precipitation  over  the  study  area  was  developed  by  creating  a  30-year 
time  sequence  of  total  annual  precipitation  for  3  locations  in  the  study  area.  Three 
weather  stations,  one  in  the  northern,  one  in  the  central,  and  one  in  the  coastal  portion  of 
the  study  area,  were  used  to  determine  the  total  annual  precipitation  in  1/100  of  an  inch 
increments  for  every  year  of  record. 
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Results  and  Conclusions 

The  minimum  and  maximum  yearly  mean  temperature  provide  two  measurements 
that  give  a  general  indication  of  the  weather  conditions  of  the  study  area.  The 
geographical  numerical  differences  were  primarily  between  positive  or  negative  1  to  2 
degrees  Fahrenheit  (Figures  B3  and  B6).  Also,  the  precipitation  time  sequence  that  was 
developed  for  three  weather  stations  in  the  study  area  did  not  have  visible  outliers  and 
appeared  to  be  compatible  with  the  trends  of  the  annual  precipitation  of  the  long-term 
period  of  record.  Therefore,  it  can  be  inferred  that  the  results  of  the  overall  dissertation 
study  should  not  be  affected  by  unusual  climate  conditions  in  reference  to  the  minimum 
and  maximum  yearly  mean  temperature  and  annual  precipitation.  More  specific 
discussion  is  listed  below. 

Figures  B1  and  B2  show  the  minimum  yearly  mean  temperature  for  the  study 
timeframe  and  the  long-term  period  of  record,  respectively.  Figure  B3  shows  the  map  of 
the  numerical  difference  between  the  mean  temperature  values  for  the  two  timeframes. 

As  shown  on  Figures  B1  and  B2,  the  mean  temperature  values  for  the  two  timeframes 
range  from  approximately  44  to  62  degrees  Fahrenheit.  As  would  be  expected,  it  is 
apparent  that  the  lowest  temperatures  are  in  the  northern  portion  of  the  study  area  with  a 
gradation  to  higher  temperatures  in  the  southern  portion  of  the  study  area.  Figure  B3 
shows  the  differences  in  the  two  timeframes  ranges  from  -3  to  6  degrees  Fahrenheit,  but 
with  the  majority  of  the  study  site  varying  only  about  positive  or  negative  1  to  2  degrees 
Fahrenheit.  When  an  average  of  all  grid  cells  depicting  the  minimum  yearly  mean 
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temperature  was  calculated  for  the  study  area,  the  value  was  the  same  for  the  study 
timeframe  and  the  long-term  period  of  record. 

Figures  B4  and  B5  show  the  maximum  yearly  mean  temperature  for  the  study 
timeframe  and  the  long-term  period  of  record,  respectively.  Figure  B6  shows  the  map  of 
the  numerical  difference  between  the  maximum  yearly  mean  temperature  values  for  the 
two  timeframes.  Figures  B4  and  B5  show  the  maximum  yearly  mean  temperature  values 
for  the  two  timeframes  range  from  about  68  to  82  degrees  Fahrenheit.  The  temperatures 
are  lowest  in  the  north  and  gradually  get  higher  in  the  southernmost  portion  of  the  study 
area.  Figure  B6  shows  the  differences  vary  from  -2  to  5  degrees  Fahrenheit.  The 
majority  of  the  study  site  varies  only  about  positive  1  to  2  degrees  or  negative  1  degree 
Fahrenheit  between  the  two  timeframes.  When  an  average  of  all  grid  cells  depicting  the 
maximum  yearly  mean  temperature  was  calculated  for  the  study  area  for  the  two 
timeframes,  the  value  was  the  same  for  the  study  timeframe  and  the  long-term  period  of 
record. 

A  plot  (Figure  B7)  was  developed  to  show  the  annual  precipitation  for  each  year 
of  record  for  the  3  stations.  The  study  timeframe  is  marked  on  Figure  B7,  illustrating  the 
pattern  of  the  study  timeframe  in  comparison  to  the  long-term  period  of  record.  It 
appears  from  Figure  8  that  the  annual  precipitation  during  the  study  timeframe  follows 
the  same  trends  as  the  long-term  period  of  record  and  that  there  are  no  exceptional  values 
that  make  the  study  timeframe  unique. 
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Figure  B1 .  Minimum  Yearly  Mean  Temperature  in  Degrees  Fahrenheit  for 
the  Study  Time  Frame. 


Figure  B2.  Minimum  Yearly  Mean  Temperature  in  Degrees  Fahrenheit  for  the 
Long-term  Period  of  Record. 
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Figure  B3.  Difference  in  Minimum  Yearly  Mean  Temperature  in  Degrees 
Fahrenheit  for  the  Study  Timeframe  and  the  Long-term  Period  of  Record. 


Figure  B5.  Maximum  Yearly  Mean  Temperature  in  Degrees  Fahrenheit  for  the 
Long-term  Period  of  Record. 


Figure  B6.  Difference  in  Maximum  Yearly  Mean  Temperature  in  Degrees 
Fahrenheit  for  the  Study  Timeframe  and  the  Long-term  Period  of  Record. 


138 


1,920  1,928  1,936  1,944  1,952  1,960  1,968  1,976  1,984  1,992  2,000 

Year 

Figure  B7.  Plot  of  annual  precipitation  for  each  year  of  record  for  3  weather  stations. 


APPENDIX  C 

LANDSCAPE  SPATIAL  PATTERN  METRICS 
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Tables  C1-C4  list  the  value  of  each  of  the  12  metrics  for  each  study  site  for 
Extent  1  -4,  respectively.  There  are  some  study  sites  that  do  not  have  a  value  for  IJI.  This 
is  due  to  the  fact  that  FRAGSTATS  cannot  calculate  UI  for  samples  with  <  3  distinct 
patch  types. 


Table  Cl.  Value  of  each  of  the  12  metrics  for  each  study  site  for  Extent  1 
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36.28 

AL 

IB 

8.82 

0 

0 

0 

81.83 

369.37 

0.22 

11.26 

2.23 

24.37 

0.32 

54.09 

AL 

52.32 

0 

0.91 

0 

45.65 

73 

0.64 

7.49 

1.73 

26.38 

0.52 

34.52 

AL 

68.51 

0 

3.57 

0 

24.87 

39.71 

0.72 

12.48 

1.58 

8.54 

0.47 

36.29 

AL 

20.73 

0 

0 

0 

77.25 

300.14 

0.26 

10.22 

1.98 

25.23 

0.36 

26.46 

AL 

m 

36.28 

0 

0 

12.91 

49.5 

136.06 

0.46 

11.75 

1.85 

21.66 

0.61 

38.57 

AL 

m 

47.26 

0 

0 

2.74 

47.35 

57.24 

0.87" 

16.68 

1.71 

14.88 

0.55 

33.36 

AL 

ga 

38.91 

0.65 

0 

4.44 

52.53 

114.91 

0.5 

17.16 

1.92 

11.31 

0.57 

35.51 

AL 

m 

24.83 

0 

2.94 

18.26 

47.91 

200.65 

0.34 

11.63 

2.26 

21.67 

0.67 

51.65 

AL 

29 

39.62 

0 

0 

0 

56.37 

159.63 

0.35 

10.19 

1.73 

21 

0.52 

33.82 

AL 

78.43 

0 

0 

0 

19.29 

13.39 

1.44 

14.44 

1.49 

1.03 

0.35 

28.74 

AL 

IB 

73.9 

0 

0 

3.99 

19.73 

19.83 

1.2 

15.82 

1.4 

1.22 

0.41 

41.71 

AL 

gg 

67.92 

0 

0 

0.95 

26.01 

18.57 

1.45 

16.81 

1.51 

3.29 

0.47 

34.75 

AL 

IB 

29.33 

0 

0 

9.24 

61.2 

163.98 

0.43" 

9.98 

1.82 

25.51 

0.53 

45.23 

AL 

g 

15.34 

0 

6.23 

55.08 

26.65 

531.73 

0.15 

10.35" 

2.99 

24.03 

0.6 

52.77 

AL 

20.14 

0 

0 

0.24 

73.05 

336.09 

0.22 

5.77 

1.96 

29.69 

0.42 

55.89 

AL 

m 

25.24 

0 

0 

265.87 

0.28 

10.49 

2.08 

23.59 

0.41 

30.79 

AL 

g 

31.81 

0 

0 

1.21 

65.55 

211.26 

0.32 

9.01 

1.8 

27.23 

0.47 

32.59 

46.32 

0 

0.18 

0.93 

47.46 

57.71 

0.84 

17.12 

1.66 

10.45 

0.56 

37.17 

AL- 

IEB 

0 

9.15 

0 

2.74 

24.2 

129.91 

0.28 

9.65 

2.32 

16.72 

0.86 

75.68 

AL 

13.74 

0 

0 

0 

81.6 

289.14" 

0.28 

8.78 

2.04 

28.37 

0.31 

41.19 

AL 

IS 

0.96 

0 

0 

0 

2.42 

12.72" 

0.19 

2.11 

1.49 

0 

58.69 

AL 

■ 

62.28 

0 

0 

0.32 

32.71 

30.41 

1.09 

14.77 

1.56 

9.17 

0.5 

AL 

\m 

37.51 

2.33 

4.07 

0.56 

19.7 

47.92 

0.56 

11.08 

1.55 

5.06 

0.78 

59.58 

AL- 

~ 

63.83 

0 

0.07 

0 

31.56 

0.47 

11.22 

1.7 

1.97 

0.49 

35.7 

aT“ 

~ 

41.65 

0 

9.55 

4.66 

27.25 

|g 

0.62 

11.21 

1.53 

15.29 

0.73 

60.65 

AL- 

9 

80.59 

0 

9.17 

0.31 

8.11 

47.17 

0.37 

1.48 

6.78 

0.34 

61.17 

Table  C2.  Value  of  each  of  the  12  metrics  for  each  study  site  for  Extent  2 


State 

Rt 

%  Ag 

% 

Fw 

% 

Dec 

%  Ev 

% 

Mx 

MPS 

PD 

IJI 

# 

(ha) 

(no/ 

(m/ha) 

(%) 

(%) 

100  ha) 

AL 

1 

38.2 

0.57 

9.8 

1.13 

0.11 

11.3 

1.62 

50.6 

0.71 

43.6 

AL 

10 

12 

45.3 

BUS 

liAlill 

5.92 

1.88 

79.9 

0.69 

39.5 

EM 

11 

20.8 

0 

58.6 

1696 

0.04 

8.45 

1.67 

70.6 

0.6 

49.5 

AL 

12 

12.6 

■3! 

9.95 

55.3 

3233 

0.03 

5.81 

1.88 

81.1 

0.64 

50.7 

AL 

13 

15.6 

0 

0.74 

11.6 

62.3 

3792 

0.02 

9.58 

2.42 

67.1 

0.57 

47.1 

AL 

15 

18.1 

0 

0 

0.74 

74.2 

5642 

0.01 

12.6 

2.42 

57.1 

0.42 

33.7 

AL 

16 

18.1 

0 

0 

2.61 

77.6 

11113 

0.01 

2.65 

68.1 

0.36 

21.5 

AL 

17 

3.81 

0 

0 

2.78 

78.6 

8136 

0.01 

9.18 

2.34 

69 

0.37 

51.8 

AL 

2 

63.9 

0.07 

28.5 

205.7 

0.15 

9.95 

1.43 

35.6 

0.51 

27.4 

AL 

22 

43.5 

1.98 

3.1 

5.31 

41.3 

1057 

0.05 

11.8 

1.69 

47.2 

0.64 

AL 

23 

m 

0.02 

1.59 

0.42 

87.3 

17433 

0.01 

2.24 

0.23 

25.7 

AL 

24 

41 

0.19 

1.08 

6.95 

50 

809.7 

0.07 

11.9 

1.67 

51.9 

0.58 

24.5 

AL 

25 

21.3 

0.11 

1.16 

2.25 

70.6 

1898 

0.04 

12.6 

1.96 

58.3 

0.45 

30.3 

AL 

26 

37.9 

1.49 

BE 

1.3 

50.3 

606.2 

0.09 

14.4 

1.87 

43 

0.6 

41.4 

AL 

28 

20.9 

0.18 

53.6 

1614 

0.04 

11.7 

1.76 

58.9 

0.65 

48.8 

AL 

29 

24.6 

0 

1.12 

3.07 

69.3 

3825 

0.02 

8 

1.8 

69.8 

0.46 

31.7 

AL 

30 

50.3 

0.52 

0.15 

1.28 

45.3 

620.3 

0.08 

16.4 

1.95 

29 

0.54 

19.1 

AL 

1.87 

2.09 

6.39 

30.8 

81.45 

0.17 

mu 

1.51 

13 

0.59 

32.1 

AL 

32 

34.1 

2.87 

0.01 

3.57 

57.6 

2063 

0.03 

16.5 

2.05 

42 

0.55 

23.9 

AL 

34 

23.8 

1.12 

0.33 

15.8 

57.6 

3026 

0.02 

8.79 

1.72 

MM 

0.59 

28.3 

AL 

2.85 

0.06 

50.8 

31.3 

10576 

0.01 

7.69 

2.25 

74.2 

0.63 

41.6 

AL 

36 

18.2 

0.21 

mg 

71.6 

4056 

0.02 

P iM 

1.89 

72.7 

0.45 

36.6 

AL 

37 

8.79 

21.6 

0 

57.2 

9166 

0.01 

6.39 

0.61 

47 

AL 

0.03 

0 

19.1 

55.4 

2919 

0.03 

8.94 

1.79 

69.3 

0.61 

37.5 

AL 

31.7 

0.05 

1.37 

2.33 

57.8 

998.4 

13.2 

1.64 

48.1 

0.56 

31.6 

AL 

41 

5.36 

1.46 

lljXjgj 

0.03 

4.14 

1.81 

38.8 

0.48 

69.1 

AL 

iretiEl 

0 

0.35 

1.23 

86.1 

20702 

0 

6.99 

3.26 

76.6 

0.25 

45.4 

AL 

44 

4.45 

0.13 

0.51 

0.89 

51.6 

1884 

0.03 

9.69 

2.07 

55.2 

0.64 

56.3 

AL 

44.9 

0 

0.04 

0.36 

53.5 

495.9 

0.11 

16 

1.65 

36.6 

0.51 

13.4 

AL 

40.4 

1.4 

3.07 

1.27 

35.5 

422.8 

0.1 

11.6 

1.55 

36.5 

0.7 

50.1 

AL 

7 

30.4 

0.04 

39.5 

0.22 

29.2 

1676 

IKES 

74.9 

0.67 

29.9 

AL 

8 

30 

1.54 

11.2 

41.3 

1560 

0.04 

8.45 

IHB 

63.5 

0.72 

45.6 

AL 

9 

60.9 

0 

12.6 

1.35 

17.7 

276.7 

0.11 

9.45 

1.56 

38.5 

0.58 

41.5 

Table  C2.  (cont’d)  Value  of  each  of  the  12  metrics  for  each  study  site  for  Extent  2 


State 

Rt 

%  Ag 

% 

Fw 

% 

Dec 

vO 

m 

< 

% 

Mx 

MPS 

PD 

ED 

MSI 

TCAI 

IJI 

# 

(ha) 

(no/ 

(m/ha 

(%) 

<%) 

100  ha) 

FL 

1 

20.7 

12.5 

0.1 

54.7 

7.62 

2184 

0.03 

9.66 

1.66 

69.5 

0.64 

39.3 

FL 

2 

0.08 

5.84 

0.05 

14.8 

1.01 

IQQ 

66 

FL 

3 

12.1 

5006 

0.02 

(jQg 

0.58 

52.5 

FL 

4 

17.2 

21 

0 

55.8 

2.04 

2920 

0.03 

9.35 

1.9 

70.9 

0.61 

39.4 

FL 

5 

2.78 

14.5 

30511 

0 

3.3 

2.43 

88 

0.43 

33.4 

FL 

6 

gag 

2.54 

0 

52.4 

4.69 

997 

0.06 

16.7 

1.91 

42.4 

0.65 

43.1 

FL 

7 

8.59 

22.5 

0.04 

41.5 

25.3 

5666 

0.02 

6.82 

2.07 

79.3 

0.71 

45.2 

FL 

8 

4.44 

20 

0 

57.1 

15.2 

11156 

0.01 

5.12 

2.17 

84 

0.61 

40.1 

FL 

9 

0.02 

7.61 

0 

24 

10.2 

3470 

0.01 

4.35 

2.28 

76.4 

0.69 

56.7 

GA 

10 

55.2 

8.15 

13.5 

11.2 

9.47 

345.4 

0.12 

14 

1.63 

31.9 

0.65 

44.2 

GA 

16 

29 

0.14 

9.9 

10.9 

45.5 

2180 

0.03 

10.8 

1.7 

59.7 

0.69 

45.6 

GA 

22 

37.2 

1.74 

12.1 

6.19 

36.4 

856 

0.07 

16.6 

1.78 

36 

0.71 

50.3 

GA 

25 

63.6 

4.26 

5.5 

17.8 

6.41 

174 

0.19 

14.6 

1.64 

21.5 

0.55 

40.9 

GA 

26 

26.2 

0 

22.7 

14.3 

32.6 

2307 

0.03 

7.69 

1.75 

70.5 

0.75 

51.4 

GA 

28 

60.1 

1.94 

0.08 

19.6 

14.2 

0.1 

11.3 

1.69 

34.8 

0.58 

38 

GA 

33 

48.9 

4.58 

10.4 

7.22 

24.2 

104 

0.09 

4.38 

1.35 

21 

0.68 

47.9 

GA 

36 

25.7 

0.06 

11.9 

8.47 

51.8 

1390 

0.05 

13.4 

1.86 

57 

0.64 

40.4 

GA 

37 

21.9 

0 

16.3 

10 

49.1 

4618 

0.02 

7.16 

1.79 

73.1 

0.67 

42.6 

GA 

38 

20.7 

4.36 

1.25 

28.3 

43.8 

5593 

0.01 

8.59 

1.83 

68.7 

0.68 

43.4 

GA 

6 

35.5 

0.07 

4.14 

1128 

0.05 

11.9 

1.58 

52.8 

0.71 

44.4 

Table  C3.  Value  of  each  of  the  12  metrics  for  each  study  site  for  Extent  3. 
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State  1 

m 

Partial 

%  Ag 

%  Fw 

%  Dec 

%Ev 

%  Mx 

MPS 

PD 

ED 

MSI 

S1DI 

IJI 

mam 

m 

Rt 

(ha) 

(no/ 

(m/ha) 

(%) 

(%) 

100  ha) 

AL 

1st 

40.04 

0 

11.61 

0 

Egg! 

75.26 

0.73 

20.09 

1.8 

5.62 

0.64 

53.29 

AL 

KH 

2nd 

4.33 

0 

23.24 

17.86 

Efl&H 

0.15 

4.23 

2.21 

33.11 

0.61 

64.62 

AL 

|Q 

3rd 

9.86 

0 

33.16 

IBQI 

6.75 

2.31 

27.33 

0.56 

82.16 

AL 

KQ 

4th 

44.02 

0 

25.95 

0 

30.03 

16.92 

1.54 

8.28 

0.65 

91.14 

AL 

5th 

4.81 

0 

15.05 

12.13 

62 

608.78 

0.15 

8.64 

2.66 

26.1 

0.57 

60.19 

AL 

1st 

29.33 

0 

0 

0 

69.99 

93.51 

0.75 

13.67 

1.59 

23.31 

0.42 

30.23 

AL 

m 

2nd 

20.59 

0 

8.92 

25.46 

31.83 

223.6 

0.3 

15.07 

14.65 

0.77 

mtffl 

AL 

hi 

3rd 

16.25 

0 

4.86 

35.41 

43.48 

282.92 

0.3 

10.3 

1.65 

22.27 

0.66 

73.43 

AL 

4th 

0 

0 

0 

47.93 

80.85 

0.59 

12.87 

1.43 

11.43 

0.5 

AL 

5th 

16.34 

0 

3.36 

18.56 

61.74 

564.49 

0.15 

4.48 

2 

31.68 

0.56 

57.51 

AL 

m 

1st 

29.12 

0 

o" 

0 

54.69 

111.58 

0.49 

21.04 

2.43 

7.17 

0.6 

51.26 

AL 

2nd 

0 

0 

0 

0 

99.93 

614.25 

0.16 

0.29 

1.86 

36.67 

0 

AL 

m 

3rd 

0 

0 

22.93 

0 

77.07 

606.17 

0.16 

0 

1.83 

37.6 

0.35 

AL 

m 

4th 

0 

0 

46.37 

2.58 

51.05 

512.43 

0.2 

0 

1.48 

46.09 

0.52 

82.19 

AL 

m 

5th 

0 

0 

4.49 

40.31 

55.21 

618.97 

0.16 

0 

1.88 

36.2 

0.53 

40.64 

AL 

|Q 

1st 

7.89 

0 

0 

50.13 

273.19 

0.3 

10.12 

1.82 

24.34 

0.63 

72.36 

AL 

|Q 

2nd 

31.07 

0 

0 

47 

0 

61.01 

0.77 

20.62 

1.75 

3.22 

0.66 

60.1 

AL 

|Q 

3rd 

46.46 

0 

0 

12.47 

40.83 

116.99 

0.46 

21.88 

2.24 

4.43 

0.6 

33.49 

AL 

KQ 

4th 

7.69 

0 

0 

0 

82.96 

90.4 

0.92 

15.59 

1.78 

23.35 

0.3 

59.79 

AL 

KQ 

5th 

40.8 

0 

0 

0 

53.68 

42.95 

1.25 

18.29 

1.32 

11.48 

0.54 

47.1 

AL 

m 

1st 

9.64 

0 

0 

0 

79.13 

266.73 

0.3 

14.81 

2.34 

20.67 

0.36 

53.78 

AL 

KQ 

2nd 

38.92 

0 

0 

0 

94. S 

0.59 

23.11 

1.95 

4.32 

0.54 

39.34 

AL 

KQ 

3rd 

43.15 

0 

0 

0 

52.77 

71.99 

0.73 

14.76 

1.47 

14.91 

0.53 

38.17 

AL 

KQ 

4th 

30.51 

0 

0 

7.38 

59.18 

149.62 

0.44 

12.24 

1.63 

19.41 

0.55 

45.64 

AL 

KQ 

5th 

26.65 

0 

0 

16.05 

57.3 

99.14 

0.74 

13.73 

1.6 

18.14 

0.57 

60.99 

AL 

KQ 

1st 

22.63 

0 

0 

0 

77.37 

139.56 

0.55 

13.01 

1.39 

21.94 

0.35 

AL 

KQ 

2nd 

48.48 

0 

0 

0 

51.52 

74.1 

0.7 

16.13 

1.55 

9.74 

0.5 

AL 

KQ 

3rd 

53.64 

0 

0 

0 

44.1 

45.2 

0.98 

13.58 

1.56 

15.22 

0.52 

36.86 

AL 

[KQ 

4th 

30.06 

0 

0 

4.79 

65.15 

83.81 

0.83 

16.19 

1.54 

16.93 

0.48 

42.41 

AL 

HQ 

5th 

43.55 

0 

o" 

7.72 

24.06 

18.95 

1.68 

19.71 

1.31 

0.32 

0.71 

64.96 

AL 

IKQ 

1st 

7.95 

0 

0 

0 

76.77 

IBS 

14.04 

2.37 

17.13 

0.38 

74.44 

AL 

KQ 

2nd 

3.83 

0 

0 

0 

84.78 

283.63 

0.3 

9.97 

1.89 

25.15 

0.27 

60.77 

AL 

IKB 

3rd 

0 

0 

0 

0 

97.15 

329.64 

0.29 

4.17 

1.52 

32.4 

0.06 

54.75 

AL 

|KQ 

4th 

18.44 

0 

0 

0 

69.27 

116.38 

0.6 

14.4 

1.67 

19.02 

0.48 

69.55 

IKB 

5th 

11.28 

0 

0 

0 

86.55 

579.08 

0.15 

10.38 

2.64 

21.41 

0.24 

AL 

1st 

53.21 

0 

0 

0 

39.84 

38.69 

1.03 

17.62 

1.47 

6 

0.56 

59.37 

AL 

P 

2nd 

36.64 

0 

9.82 

0 

25.79 

47.63 

0.75 

18.39 

1.47 

1.3 

AL 

P 

3rd 

83.13 

0 

0 

0 

12.02 

9.16 

1.31 

10.51 

1.34 

0 

0.29 

77.42 

■ 

4th 

54.5 

0 

0.98 

0 

41.51 

58.01 

0.73 

15.38 

1.66 

6.08 

0.53 

40.14 

AL 

H 

0 

1.3 

0 

42.93 

151.79 

0.29 

12.14 

1.92 

15.64 

0.63 

AL 

\m 

1st 

92.98 

0 

0 

0 

4.81 

11.07 

0.43 

4.08 

1.37 

0 

0.13 

58.42 

AL 

m 

2nd 

61.64 

0 

IH 

32.23 

31.29 

1.17 

15.01 

0.51 

63.43 

Table  C3.  (cont’d)  Value  of  each  of  the  12  metrics  for  each  study  site  for  Extent  3. 
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State  1 

m\ 

Partial 

%  Ag 

%  Fw 

%  Dec 

%  Ev 

%  Mx 

MPS 

PD 

HI 

SIDI 

tm\ 

HI 

m 

Rt 

(ha) 

(no/ 

(m/ha) 

(%) 

(%) 

■ 

100  ha) 

AL 

m 

5th 

90.03 

0 

0 

0 

6.33 

6.3 

1 

7.67 

1.41 

0 

0.18 

71.83 

AL 

is 

1st 

72.08 

0 

0 

2.1 

22.69 

20.75 

1.19 

15.68 

1.44 

4.56 

0.43 

51.66 

AL 

gg 

2nd 

78.59 

0 

0 

9.84 

11.56 

16.19 

1.32 

14.26 

1.45 

0 

0.36 

72.43 

AL 

E9 

3rd 

81.45 

0 

0 

6.81 

9.98 

28.11 

0.6 

11.55 

1.37 

0 

0.32 

52.51 

AL 

m 

4th 

80.24 

0 

0 

0 

19.76 

Qggg 

1.47 

14 

1.35 

0 

0.32 

AL 

E9 

5th 

59.45 

0 

0 

0 

34.25 

23.22 

1.48 

22.29 

1.36 

0.69 

48.7 

AL 

gg 

1st 

60.74 

0 

0 

0 

22.36 

13.47 

1.66 

1.57 

0 

0.57 

47.8 

AL 

2nd 

72.71 

0 

0 

0 

26.4 

14.63 

1.8 

19.42 

1.5 

0 

0.4 

33.51 

AL 

gg 

3rd 

54.13 

0 

0 

3.88 

40.75 

36.49 

1.22 

17.2 

1.56 

9.43 

0.54 

41.35 

AL 

4th 

72.52 

0 

0 

0.6 

23.21 

12.2 

1.95 

19.15 

■gg 

0 

0.42 

43.59 

AL 

5th 

79.34 

0 

0 

0 

19.09 

13.88 

1.38 

14.59 

1.57 

0 

0.33 

33.68 

AL 

gg 

1st 

47.73 

0 

0 

o" 

51.21 

59 

0.87 

19.46 

1.66 

7.86 

0.51 

19.09 

AL 

ES 

2nd 

38.76 

0 

0 

16.62 

44.61 

70.47 

0.87 

8.14 

1.51 

26.11 

0.62 

91.73 

AL 

E3 

3rd 

18.79 

0 

0 

26.07 

55.13 

563.01 

0.14 

9.13 

2.61 

22.36 

0.59 

75.16 

AL 

ES 

4th 

24.44 

0 

0 

0.43 

75.13 

174.64 

0.43 

5.74 

1.72 

27.84 

0.38 

32.35 

AL 

ES 

5th 

13.78 

0 

0 

6.21 

80.01 

298.8 

0.29 

5.71 

1.77 

32.98 

0.34 

63.07 

AL 

1st 

22.87 

0 

0 

77.13 

0 

265.66 

0.29 

11.07 

2.58 

24.18 

0.35 

AL 

g 

2nd 

25.49 

0 

0 

70.66 

0 

161.74 

0.44 

10.46 

2.17 

22.75 

0.43 

41.15 

AL 

gg 

3rd 

4.32 

0 

0 

95.3 

0 

657.53 

0.14 

3.69 

2.28 

31.99 

0.09 

39.69 

AL 

4th 

11.46 

0 

0 

66.68 

285.27 

0.29 

13.09 

2.93 

20.42 

0.51 

60.1 

AL 

gg 

5th 

18.44 

0 

1.09 

13.86 

61.33 

263.47 

0.29 

15.53 

2.15 

14.84 

0.57 

59.15 

AL 

gg 

1st 

15.45 

0 

0 

0 

77.89 

529.68 

0.15 

13.89 

2.93 

15.27 

0.37 

45.38 

AL 

2nd 

2.74 

0 

0 

0 

95.15 

651.63 

0.15 

4.05 

2.25 

29.3 

0.09 

69.2 

AL 

gj 

3rd 

0 

0 

0 

0 

100 

684.62 

0.15 

0 

1.97 

36.23 

0 

AL 

gg 

4th 

14.63 

0 

0 

0 

85.08 

581.62 

0.15 

2.31 

1.96 

34.41 

0.25 

43.37 

AL 

5th 

69.38 

0 

0 

1.15 

7.22 

9.48 

0.88 

6.96 

1.58 

0 

0.48 

62.89 

AL 

m 

1st 

14.96 

0 

0 

3.01 

82.04 

290.63 

0.29 

7.63 

1.89 

27.46 

0.3 

44.13 

AL 

m 

2nd 

5.84 

0 

0 

0 

94.16 

641.96 

0.15 

4.23 

2.27 

31.4 

0.11 

— 

AL 

m 

3rd 

23.73 

0 

o" 

0 

76.27 

173.8 

0.44 

12.97 

1.76 

19.97 

0.36 

AL 

m 

4th 

44.14 

0 

0 

1.08 

53.79 

93.33 

0.59 

17.3 

2.1 

13.05 

0.52 

33.45 

AL 

IES 

40.62 

0 

0 

6.69 

52.69 

101.38 

0.59 

12.8 

1.57 

14.2 

0.55 

56.49 

AL 

ies 

1st 

20.15 

0 

0 

0 

79 

178.98 

0.44 

12.39 

1.82 

24.77 

0.34 

44.41 

AL 

IES 

2nd 

47.09 

0 

0 

0 

50.68 

85.94 

0.59 

17.18 

1.79 

11.74 

0.52 

36.28 

AL 

IES 

3rd 

1.73 

0 

0 

4.93 

90.82 

650.37 

0.15 

4.08 

2.16 

32.57 

0.17 

53.04 

AL 

IES 

4th 

6.75 

0 

0 

0.72 

0.15 

4.03 

2.09 

32.39 

0.14 

57.91 

AL 

IES 

5th 

83.7 

0 

0 

0.42 

14.83 

20.21 

1  Hi 

8.38 

1.38 

0 

0.28 

38.19 

AL 

IQ 

1st 

0 

14.32 

0 

0.26 

0 

43.4 

0.34 

6.01 

4.34 

0 

0.82 

74.21 

AL 

IES 

2nd 

0 

18.51 

0 

8.55 

0 

83.1 

0.33 

9.33 

0.11 

0.84 

79.06 

AL 

IES 

3rd 

0 

0 

0 

0 

48.24 

66.21 

0.73 

14.09 

1.47 

19.44 

0.69 

70.04 

AL 

IES 

4th 

0 

18.34 

0 

0 

55.01 

507.97 

0.14 

10.9 

2.43 

22.04 

0.63 

71.41 

AL 

IES 

5th 

0 

0 

0 

5.24 

23.82 

100.34 

0.29 

9.53 

1.64 

13.92 

0.69 

67.2 

AL 

IES 

1st 

8.37 

0 

0 

0 

73.92 

251.92 

0.29 

12.97 

IBS 

0.43 

58.41 

Table  C3.  (cont’d)  Value  of  each  of  the  12  metrics  for  each  study  site  for  Extent  3. 
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State 


AL¬ 


AI- 

AT 

al" 

JT 

AL" 

A T 
A r 
aT 
a r 
aT 
al" 
al" 
al" 
al" 
aT 
al" 
al" 
al" 
al" 
al” 
al” 
al” 
aT 
al“ 
al" 
7T” 
al" 
al” 
7T“ 
al” 
al" 
al" 
al" 
aT 
al" 
al” 
71“ 
al" 
fl” 
fl“ 
fl” 


Rt 


43 


43 


43 


43 


44 


44 


44 


44 


44 


Partial 


Rt 


2nd 


3rd 


4th 


5th 


1st 


2nd 


3rd 


4th 


5th 


1st 


2nd 


3rd 


4th 


5th 


1st 


2nd 


3rd 


4th 


5th 


1st 


2nd 


3rd 


4th 


5th 


1st 


2nd 


3rd 


4th 


5th 


1st 


2nd 


3rd 


4th 


5th 


1st 


2nd 


3rd 


4th 


5th 


1st 


2nd 


1  3rd 


%  Ag 


39.62 


22.38 


0 


4.45 


30.86 


31.72 


65.88 


65.01 


36.07 


37.09 


69.44 


69.96 


70.72 


64.95 


49.28 


52.24 


68.99 


6.04 


9.56 


59.41 


63.6 


73.75 


74.79 


47.02 


65.58 


52.37 


42.38 


20.13 


21.19 


73.28 


79.34 


91.28 


95.72 


61.95 


27.59 


37.58 


2.27 


%  Fw 


0 


11.17 


0 


0 


13.82 


%  Dec 


0.87 


19.5 


0.32 


0 


12.58 


0.06 


15.71 


24.13 


0.83 


2.79 


7.25 


32.28 


0 


%  Ev 


1.82 


0.49 


2.57 


1.53 


2.66 


6.59 


4.67 


6.32 


8.3 


1.47 


0.03 


0 


65.38 


59.45 


76.79 


%  Mx 


60.38 


76.86 


98.88 


97.01 


8.03 


3.01 


2.28 


0.95 


48.63 


63.18 


30.1 


28.21 


60.08 


62.07 


27.77 


20.59 


29.28 


23.39 


50.04 


39.6 


9.07 


0 


37.73 


36.08 


22.75 


24.76 


37.99 


25.32 


30.38 


21.85 


9.08 


46.37 


24.25 


16.74 


0 


4.25 


3.88 


2.95 


4.41 


MPS 


(ha) 


82.67 


103.96 


332.66 


657.67 


27.13 


20.12 


15.34 


2.16 


65.96 


88.28 


34.92 


22.15 


106.87 


70.88 


47.57 


12.77 


20.08 


18.97 


56.98 


42.15 


18.13 


48.28 


127.16 


124.2 


22.17 


27.9 


128.95 


27.29 


81.27 


31.97 


56.13 


267.06 


34.21 


44.64 


19.93 


14.61 


220.67 


118.71 


85.73 


647.26 


PD 


(no/ 


100  ha) 


0.73 


0.74 


0.3 


0.15 


0.3 


0.15 


0.15 


0.44 


0.74 


0.74 


0.88 


1.31 


0.59 


0.88 


0.58 


1.61 


1.46 


1.31 


0.88 


0.5 


0.64 


0.3 


0.29 


1.03 


0.89 


0.29 


1.17 


0.59 


0.88 


0.44 


0.3 


0.73 


0.44 


0.44 


0.29 


0.15 


0.58 


0.73 


0.15 


ED 


(m/ha) 


15.71 


12.32 


3.26 


2.82 


6.55 


2.86 


1.82 


1.22 


14.05 


17.35 


18.01 


16.06 


15.46 


12.23 


13.93 


14.1 


16.92 


15.29 


20.26 


17.39 


4.96 


13.18 


11.38 


11.08 


12.19 


12.47 


10.16 


18.68 


11.16 


10.07 


7.89 


6.13 


11.44 


7.64 


5.51 


2.47 


5.37 


16.2 


16.58 


4.21 


MSI 


1.85 


1.53 


1.53 


2.17 


1.48 


1.33 


1.13 


1.62 


1.6 


1.57 


1.89 


1.6 


1.6 


1.76 


2.12 


1.37 


1.62 


1.35 


1.67 


1.47 


1.32 


1.49 


1.8 


1.85 


1.38 


1.68 


1.72 


1.52 


1.33 


1.48 


1.45 


1.58 


1.65 


1.22 


1.33 


1.21 


1.81 


1.59 


1.68 


2.21 


TCAI 


(%) 


14.26 


26.05 


34.34 


32.45 


0 


17.53 


13.37 


0.02 


0.66 


10.72 


20.6 


3.37 


4.5 


0.47 


6.07 


7.33 


2.45 


0.09 


0.54 


0.09 


7.1 


0.1 


14.04 


8.43 


5.02 


27.45 


5.48 


7.37 


0 


11.37 


11.66 


14.37 


30.91 


SIDI 


0.48 


0.36 


0.02 


0.06 


0.57 


0.57 


0.43 


0.34 


0.61 


0.65 


0.5 


0.47 


0.5 


0.51 


0.48 


0.44 


0.46 


0.41 


0.52 


0.51 


0.57 


0.49 


0.78 


0.69 


0.5 


0.47 


0.4 


0.38 


0.63 


0.5 


0.62 


0.74 


0.71 


0.67 


0.4 


0.34 


0.16! 


0.08 


0.51 


0.49 


0.5 


0.39 


IJI 


(%) 


41.12 


26.41 


40.02 


68.62 


71.54 


68.98 


55.66 


63.89 


74.78 


51.36 


31.42 


49.11 


39.18 


31.88 


37.68 


63.26 


43.82 


22.75 


43.84 


63.4 


70.95 


76.82: 


53.75! 


19.86 


71.15 


22.29 


65.82 


41.37 


74.61 


63.81 


82.76 


80.04 


42.98 


52.9 


51.08! 


24.88 


82.72 


56.27 


32.66 


50.63 


Table  C3.  (cont’d)  Value  of  each  of  the  12  metrics  for  each  study  site  for  Extent  3. 
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State  1 

E9 

Partial 

%  Ag 

%  Fw 

%  Dec 

%  Ev 

%Mx 

MPS 

PD 

ED 

2SH 

SIDI 

IJI  | 

m 

Rt 

(ha) 

(no/ 

(m/ha) 

(%) 

(%) 

m 

100  ha) 

FL 

1 

4th 

3.38 

9.23 

0 

55.76 

12.16 

128.96 

0.6 

16.79 

1.88 

17.54 

0.65 

72.57 

FL 

“7 

5th 

0 

0 

37.01 

6.74 

33.23 

1.32 

23.7 

2.1 

1.56 

0.73 

60.75 

FL 

m 

1st 

8.13 

0 

61.97 

0 

67.24 

1.04 

10.78 

1.7 

24.23 

0.58 

63.42 

FL 

m 

2nd 

0 

0 

0 

32.56 

0 

22.91 

1.42 

24.59 

1.61 

0 

0./6 

74.35 

FL 

~ 

3rd 

0 

0 

0 

19.06 

0 

27.91 

0.68 

14.76 

2.12 

0 

0.82 

85.61 

FL 

~ 

4th 

0 

0 

0 

17.41 

0 

110.39 

0.16 

7.82 

1.77 

1.55 

0.71 

72.73 

FL 

~2 

5th 

0 

10.48 

0 

0.22 

mm 

0.86 

9.92 

1.47 

0 

0.75 

63.57 

FL 

m 

1st 

0 

0 

o" 

66.96 

12.95 

279.14 

0.29 

8.54 

1.74 

27.67 

0.52 

65.91 

FL 

m 

2nd 

0 

0 

0 

88.72 

4.76 

653.79 

0.14 

4.23 

2.22 

29.32 

0.21 

66.45 

FL 

m 

3rd 

0 

6.38 

0 

83.99 

0 

632.88 

0.14 

5.67 

2.32 

26.37 

0.28 

55.51 

FL 

m 

4th 

0 

0 

0 

74.2 

0.02 

129.68 

0.57 

mm 

1.69 

26.89 

0.43 

65.21 

FL 

m 

5th 

21.21 

2.67 

o" 

52.08 

mm 

83.37 

0.71 

15.99 

1.72 

18.03 

0.67 

69.79 

FL 

4 

1st 

40.92 

6.25 

0 

42.47 

8.65 

65.65 

0.87 

19.64 

1.5 

9.39 

0.64 

63.93 

FL 

4 

2nd 

44.42 

EES 

0 

42.85 

7.92 

62.2 

0.88 

25.36 

2.13 

B8! 

0.61 

50.89 

FL 

4 

3rd 

42.46 

2.63 

0 

53.1 

0 

126.97 

0.44 

14.64 

1.7 

8.45 

0.54 

46.02 

FL 

4 

4th 

9.26 

10.11 

0 

75.35 

0 

292.21 

0.29 

11.58 

2.02 

17.07 
— — 

0.41 

64.18 

FL 

4 

5th 

8.12 

4.81 

o" 

87.07 

0 

314.49 

0.29 

4.01 

1.94 

32.37 

0.23 

62.89 

FL 

m 

1st 

2.13 

10.36 

0 

87.51 

0 

663.54 

0.15 

1.61 

2.03 

33.96 

FL 

m 

2nd 

WBk 

14.54 

0 

80.76 

0 

648.61 

0.15 

2.19 

2.09 

33.59 

0.32 

35.47 

FL 

m 

3rd 

24.55 

2.23 

0 

67.97 

4.21 

101.52 

0.73 

12.78 

2.05 

26.16 

0.48 

56.25 

FL 

~ 

4th 

18.16 

0 

71.06 

2.25 

522.44 

0.15 

16.82 

3.18 

9.79 

0.46 

53.56 

FL 

~ 

5th 

2.52 

0 

52.36 

9.69 

245.09 

0.29 

10.82 

1.91 

18.22 

0.68 

80.71 

FL 

~ 

1st 

36.79 

2.05 

0 

51.66 

1.57 

63.2 

0.87 

27.53 

1.87 

2.28 

0.59 

45.31 

FL 

~ 

2nd 

29.01 

0 

0 

66.47 

3.45 

95.54 

0.73 

1.93 

17.02 

0.47 

39.58 

FL 

m 

3rd 

41.76 

0 

0 

50.48 

0.65 

67.45 

0.76 

24.55 

1.83 

0.21 

0.57 

36.86 

FL 

m 

4th 

12.06 

0 

0 

50.88 

0 

86.9 

0.59 

18.9 

1.98 

1.59 

0.65 

73.52 

FL 

m 

5th 

7.87 

0.67 

0 

89.97 

0.19 

308.49 

0.29 

8.51 

1.79 

25.52 

0.18 

47.65 

FL 

7 

1st 

0.05 

17.25 

0 

74.96 

5.81 

676.07 

0.14 

2.47 

2.11 

33.17 

0.4 

34.93 

FL 

~ 

2nd 

2.37 

0 

0 

57.43 

26.98 

292.34 

0.29 

12.86 

2.46 

21.52 

0.59 

62.57 

FL 

m 

IBSEI 

0 

0.01 

3.19 

46.26 

46.19 

329.12 

0.29 

5.22 

1.91 

30.98 

0.57 

40.04 

FL 

m 

4th 

22.27 

12.65 

0 

1.84 

59.84 

100.29 

0.74 

16.07 

1.57 

20.78 

0.58 

54.26 

FL 

~ 

5th 

29.41 

8.87 

0 

11.48 

24.41 

1.76 

8.74 

0.66 

55.34 

FL 

"8 

1st 

56.93 

2.94 

0 

15.97 

22.89 

40.19 

1.04 

28.48 

1.87 

0.11 

0.6 

51.71 

FL 

2nd 

32.59 

7.12 

0 

25.89 

25.61 

65.86 

0.89 

27.47 

1.74 

3.15 

0.75 

57.22 

■ 

3rd 

0 

0.41 

0 

80.93 

14.18 

650.07 

0.15 

2.44 

2.06 

33.68 

0.32 

37.99 

FL 

■ 

4th 

0 

15.24 

0 

77.86 

IBS] 

684.49 

0.15 

0 

1.97 

36.14 

0.37 

86.99 

FL 

n 

5th 

0.08 

0 

0 

96.18 

3.75 

684.44 

0.15 

0.33 

1.97 

36.12 

0.07 

e 

1st 

0 

0 

0 

77.15 

0 

0.15 

18.56 

2.77 

14.65 

0.38 

61.25 

FL 

m 

2nd 

0 

0 

0 

61.54 

0 

132.81 

0.46 

19.44 

1.69 

12.15 

0.58 

71.18 

FL 

E 

0 

0 

82.9 

0 

0.15 

13.31 

2.24 

22.63 

0.29 

29.63 

FL 

E 

4th 

0 

0 

0 

75.59 

8.2 

551.32 

0.15 

13.91 

2.24 

22.71 

0.41 

65.33 

FL 

O 

5th 

0 

1.86 

0 

84.38 

10.26 

0.15 

5.84 

2.41 

25.78 

0.28 

53.72 

Table  C3.  (conf  d)  Value  of  each  of  the  12  metrics  for  each  study  site  for  Extent  3 
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Table  C3.  (cont’d)  Value  of  each  of  the  12  metrics  for  each  study  site  for  Extent  3. 


Table  C4.  Value  of  each  of  the  12  metrics  for  each  study  site  for  Extent  4. 
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Table  C4.  (cont’d)  Value  of  each  of  the  12  metrics  for  each  study  site  for  Extent  4. 
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